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ABSTRACT 


A  computer  analysis  has  been  performed  to  evaluate  the  energy  dependent  re¬ 
sponse  of  a  capture  tank  to  the  gamma-ray  cascades  emitted  from  excited  The 

GAMINT  code  was  developed  to  simulate  the  decay  of  the  '^®U  nucleus,  formed  in 
the  ^^®U(n,7  reaction,  in  order  to  provide  the  source  spectrum  for  the  complete 
analysis  of  the  capture  tank  efficiency.  This  model  determines  the  energies  of  the 
gamma-ray  cascades,  the  order  of  emission  of  the  gamma  rays  in  a  cascade,  and 
the  gamma-  ray  multiplicities  by  Monte  Carlo  techniques.  A  gamma-ray  emission 
spectrum  for  the  excited  nucleus  is  generated. 

In  the  GAMINT  code,  known  level  data  for  is  used  below  1  MeV.  A 
statistical  approach  based  on  the  back-shifted  Fermi  gas  model  is  used  for  the  con¬ 
tinuum  level  density.  A  single- particle  model  description  for  transition  rates,  with 
hindrance  factors  applied,  is  used  to  determine  the  gamma  ray  transition  proba¬ 
bilities.  Internal  conversion  probabilities  are  determined  and  the  inclusion  of  this 
competing  process  suppresses  the  low  energy  portion  of  the  gamma  spectrum. 

A  capture  tank  responds  to  the  combined  effect  of  the  gamma  rays  of  various 
energies  from  a  cascade,  after  being  transported  through  the  matenal  between  the 
sample  and  the  capture  tank.  Examined  is  the  energy  deposition,  in  a  capture  tank, 
by  the  cascades  generated  from  resonant  and  off-resonant  capture  in  a  sample. 
Internal  conversion  has  a  negligible  effect  on  the  average  cascade  energy  deposited 
in  the  tank.  Off- resonant  (volumetric)  capture  deposits,  on  the  average,  less  energy 
than  resonant  (surface)  capture  in  the  capture  tank  as  a  result  of  self- shielding  of 
the  gamma  rays  in  the  capture  sample. 
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CHAPTER  1 

INTRODUCTION  AND  HISTORICAL  REVIEW 


1.1  Introduction 

Gamma  rays  are  produced  m  many  nuclear  reactions.  In  neutron  capture 
reactions,  a  nucleus  becomes  excited  by  the  addition  of  the  captured  neutron  and 
then  decays  to  the  ground  state  by  the  emission  of  sequential  gamma  rays.  This 
series  of  gamma  rays  is  called  a  gamma  cascade.  In  this  work,  the  computer  code 
GAMINT  is  developed  to  simulate  the  gamma  cascades  formed  in  the  ■“U(n,7)^^®U 
reaction. 

The  study  of  the  sequential  emission  of  gamma  rays  may  be  performed  ex¬ 
perimentally.  The  current  generation  of  expenments  resolve  the  energies  of  one  or 
two  of  the  gamma  rays  in  a  cascade.  The  determination  of  the  angular  correlation 
between  a  gamma  emitted  in  some  preferred  direction  and  the  subsequent  gamma 
ray  is  also  bmited  to  two  gammas.  Future  experiments  Ti  will  be  able  to  simul¬ 
taneously  measure  the  energies  of  up  to  a  half  dozen  gamma  rays  produced  in  the 
step  by  step  deexcitation  of  a  single  excited  nucleus.  This  ability  to  measure  the 
multifold  coincidence  of  gamma  rays  wiD  not  be  available  until  1994  and  will  not  be 
without  tremendous  effort  and  expense. 

Internal  conversion  is  an  additional  decay  process  which  competes  with  gamma 
decay.  For  low  energy  transitions  in  nuclei  of  high  mass  number,  the  internal  conver¬ 
sion  process  dominates  over  the  gamma  emission  process.  It  is  therefore  necessary 
to  include  internal  conversion  in  the  cascade  model.  This  work,  for  the  first  time, 
includes  internal  conversion  in  a  decay  model. 

Developing  a  computer  code  to  simulate  the  decay  of  the  excited  nucleus 
allows  one  to  study  the  character  of  the  gammas  emitted  in  the  cascade,  without  the 
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expense  of  repeated  experimental  procedure.  The  computer  model  determines  the 
energy,  the  multipolarity,  order  of  emission  of  a  gamma  ray  in  a  cascade,  and  the 
number  of  gamma  rays  emitted  in  a  cascade.  The  angular  distribution  between  the 
first  gamma  ray  in  a  cascade  and  any  subsequent  gamma  ray  in  the  same  cascade  is 
also  determined.  In  addition,  the  computer  code  generates  the  source  spectrum  for 
the  ■^U(n.7)^^®U  reaction.  This  spectrum  is  used  in  this  work  for  the  determina¬ 
tion  of  the  energy  dependent  efficiency  for  gamma-ray  detection  in  neutron  capture 
tank  experiments.  This  is  important  so  that  the  data  collection  efforts  of  future 
experiments  may  be  improved  and  the  results  of  past  experiments  may  be  verified. 

The  ^^*U(n,7)^^®U  reaction  is  of  importance  in  the  design  of  fast  reactors, 
and  therefore  must  be  fully  understood  in  order  to  optimize  fuel  development  in 
future  reactor  programs  utihzing  fast  reactor  technology.  Experiments  have  been 
performed  at  Rensselaer  Polytechnic  Institute  2:  and  recently  elsewhere  3]  4!  to 
study  the  ■^U(n,7)^^®U  reaction,  using  a  large  capture  tank.  An  analysis  of  the 
data  collected  in  the  Rensselaer  experiments  '5i  by  Harris  et  al.  6'  examined  the 
capture  tank  efficiencies  for  average  s-wave  and  p-wave  neutron  capture.  However, 
the  capture  tank  responds  to  the  combined  effect  of  the  gamma  rays  of  various 
energies  from  a  cascade.  Therefore,  the  cascade  structure  has  a  profound  effect  on 
the  response  of  the  capture  tank  to  individual  cascades.  Further  analysis  is  necessary" 
to  address  the  efficiency  of  the  capture  tank  for  detecting  gamma  rays  emitted  in 
neutron  capture  experiments.  This  work,  for  the  first  time,  analyzes  the  multiple 
gammas  in  the  individual  cascades  of  the  ■“U(n,7)^^®U  reaction,  and  the  resulting 
energy  deposition  in  the  capture  trink. 

Previous  work  in  gamma  cascade  modeling  has  been  performed  by  Booth  et 
al.  7!,  Takahashi  [8],  and  Yamamuro  et  ed.  ^9].  Experimental  results  are  found  in 
the  work  of  Campion  et  al.  10]  and  John  et  al.  111.  This  chapter  summarizes  the 
hterature  applicable  to  the  gamma  cascade  model  developed  in  this  work. 


1.2 


Objectives 


The  objectives  of  this  thesis  are: 

•  To  develop  a  computer  simulation  of  the  gamma-ray  cascade  process  in  the 
■^*U(n,7)^^U  reaction  in  order  to  provide  the  source  cascades  in  a  capture 
tank  model. 

•  To  transport  the  cascade  gamma  rays  from  the  source  to  the  capture  tank, 
through  the  intervening  material,  using  a  point-kernel  calculation  code. 

•  To  examine  the  energy  deposition  in  the  capture  tank  for  the  following  specific 
cases: 

1.  Resonant  (surface)  capture  in  the  sample. 

2.  Off- resonant  (volume)  capture  in  the  sample. 

3.  Cascades  in  which  internal  conversion  is  considered. 

1.  Cascades  in  which  internal  conversion  is  not  considered. 

1.3  Calculations  of  Booth  et  al. 

Booth  et.al  7l  used  a  calculational  technique  based  on  the  works  of  Yost  et 
al.  [12]  T3],  in  which  the  capture  gamma-ray  yields  for  neutron  capture  in  ‘^U 
were  calculated  by  a  technique  which  incorporates  existing  experimental  results. 
The  emphasis  of  this  work  was  on  implementation  of  the  Ccilculation  on  ’“U  and 
comparison  on  an  absolute  basis,  the  resTilts  with  differential  measurements  of  cap¬ 
ture  gamma-ray  yields.  .\n  advantage  of  this  calculation  is  that  gamma-ray  group 
boundaries  may  be  altered  and  new  results  determined  at  a  fraction  of  the  cost 
and  time  associated  with  the  experimental  work.  The  calculational  technique  uses 
several  nuclear  models  with  the  constraint  of  conservation  of  binding  energy.  .\lso 
employed  is  a  simple  integral  gamma-ray  yield  measurement  to  obtain  differential 
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yield  data.  These  experiments  were  conducted  at  the  Tower  Shielding  Facility  at 
Oak  Ridge  National  Laboratory.  These  experiments  were  used  to  fix  certain  pa¬ 
rameters  by  requiring  the  model  to  predict  the  experimental  results  on  an  absolute 
basis. 

In  the  generation  of  the  neutron-energy-dependent  capture  gamma-ra\'  yields. 
Booth  assumed  that  the  composite  epithermal  spectrum  was  dominated  by  the  cap¬ 
ture  of  s,  p,  and  d-wave  neutrons.  The  calculation  of  the  gamma-ray  spectrum 
corresponding  to  each  of  the  possible  capture  states  was  accomplished  by  using  the 
DUCAL  code  of  Yost.  In  this  code,  the  gamma-ray  cascade  process  is  described 
in  terms  of  (1)  parameters  which  define  relative  electric  and  magnetic  multipole 
transition  strengths  and  which  incorporate  the  appropriate  selection  rules,  (2)  sta¬ 
tistical  assumptions  regarding  energy-level  density  and  spin  and  parity  distributions 
in  the  energy  continuum,  (3)  known  energy  levels  in  the  lower  energy  region,  and 
(4)  energy  levels,  gamma-ray  transition  probabilities,  and/or  estimates  of  spin  and 
parity  distributions  generated  with  an  appropriate  nuclear  model  over  the  range  of 
energies  for  which  the  nuclear  model  has  been  shown  to  agree  better  than  the  model 
using  statistical  assumptions. 

The  DUCAL  code  of  Yost  incorporates  some  parameters  which  are  propor¬ 
tional  to  the  matrix  elements  of  the  radiative  transition  process.  In  generating  the 
spectrum  calculations.  Booth  relied  on  the  accuracy  of  the  thenncil- cap  lure  primary 
line  intensities  found  in  the  literature  in  order  to  fit  the  parameters  of  the  model. 
The  measured  thermal  line  intensities  were  reproduced  within  10  percent.  Based  on 
this  fact,  Booth  applied  the  parameters  which  described  the  thermal  spectrum  to 
compound  nucleus  capture  states  at  higher  energies. 

The  statistical  description  of  the  nuclear  energy  levels  in  the  continuum  region 
used  in  the  DUCAL  code  is  based  on  the  Gilbert  and  Cameron  !24l  formulation. 
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which  predicts  an  energy-dependent  observable  level  spacing,  D(U)  given  by: 


D(U)  = 


Texp[-{U  -  EJ/T’,  E<E, 
12v^  a^t^exp(  —  \/aE),  E  >  E* 


where  U  =  E  -  A  and  E,  =  2.5  —  150/.\  —  the  pairing  energy  in  Me\’.  The  parameter 
A  is  the  fictive  ground  state  level.  In  Booth  s  cadculations.  the  parameters  relevant 
to  the  Gilbert  and  Cameron  formulation  shown  above  are  determined  in  part  from 
the  resonance  data  as  given  in  BN'L-325  14  as  well  as  the  known  discrete 
energy  levels. 

The  DUC.A.L  code  requires  the  energy,  spin,  and  parity  of  the  energy  levels  in 
order  to  determine  the  gamma-ray  spectrum  for  a  given  capture  state.  The  required 
properties  are  only  known  for  the  energy  levels  of  up  to  about  1  MeV.  A 
level  scheme  was  compiled  bcised  on  existing  experimental  data  and  the  results  of 
a  nuclear  model  calculation  suitable  for  deformed  odd-A  nuclei.  In  this  composite 
level  scheme,  level  properties  from  measurements  and  detailed  in  the  Nuclear  Data 
Tables  [15]  were  used  in  the  energy  range  0.0  <  E  <  0.83  MeV,  and  level  properties 
from  the  nuclear  model  were  used  in  the  range  0.83  <  E  <  2.0  MeV.  Individual  levels 
in  the  range  1.23  <  E  <  2.0  MeV  were  weighted  by  the  ratio  of  the  density  of  levels 
as  suggested  by  the  Gilbert  and  Cameron  formulation  to  that  of  the  nuclear  model. 
Thus  the  nuclear  model  was  used  only  to  predict  the  spin  and  parity  distributions 
in  the  1.23  <  E  <  2.0  MeV  range.  The  primary  line  intensities  used  in  Booth’s 
calculation  were  taken  from  BNL-325  T4I. 

In  Booth’s  work,  comparison  of  the  calculated  energy  spectrum  to  experimen¬ 
tally  determined  spectrum  was  not  done  in  the  region  below  3.0  MeV,  because  of 
relatively  large  errors  in  the  experimental  results.  The  measured  discrete  photon 
intensity  at  4.06  MeV  was  about  11%  below  the  intensity  measured  in  experiments 
performed  by  Sheline  et  al.  :16i  as  stated  in  Booth  s  paper.  .4s  an  additional  check 
on  the  calculated  gamma-ray  vdelds,  the  percent  of  the  binding  energy,  which  must 


6 


be  conserved  in  the  calculation,  is  given.  This  figure  in  Booth's  calculation  for  the 
capture  of  thermal  neutrons  is  given  as  89.0%. 

1.4  Calculations  of  Takahashi 

In  Takahashi's  work  i8l.  the  gamma- ray  production  cross  sections  of '“U(n.7) 
reactions  are  evaluated  for  thermeil  and  epithermal  energies  primarily  using  a  sta¬ 
tistical  model  for  the  gamma  cascade  process.  Also  in  this  work,  the  gamma- rav 
spectra  is  determined  by  using  a  model  which  takes  into  account  the  known  discrete 
energy  levels,  and  uses  a  statistical  theory  for  the  continuum  energy  levels.  This  is 
done  using  a  modified  version  of  the  DUCAL  code  of  Yost  called  DUCAL- AV'^T. 

In  the  process  of  calculating  the  gamma  cascade  using  discrete  and  continuous 
energy  levels,  an  important  consideration  is  the  matching  of  the  discrete  level  density 
to  the  continuous  level  density  at  the  junction  energy  (Eg).  In  the  Takahashi  work, 
the  formulation  of  Gilbert  and  Cameron  is  used  to  provide  the  level  densities  in 
the  continuum.  Since  the  discrete  levels  formed  in  the  (n,7)  reaction  are  not  aU  of 
the  discrete  levels  described  by  Gilbert  and  Cameron’s  formula,  and  some  of  the 
discrete  levels  are  not  populated  in  the  (n,7)  reaction,  normahzation  of  the  discrete 
level  density  is  earned  out.  In  the  DUC.AL-.AVT  code,  this  normadization  is  carried 
out  in  such  a  way  that  the  level  density  of  the  discrete  levels  matches  the  level 
density  of  the  continuum  at  the  junction  energy. 

The  branching  ratios  for  the  transitions  to  and  from  the  levels  in  the  continuum 
region  are  determined  using  the  formulation  of  Blatt  and  Weisskopf  171  with  a 
hindrance  factor  used  to  adjust  the  calculated  ratio.  These  values  are  shown  in 
Table  1.1. 

In  Takahashi’s  work,  the  branching  ratios  are  assumed  to  be  independent  of 
the  initial  excitation  energy.  The  transition  probabilities  are  determined  by  a  code 
called  NGAMMA.  The  branching  ratios  do  not  depend  strongly  on  the  energy  of  the 


Table  1.1:  Branching  ratios  used  in  Takahashi^s  calculation 


El 

Ml 

E2 

M2 

Branching  Ratio 

8.0  X  10"^^ 

1.99  X  10"^ 

1.0  X  10-' 

9.9  X  10-^ 

Energy  (MeV) 

Transition  Probabilities  (sec 

0.1 

2.4 

4.8 

bihbb 

2.8  X  10'- 

initial  states;  therefore,  Takahashi  assumes  that  the  energy-independent  branching 
ratios  used  in  the  DUCAL-AVT  code  are  a  good  approximation. 

Since  the  intensities  of  the  primary'  transitions  in  the  low  energy  discrete  states 
were  available  to  Takahashi  in  the  Nuclear  Data  Tables  15',  the  transition  proba¬ 
bilities  for  primary  transitions  to  these  states  were  calculated  from  these  data.  The 
values  for  the  transition  probabilities  used  in  Takahashi’s  work  are  also  found  in  Ta¬ 
ble  1.1.  In  Takahashi’s  work,  he  found  that  if  he  calculated  the  transition  probability 
for  the  primary  transitions  in  the  discrete  region,  the  intensities  for  those  transitions 
were  about  four  times  smaller  than  the  experimentally  observed  intensities.  This  is 
a  result  of  the  fact  that  the  statistical  model  treats  the  probabilities  of  transition  to 
discrete  states  in  the  same  way  as  transition  probabilities  to  continuous  states,  and 
does  not  take  into  account  non-statistical  capture  processes.  Therefore,  Takahashi 
takes  the  transition  probabilities  to  the  discrete  levels  to  be  four  times  larger  than 
the  transition  probabilities  to  the  continuous  states. 

In  Takahashi’s  work,  his  calculated  spectrum  is  compared  to  the  experimen¬ 
tally  determined  spectrum  of  John  et  ai.  11’,.  The  spectrum  differs  from  the  ex¬ 
perimental  results  of  John  et  al.,  specifically  in  the  position  of  the  maximum.  The 
maximum  of  the  experimental  spectrum  lies  at  1.8  MeV,  while  that  of  the  calculated 
spectrum  hes  around  1.2  MeV.  Takahcishi  was  unable  to  shift  this  peak  by  adjusting 
the  parameters  in  the  Gilbert  and  Cameron  formulation. 
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Table  1.2:  Comparison  of  the  calculated  gamma  yields  from  neutron 
capture  in  *“U 


Energy  (MeV) 

Intensity  (gamma  rays  per  100  captures) 

Booth  et  al. 

Takahashi 

0.0  -  0.25 

- 

124.12 

0.25  -  0.50 

- 

13.47 

0.50  -  0.75 

- 

43.28 

0.75  -  1.00 

- 

28.12 

1.00  -  1.25 

29.9 

34.23 

1.25  -  1.50 

25.5 

33.20 

1.50  -  1.75 

25.7 

29.80 

1.75  -  2.00 

25.4 

25.63 

2.00  -  2.25 

21.3 

21.30 

2.25  -  2.50 

17.2 

17.07 

2.50  -  2.75 

13.8 

13.09 

2.75  -  3.00 

10.6 

9.49 

3.00  -  3.25 

5.4 

6.43 

3.25  -  3.50 

3.5 

4.06 

3.50  -  3.75 

5.3 

2.34 

3.75  -  4.00 

3.9 

1.014 

4.00  -  4.25 

10.8 

5.72 

4.25  -  4.50 

0.024 

0.0514 

4.50  -  4.75 

0.43 

0.203 

An  examination  of  the  conservation  of  binding  energy  in  the  experimental  and 
calculated  spectrum  is  accomphshed  by  analyzing  the  integral: 

I(E,)E,dE,.  (1.1) 

JlMeV 

The  experimental  data  covers  95%  of  the  binding  energy  (Eg)  and  Takahashi’s 
calculation  covers  80%  of  the  binding  energy,  for  this  integrated  region.  Table  1.2 
compares  the  calculated  neutron  capture  yields  for  the  work  of  Booth  et  al.  and 
Takahashi. 


1.5  Calculations  of  Yamamuro  et  al. 


In  Yamamuro’s  work  :9l,  the  capture  cross  sections  and  gamma- ray  spectra  for 
neutron  induced  reactions  in  ^®^Ta  and  ^^’^Au  were  calculated  using  a  low-lying  level 
fitting  method  in  which  the  constant  temperature  partial-level  density  formula  and 
the  spin-selected  level  cumulation  were  used.  The  calculation  of  the  capture  cross 
section  and  the  gamma-ray  spectra  uses  the  Brink- Axel  form  to  represent  the  giant 
resonance  of  photon  absorption  for  the  gamma-ray  profile  function. 

The  nuclear  temperature  is  normally  determined  by  fitting  the  cumulative 
low-lying  levels  with  the  function  describing  the  the  number  of  levels.  The  low-lying 
level  cumulation  for  heavy  nuclei  is  usually  below  the  expected  exponential  increase 
as  a  function  of  energy  as  a  result  of  levels  missing  from  the  experimental  data. 
Therefore,  the  temperature  cannot  be  uniquely  determined  from  the  cumulative 
plot  for  an  element.  More  than  one  temperature  may  therefore  be  used  in  the 
calculation  of  the  gamma-ray  spectrum  instead  of  the  constant  temperature  model. 

In  Yamamuro’s  work,  the  statistical  model  called  GNASH  was  used  for  the 
calculation  of  the  gzunma-ray  spectrum  for  ^®*Ta  and  *®^Au.  This  calculation  was 
performed  for  incident  neutron  energies  in  the  range  of  10  keV  to  4.5  MeV. 

Since  the  average  gamma-ray  transition  probabilities  depend  strongly  on  the 
nuclear  level  density,  the  level  density  p(E,J)  is  an  important  quantity  in  the  calcu¬ 
lation  of  the  gamma-ray  spectrum.  In  the  Yamamuro  paper,  the  Fermi-gas  model 
formula  is  used  and  is  referred  to  as  the  partial  level  density  formula,  which  is 
essentially  the  density  of  levels  at  energy  E  of  each  spin  J,  not  of  all  observable  spin. 


where 


C„  =  24's/2(0.146)-JaA. 
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and  the  spin  tenr  R(J,E)  is  defined  by  the  relationship 

R(J.E)=(2J-l)exp{^i^'}. 

where  c'’(E)  is  defined  as  the  spin  cutoff  factor,  and  is  expressed  by 

o-^(E)  =  0.146{a(E  - 

In  Yamamuro  s  work,  the  constant  temperature  formula  is  used  in  the  lower 
excitation  energy  region,  and  is  joined  to  the  Fermi  gas  level  density  formula  at 
energy  E*.  The  constant  temperature  formula  is: 

p(E,J)  =^(E.)exp{^^;^}R(J,E),  (1.3) 

for  E  <  Ej,,  and  where  p(E,)  is  the  level  density  at  energy  E,,.  In  the  constant  tem¬ 
perature  formula,  the  nuclear  temperature  is  related  to  the  level  density  parameter 
(a),  the  smooth  joining  energy  (E*),  and  the  pairing  correction  (A),  by  the  formula; 

E.  -  A 

T  =  - .  (1.4) 

^a(E.  -  A)  -  2 

In  the  lower-lying  energy  region  in  which  known  levels  exist,  Yamamuro  uses  these 
known  levels  in  the  statisticjd-model  calculations. 

In  the  statistical- model  code  GNASH,  which  was  used  in  Yamamuro's  work 
to  calculate  the  gamma-ray  spectrum,  transitions  of  the  El,  E2,  Ml,  and  M2  type 
were  permitted  for  the  gamma  rays.  The  Bnnk-.Yxel  strength  function  was  used  for 
the  El  and  Ml  strength.  The  Weisskopf  formulation  was  used  for  the  E2  and  M2. 
No  higher  multipolarities  were  considered. 

1.6  Review  of  Primary  Experimental  Work 

The  primary  literature  on  experimentally  determined  gamma- ray  energy  spec¬ 
tra  for  neutron  capture  in  can  be  attributed  to  the  work  of  Campion  et  al.  lOi  , 
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and  John  et  al.  111.  The  work  of  John  et  ai.  is  used  in  the  literature  as  a  benchmark 
for  the  comparison  of  ccilculated  spectra.  These  works  are  briefly  summarized  here 
in  order  to  allow  the  reader  to  gain  some  insight  into  the  methods  used  to  collect 
the  experimental  data. 

1.6.1  Campion  et  al. 

In  1959.  Campion  et  al.  lOi  examined  the  gamma  spectrum  resulting  from 
the  thermal  neutron  bombardment  of  uranium  metal,  depleted  to  0.034%  “^^U. 
Figure  1.1  llO  shows  the  smoothed  gamma  ray  spectrum  due  to  neutron  capture 
by  obtained  by  Campion.  The  gamma  rays  in  the  energy  region  between  0.14 
MeV  and  3.4  MeV  were  examined  with  a  crystal  spectrometer.  The  gamma  rays 
in  the  3.4  MeV  to  4.2  MeV  range  were  studied  using  a  pair  spectrometer  with  a 
resolution  of  about  1%.  A  survey  was  conducted  at  about  3%  resolution  for  the 
gamma-energy  range  from  3.0  MeV  to  7.7  MeV  also  with  the  pmr  spectrometer. 

The  absolute  intensities  for  gamma  rays  above  3.4  MeV  were  determined  by 
the  use  of  a  nickel  comparison  method  attributed  to  ICinsey  and  Bartholomew  19l. 
The  absolute  intensities  for  the  gamma  rays  in  the  3.4  MeV  and  below  region  were 
determined  by  direct  comparison  to  the  gamma  rays  of  known  intensity  emitted  from 
^^®Np  following  the  beta  decay,  which  are  in  equilibrium  with  the  (n,7)  reaction. 
Table  1.3  lists  the  energies  and  the  absolute  intensities  of  the  resolved  gamma- rays 
emitted  in  the  thermal  neutron  capture  in  ^U.  Note  that  the  intensity  is  reported 
on  a  number  per  single  capture  event  basis. 

In  Campion’s  work,  the  energy  radiated  per  capture  in  the  thermal  neutron 
capture  in  is  found  to  be  4.3  MeV.  No  gamma  rays  were  observed  above  4.062 
MeV.  In  this  work,  an  attempt  was  made  to  find  a  gamma-ray  in  coincidence  with 
the  4.062  MeV  gamma  ray.  Two  Nal(Tl)  scintillation  spectrometers  were  used  two 
observe  the  spectrum  from  the  bombardment  of  metal  with  thermal  neutrons. 
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Table  1.3:  Resolved  neutron  capture  gamma  radiation  from  neutron 
capture  in  from  Campion  et  al. 


Energy  (MeV) 

Intensity  (gamma  rays  per  capture) 

0.5401  =  0.004 

0.09 

0.5490  -  0.004 

0.09 

0.5772  *  0.004 

0.09 

0.5885  -  0.004 

0.05 

0.6093  -  0.004 

0.06 

0.6275  r:  0.004 

0.09 

0.6375  -  0.004 

0.03 

3.576  =  0.010 

0.02 

3.662  :r  0.018 

0.02 

4.062  r:  0.010 

0.07 

The  resulting  spectrum  indicated  that  the  4.062  MeV  gamma  ray  was  not  to  the 
ground  state. 


1.6.2  John,  et  al. 

In  1970,  J.  John,  V.J,  Orphan,  and  C.G.  Hoot  20]  reported  experimental 
results  on  the  epithermal- neutron  capture  gamma  rays  from  shield  materials,  par¬ 
ticularly  "^U.  They  measured  the  gamma  rays  from  resonant  capture  in  of 
neutrons  up  to  100  keV.  An  electron  linear  accelerator  was  used  to  produce  a  beam 
of  pulsed  neutrons  from  a  cylindrical  uranium  target.  The  6.0  inch  diameter  depleted 
uranium  target  wzu  located  at  the  end  on  a  16.0  meter  flight  tube.  The  neutrons 
incident  on  the  target  were  moderated  by  a  1.0  inch  thick  slab  of  polyethylene.  The 
gamma  rays  were  detected  using  a  Ge(LI)-NaI(Tl)  spectrometer.  This  detector  was 
operated  simultaneously  as  a  three-crystal  pair  spectrometer  and  an  anti-Compton 
spectrometer. 

Gamma-ray  pulse-height  spectra  were  generated  from  the  data  for  15  neutron 
energy  groups.  The  gamma-ray  intensities  were  determined  from  the  experimentally 
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Table  1.4;  Neutron  capture  gamma  radiation  from  neutron  capture 
in  from  John  et  ai. 


Energy  (MeV) 

Intensity  (gamma  rays  per  100  captures) 

1.00 

28.0 

1.25 

33.0 

1.50 

37.0 

1.75 

40.0 

2.00 

36.0 

2.25 

29.0 

2.50 

16.0 

2.75 

10.0 

3.00 

6.0 

3.25 

4.2 

3.50 

1.3 

3.75 

1.5 

4.00 

4.6 

4.25 

0.16 

4.50 

0.03 

4.75 

0.12 

collected  spectra.  The  intensities  were  reported  in  0.25  MeV  energy  bins.  No  exper¬ 
imental  data  is  given  for  below  1.0  MeV.  The  results  for  the  lowest  neutron  energy 
capture  (5.6  eV  <  E„  <  7.8  eV)  is  given  in  Table  1.4.  Note  that  the  intensity  is 
reported  on  a  number  per  one  hundred  capture  events  basis. 


CHAPTER  2 

ENERGY  LEVELS  OF  THE  NUCLEUS 

2.1  Nuclear  Levels 

For  a  potential  well  of  finite  depth,  the  levels  below  the  zero  energy  state  are 
referred  to  as  bound  states.  The  energy  levels  above  the  zero  energy  state  are  referred 
to  as  idrtuai  states.  In  the  single  particle  shell  model  of  the  nucleus,  protons  and 
neutrons  fill  the  bound  levels  of  the  nucleus  in  accordance  with  the  Pauh  exclusion 
principle,  which  states  that  two  protons  or  two  neutrons  in  the  same  nucleus  can 
not  have  the  same  set  of  quantum  numbers.  The  lower  energy  levels  are  fiUed  first. 
The  excited  states  of  the  nucleus  can  be  formed  by  the  outermost  nucleon  moving 
up  to  a  higher  energy  level.  Deexatation  of  a  bound  level  can  only  occur  by  gamma 
emission  or  by  the  competing  process  of  internal  conversion.  The  excitation  energy 
is  not  high  enough  in  these  bound  states  to  allow  for  particle  emission. 

In  the  lower-lying  levels,  the  energy  levels  are  spaced  far  enough  apart  to  be 
considered  discrete  levels,  and  their  width  is  small  enough  not  to  overlap  one  another. 
However,  in  the  high  energy  levels,  above  the  2  MeV  range,  the  single  particle  shell 
model  is  insufficient  in  detailing  the  energy  levels  and  widths.  .A.t  these  energies, 
the  levels  become  very  close  together,  and  the  widths  become  relatively  large.  The 
density  of  energy  levels  increa.ses  as  the  excitation  energy  is  increased.  Also,  at  these 
higher  energy  levels,  there  are  more  competing  modes  of  decay.  The  model  used  to 
describe  the  behavior  of  the  nucleus  at  these  higher  energies  is  based  on  statistical 
mechanics  and  thermodynamics.  ;21| 

Table  2.1  displays  some  discrete  energy  levels  and  the  spin-parity  assignments 
of  the  excited  nucleus,  as  shown  in  Booth  et  al.  and  the  Nuclear  Data 

Sheets. 


15 


6 


Table  2.1:  Discrete  energy  levels  in  up  to  0.950  MeV 


Level  Booth  et  ai. 


0.0  (2.5-) 
.043  (3.5-) 
.097  (4.5-f) 
.1331  (.5^) 
1428  (1.5-) 
.173  (3.5^) 
1915  (2.5-) 


.2:':  91  (4.5^) 
.3072  (5.5+) 
.683  (  5— ) 
.713  (1.5+) 
.7423  (2.5-) 
.7428  (1.5-) 
.750  (.5-) 
.7986  (3.5-) 
.811  (.5-) 
.8191  (.5-) 
.823  (1.5-) 
.860  (2.5-) 
.943  (4.5-) 


Nuclear  Data  Sheets 


0.0  (2.5-) 
.1338  (.5,1.5-) 
.1456  (.5.1.5—) 
.1839  (2.5-) 
.6877  (.5,1.5-) 
.7157  (.5,1.5—) 
".'00  (.5.1.5—) 
.7347  — 
.7392  (.5,1.5-) 
.7459  (.5,1.5-) 
.7567  (2.5-) 
~.800  — 
.8150  (.5,1.5-) 
.8236  (.5.1.5-) 
.8530  (.5,1.5—) 
.8882  (2.5-) 
.8933  (2.5-) 
9325  (.5,1.5-) 


The  usuai  convention  employed  in  dealing  with  level  densities  is  to  discuss  the 
density  of  levels  instead  of  the  density  of  states.  The  level  density  is  the  number 
of  levels  per  MeV  of  energy.  In  the  lower  lying  energy  region,  the  level  density  is 
relatively  small,  cind  in  the  higher  energy  region,  the  level  density  is  relatively  large. 
The  energy  level  densities  up  to  a  few  MeV  can  be  obtained  experimentally  from  the 
determination  of  the  population  of  levels  in  nuclear  reactions  or  radioactive  decay. 
Slow-neutron  resonances  can  reveal  information  concerning  the  level  densities  in  the 
6-8  MeV  range.  Higher  levels  require  the  use  of  statistical  mechanics  to  determine 
the  level  densities. 

There  is  a  degeneracy  of  2J  -i-  1  states  in  a  level  with  spin  J.  The  physical 
observable  is  the  energy  level  spacing  (D),  which  is  the  reciprocal  of  the  energy  level 
density  (p).  The  evaporation  theory  of  nuclear  decay  predicts  that  the  level  density 
as  a  function  of  energy  behaves  as: 

P(E)  =  ^  =  Po  *  exp  ;aE‘2j  ,  (2.1) 

where: 

•  p{E)  =  level  density  at  energy  E 

•  D  =  spacing  between  levels  (MeV) 

•  Pg  =  I  to  0.01  (MeV)~^  for  A  =  20  to  200 

•  a  =  1  to  7  (MeV)“:  for  .4  =  20  to  200 

The  constant  (a)  is  called  the  level  density  parameter  and  is  proportional  to 
the  atomic  weight.  The  value  of  p^  is  determined  by  the  level  density  at  a  low  (about 
1  Mev)  energy  level  and  the  above  expression  is  adjusted  for  the  higher  energy  levels 
by  the  level  density  parameter.  The  nucleus  can  also  be  characterized  221  by  its 
nuclear  temperature  (t),  which  is  related  to  the  excitation  energy  of  the  nucleus  by 
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the  relation; 


E  =  at^ 


(2.2) 


2.2  Bethe’s  Free  Gas  Model 

In  the  Bethe  free  gas  model  ;23i,  the  nucleus  is  considered  to  be  a  Fermi  gas  of 
free  neutrons  and  protons.  These  particles  are  confined  to  the  volume  of  the  nucleus. 
In  this  model,  the  level  spacing  (inverse  of  the  level  density)  is  found  to  depend  on 
the  mass  number  (A)  and  the  excitation  energy  (E). 

The  spacing  between  two  levels  of  a  given  spin  I,  in  the  Bethe  free  gas  model, 
is  given  by  the  relationship: 
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D  «  121_2  ’  C 


21  -  1 


(2.3) 


where; 


•  c  =G(1  -(’rVl2)(^Co)-^  *  •••). 

•  Co  =  (N/C)1. 

•  5“^  =  average  excitation  energy  of  the  individual  particles. 

•  N  is  the  total  number  of  particles. 

.  c  =  (|i)(^)i 

•  X  =  (AE/2.20)^  for  exatation  energy  E  in  MeV, 


•  M  =  particle  meiss. 

•  R  =  nuclear  radius. 

Since  the  level  density  is  the  reciprocal  of  the  level  spacing  we  can  state  that  the 
level  density  in  the  Bethe  formulation  is  simply: 

5 


>  D  432  Vlog2 


(2/  - 


(2.4) 
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Table  2.2:  Spacing  of  nuclear  energy  levels  of  zero  angular  momen¬ 
tum  in  eV 


E*A  (MeV) 

spacing  (MeV) 

100 

1x10' 

200 

2.4x10® 

400 

1.9x10® 

600 

2.1x10'' 

800 

2800 

1000 

450 

1200 

85 

1500 

8.5 

1800 

1.0 

The  spacing  D^,  is  the  spacing  between  levels  of  zero  angular  momentum  and  is 
simply  the  value  of  equation  2.3  times  (21—1).  Table  2.2  displays  some  representative 
values  of  the  level  spacing  for  given  values  of  the  product  of  the  excitation  energy  (E) 
and  the  mass  number  (A).  As  seen  in  the  table,  the  level  spacing  between  successive 
levels  decreases  rapidly  with  excitation  energy  for  a  given  mass  number.  If  one 
compares  different  mass  number  elements,  then  for  a  given  excitation  energy,  the 
level  spacing  decreases  rapidly  with  increasing  mass  number. 

For  many  years  experimental  observations  of  nuclear  level  densities  have  been 
interpreted  in  the  framework  of  the  Fermi  gas  model,  and  most  statistical  model 
calculations  have  employed  level  density  formulas  of  the  type  introduced  by  Bethe, 
where  the  level  density  exponentially  increases  with  both  the  excitation  energy  and 
the  mass  number.  This  simple  model  does  not  take  into  account  the  odd-even 
and  shell  effects  so  modifications  to  the  original  level  density  formula  have  been 
proposed.  More  modem  level  density  formulations  are  now  discussed  and  their 
details  examined  in  the  context  of  the  present  gamma  cascade  model. 
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2.3  Gilbert  and  Cameron  Formulation 


In  1965,  Gilbert  and  Cameron  24'  developed  a  composite  nuclear-level  density 
formula  which  contained  shell  corrections.  In  this  formulation,  a  constant  nuclear 
temperature  model  is  used  for  the  lower  lying  nuclear  level  densities.  The  regular 
Fermi  gas  formula  is  used  at  the  higher  excitation  energies.  The  two  methods  are 
joined  by  fitting  the  curves  tangentially. 

The  level-density  relation  adopted  for  the  higher  energy  levels  is  essentially 
similar  to  Bethe’s  in  that  there  exists  an  dependence  on  the  level  density. 


This  level  density  formula  is  given  by  the  relationship: 


p(U,J)  = 


_  v^exp(2vTU)  (2J  -r  l)exp(-(J  --|)^/2g-^) 


12  aU’ 


(MeV-^),  (2.5) 


where: 


/j(U,J)  =  level  density  (in  MeV  '). 


•  J  =  angular  momentum. 


U  =  energy  above  fully  degenerate  state  (which  is  at  excitation  energy  Uo). 


•  a  —  spin-dependence  parameter. 


Note  that  this  representation  accounts  for  the  spin  J.  The  observable  level 


density  is: 


.(U)  =  '  (M=V-'). 

/  a«  U  <  V  ^TTcr 


This  observable  level  density  is  not  the  same  tis  the  total  level  density  which  also 
includes  those  levels  degenerate  in  the  magnetic  quantum  number,  so  the  total  level 
density  W(U)  is  a  factor  of  (2J-t-l)  times  larger  than  the  observable  level  density 
p(U).  This  total  level  density  is  given  by: 

w{u)  =  y;(2j  -  1WU.J)  =  (2.7) 
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In  order  to  use  this  expression,  one  must  know  the  pairing  corrected  energy 
(U),  the  spin  dependence  parameter  (o-),  and  the  level  density  parameter  (a).  Note 
that  in  the  Gilbert  and  Cameron  formulation,  the  energy  U  is  defined  by  the  total 
excitation  energy  E,  corrected  for  pairing  effects  by  subtracting  P(N)  and  P(Z), 
which  are  the  pairing  energies  for  neutrons  and  protons  respectively. 

U  =  E  -  P(N)  -  P(Z).  (2.8) 

The  spin  dependence  parameter  is  given  by  the  relation; 

=  0.0888 -a -t -A^  (2.9) 

where  t  =  ^/u/a  and  is  the  thermodynamic  temperature  of  the  nucleus.  Therefore, 
a  varies  as  ,  which  is  a  slow  change  with  an  increase  of  excitation  energy. 

The  only  free  parameter  left  in  the  formula  is  the  level  density  parameter.  The 
value  used  for  the  parameter  for  undeformed  nuclei  is  given  by: 

4  =  0.00917 -S  -  0.142,  (2.10) 

A 

and  for  deformed  nuclei  by; 

4  =  0.00917-5  -  0.120,  (2.11) 

A 

where  S  =  S(N)  —  S(Z),  the  total  shell  correction  (in  MeV),  where  S(N)  represents 
the  shell  correction  for  neutrons  and  S{Z)  represents  the  shell  corrections  for  protons. 
As  can  be  seen  by  the  given  relationship,  the  ratio  of  the  level  density  parameter  to 
the  atomic  weight  (a/A)  is  not  a  constant  as  in  the  earlier  models  such  as  the  Bethe 
model,  which  did  not  account  for  the  shell  effects. 

The  lower  levels  of  the  nuclei  are  generally  known.  In  Gilbert  and  Cameron’s 
work,  they  fitted  the  lower  lying  levels  to: 


N(E)  =  expi(E  -  E„)/T;. 


(2.12) 
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where  Eo  and  T  are  parameters  determined  by  the  fitting  of  the  upper  and  lower 
level  density  formulas  at  the  point  of  tangency. 

The  density  of  the  lower  levels  may  be  represented  by: 

N(E) 

(2.13) 

and  the  density  of  the  higher  levels  is  given  by  Equation  2.6.  The  two  curves  must 
be  fitted  tangenticilly.  The  conditions  at  the  point  at  which  the  two  densities  meet 
are 


*  P lower  P upper- 

•  T  =  nuclear  temperature. 

This  two  level  formulatioa  is  termed  a  composite  level  density  formula.  There 
are  some  problems  with  it,  however.  Some  nuclei  will  not  have  a  good  fit  between 
the  upper  and  lower  densities.  This  may  be  a  result  of  the  absence  of  information 
about  the  low-lying  level  structure.  The  nuclei  with  A  <  40  and  the  nuclei  with  X 
or  Z  within  2  units  of  a  magic  number  tend  to  give  this  composite  formula  problems. 
Nuclei  near  the  magic  numbers  have  a  much  smaller  number  of  lower  levels  up  to 
a  given  energy  than  do  those  nuclei  far  awqy  from  the  magic  numbers.  Therefore 
fitting  the  upper  and  lower  level  densities  at  the  tangency  point  for  the  magic  nuclei 
is  not  as  simple  as  for  the  others.  Gilbert  and  Cameron  leave  the  reader  with  the 
warning  that  the  error  in  using  this  composite  level  formula  is  much  greater  near 
the  magic  number  nuclei.  24( 

2.4  The  Back-Shifted  Fermi-Gas  Model  of  Dilg  et  al. 

Dilg  et  aJ.  ,251  proposed  a  parameterization  of  the  back-shifted  Fermi  gas 
model  proposed  by  Vonach  and  Graves  261  in  the  whole  meiss  range.  In  this  ap¬ 
proach,  the  Fermi  gas  model  is  adopted  with  adjustable  parameters  representing 
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the  Active  ground  state  position  (A)  and  the  level  density  parameter  (a).  These 
parameters  are  adjusted  by  the  fitting  of  low  energy  and  resonance  data.  This  pa¬ 
rameterization  provides  a  reasonable  description  of  level  densities  up  to  10-20  MeV. 

In  the  Dilg  approach  the  following  forms  of  the  spin-dependent  and  total  level 
densities  are  adopted  from  Lang  and  Le  Couteur.  27' 

a-  II  -  1  (2J-l)exp:2(a(U  -A))-^  -  J(J  -  I)/2^^ 

’  24v^  (U  -  A  -  tj'i 

and 

1  1  expi2(a(U  —  A))"2 

12v^o-a^  (U-A-t)i  ’ 
where  the  thermodynamic  temperature  is  defined  by: 


(2.14) 

(2.15) 


U  -  A  =  at^  -  t, 

where  t  is  the  nuclear  temperature  and  a  and  A  are  the  fitted  parameters.  The  Active 
ground  state  position  is  shifted  by  the  amount  A  because  two-body  interactions 
between  the  nucleons  have  been  ignored. 

These  level  density  expressions  are  derived  from  the  total  density  of  states 
assuming  random  couphng  of  angular  momenta.  In  this  treatment,  an  additional 
parameter  is  introduced.  This  spin  cut-off  parameter  is  used  to  describe  the  spin 
distribution.  At  high  excitation,  the  value  of  the  spin  cut-off  parameter  is  expected 
to  approach: 

<^«^~  00150A-n.  (2.16) 


The  back-shifted  Fermi  expression.  Equation  2.14  can  be  represented  by  a 
simpler  expression  for  28  i 

exp(-J(J  l)/2<r|)exp(2a’j(U  -  A)‘2) 


p(U,J)  =  G(2J  -  1)- 


(U  -  A)^ 


860 


(2.17) 


where  the  value  of  G  =  0.00010091,  a  =  9.20212,  A  =  -  0.58  and  a  =  3.7283. 
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2.5  Comparison  of  the  Level  Density  Models 

The  Bethe  formulation  has  an  exponential  incretise  in  the  level  density  with 
both  the  excitation  energy  and  the  mtiss  number.  This,  however,  does  not  represent 
the  low  energy  levels  well.  The  constant  temperature  model  predicts  the  low  energy 
levels  well,  but  does  not  represent  the  high  energy  levels  well.  The  Dilg  back-shifted 
model  takes  into  account  the  odd-even  and  shell  effects,  resulting  in  a  more  realistic 
prediction  of  the  level  densities  in  the  low  and  high  energy  regions.  The  model  used 
in  this  work  is  a  modification  of  the  Dilg  back  -shifted  model,  and  agrees  well  with 
the  other  models. 

Figure  2.1  through  Figure  2.3  display  level  density  versus  energy  for  the  present 
work,  the  Dilg  back-shifted  model,  and  the  Fermi  gas  model.  The  level  density  is 
shown  in  energy  levels  per  MeV.  Figure  2.1  displays  the  level  density  for  J  =  l/2. 
Figure  2.2  displays  the  level  density  for  J=3/2.  Figure  2.3  displays  the  level  density 
for  J=5/2.  Note  that  the  level  density  scale  is  logarithmic. 
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level  OENSin  vs.  ENE-=Gt,  J-1/2 
CCr'.PRRISON  QF  OIFFE.RENT  rORMULRT'CNS 


Figure  2.1:  Level  density  versus  energy  for  J=l/2 
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level  density  vs.  energy,  J-3/2 


Figure  2.2:  Level  density  versus  energy  for  J=3/2 
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'.EVEL  OENSIT'f  VS.  EME='?'.  J-5/2 


Figure  2.3:  Level  density  versus  energy  for  J=5/2 


CHAPTER  3 

GAMMA  TRANSITIONS 


A  nucleus  can  be  excited  in  the  neutron  capture  reaction.  For  example,  in  the 
^^U(n.7)^^®U  reaction,  and  a  neutron  will  form  an  excited  nucleus. 

238^-  ,  _239 

This  excited  nucleus  can  then  deexcite  through  gamma  decay,  or  by  the  competing 
process  of  internal  conversion.  Gamma  decay  is  the  process  by  which  a  nucleus 
gives  up  its  excess  energy  by  the  release  of  electromagnetic  radiation.  The  internal 
conversion  process  is  the  result  of  the  interaction  of  the  nucleus  and  orbital  electrons, 
leading  to  the  emission  of  an  electron  with  kinetic  energy  equal  to  the  excitation 
energy  of  the  nucleus  minus  the  binding  energy  of  the  electron  in  the  atom. 

Ee-  =  Eg*  —  Efeitui,  (3.1) 

where  E*-  is  the  electron’s  kinetic  energy  upon  emission  from  the  atom,  Eg*  is  the 
excitation  energy  of  the  nucleus,  and  E{„n<j  is  the  binding  energy  of  the  emitted 
electron  in  the  atom. 

3.1  Energetics  Of  Gamma  Decay 

By  using  conservation  of  momentum  and  conservation  of  energy,  one  may 
compare  the  relative  amounts  of  energy  shared  among  the  products  of  a  decaying 
nucleus.  First,  assume  that  the  rest  mass  of  the  excited  nucleus  is  m^  and  the  rest 
mass  of  the  final  state  is  m,,.  Then  by  conservation  of  energy: 

rajc"  =  m^c^  -  E,  -  T/,  (3.2) 

and  by  conservation  of  momentum: 

p,  -  pf  =  0.  (3.3) 
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where: 

•  =  Gamma-ray  energy. 

•  T I  =  Energy  of  the  recoil  nucleus. 

•  =  Momentum  of  the  gamma  ray. 

•  P/  =  Momentum  of  the  recoil  nucleus. 


and  since 


T  f  —  -m„v^ 


*> 

P/ 

2m„’ 


P/  =  P->, 


(3.4) 


(3.5) 


the  final  energy  of  the  recoil  nucleus  may  be  expressed  in  terms  of  the  gamma-ray 
momentum, 


T,  = 


2m„ 


(3.6) 


Here,  the  relationship  between  momentum  and  energy  may  be  used  to  put  this 
relationship  in  terms  of  the  gamma- ray  energy. 


T/ 


2m<,c^ 


(3.7) 


Using  realistic  values  of  gamma- ray  energy  (E  =  2  MeV)  and  the  atomic  mass 
number  (A  =  239)  for  the  calculated  value  of  the  nucleus  recoil  energy  is 

approximately  9.0  eV.  This  is  only  0.0005%  of  the  value  of  the  gamma-ray  energy, 
and  may  be  considered  negligible.  Therefore,  the  gamma-ray  may  be  considered  to 
carry  off  all  of  the  decay  energy  of  the  transition.  ,29] 
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3.2  Level  Width 


The  level  width  is  related  to  the  mean  lifetime  of  the  of  the  level.  Planck’s 
constant  divided  by  the  hfetime  of  the  level  is  defined  as  the  width  of  the  level. 


r 


(3.8) 


The  uncerlairty  principle  requires  the  energy  level  to  have  some  finite  width.  This 
is  a  result  of  not  being  able  to  precisely  know  the  energy  and  time  at  the  same 
moment.  The  level  width  (r)  is  equal  to: 


r 


0.66  X  10~^^ 

mean  lifetime  of  width 


(eV). 


(3.9) 


In  nuclei  which  have  competing  modes  of  decay,  the  width  of  the  level  is  simply  the 
sum  of  the  different  level  widths  for  each  mode  of  decay. 


3.3  Gamma  Transitions 

When  a  gamma-ray  is  emitted  from  a  nucleus,  it  carries  with  it  some  amount 
of  angular  momentum.  The  value  of  this  angular  momentum  is  quantized,  and  is  an 
integral  multiple  of  Planck’s  constant  (L^,  where  L  is  an  integer  greater  than  or  equrd 
to  1).  For  each  ^'alue  of  L,  there  exists  two  different  characteristic  gamma  radiations. 
There  is  an  electric  multipole  and  a  magnetic  multipole,  each  with  different  parities. 

When  the  charge  in  the  nucleus  oscillates,  a  characteristic  radiation  is  emitted. 
The  oscillation  of  an  electric  dipole  (two  charges  of  equal  magnitude  separated  by 
a  distance)  results  in  the  emission  of  a  specific  pattern  of  radiation,  given  the  title 
of  electric  dipole  radiation.  A  circulating  charge  (current  loop)  will  emit  another 
characteristic  pattern,  given  the  title  of  magnetic  dipole  radiation. 

One  major  difference  between  the  radiation  emitted  in  the  classical  sense  and 
the  radiation  emitted  in  nuclear  decay  is  the  parity  operation  in  going  from  one  nu¬ 
clear  level  to  another.  In  gamma  decay,  each  multipole  moment  of  order  L  produces 
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radiation,  taking  with  it  some  angular  momentum  (I).  From  the  law  of  conservation 
of  angular  momentum,  it  is  required  that: 

L  =  (3.10) 

where  the  subscripts  indicate  the  initial  and  final  states  of  the  nucleus.  The  parity 
change  is  related  to  the  angular  momentum  carried  off.  as  shown  by  the  relationships; 

•  For  electric  multipole  radiation,  parity  change  All  =  (  —  1)'^. 
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For  the  first  four  electric  and  magnetic  multipoles,  the  transition  probabilities 
can  be  shown  to  be: 

A(E1)  =  1.0x10^\v'e^ 

A(E2)  =  7.3x10^V3E= 

A(E3)  =  34A2E" 

A(E4)  =  1.1  X  10~^A.'3E®  (3.16) 

A(M1)  =  5.6  X  lO'^E"* 

A(M2)  =  3.5x10^4'3E“ 

A  (M3)  =  16A3E^ 

A(M4)  =  4.5  X  10-®A-E® 

As  the  multipole  order  increases,  the  decay  constant  (A)  decreases,  and  this  can 
be  interpreted  as  an  increase  in  the  half-life,  or  a  smaller  probability  of  the  nuclear 
transition  decaying  by  that  mode.  The  probability  for  an  electric  multipole  radiation 
can  be  shown  to  be  greater  than  for  a  magnetic  multipole  of  the  same  order  by  simply 
dividing  the  two  disintegration  constants.  X  graph  for  the  half-hfe  versus  gamma- 
ray  energy  is  shown  for  X  =  55.  This  graph  summarizes  the  relative  likelihood  of 
electromagnetic  radiation  being  emitted  in  a  particular  multipole  transition.  For 
a  given  order  multipole,  it  is  seen  from  the  graph  that  the  relative  probability  of 
having  an  electric  transition  is  an  order  of  magnitude  larger  (the  magnitude  of  the 
half-life  is  an  order  of  magnitude  smaller)  than  the  probability  of  having  a  magnetic 
transition. 

Note  that  this  graph  is  a  log-log  plot.  29] 

3.4  Selection  Rules 

Classically,  the  electromagnetic  field  produced  by  an  oscillating  charge  (or 
current)  transmits  energy  and  angular  momentum.  .4  multipole  of  order  L  transfers 
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an  angular  momentum  of  Lft  per  photon.  Consider  a  nucleus  in  some  initial  state 
(li,  and  TT,)  which  decays  to  a  final  stale  (1^.  and  tt j).  Conservation  of  total  angular 
momentum  states  that  the  initial  angular  momentum  of  the  nucleus  must  be  equal 
to  the  final  cingular  momentum  of  the  nucleus  plus  that  angular  momentum  carried 
off  by  the  gamma  ray. 

I.  =  L-I^.  (3.17) 

The  possible  values  of  L  are  restricted. 

-I/!<L<I.^I/.  (3.18) 

The  dominant  transition  will  have  L  =}  I,-  —  1/  i  because  of  the  strong  L  dependence 
of  the  transition  rate. 

For  example,  if  T  =  |  and  ly  =  ',  then  L  can  be  1,  2,  3,  or  4  as  determined 
by: 


3 

2 


5 

2 


5 

2' 


(3.19) 


The  type  of  dominant  radiation  emitted  (electric  or  magnetic)  depends  on  the 
value  of  L  (odd  or  even)  and  the  relative  parity  of  the  initial  and  final  states  of  the 
nucleus. 

For  electnc  transitions; 


AH  =  NO  ^  radiation  will  have  even  L. 
AH  =  YES  =>  radiation  will  have  odd  L. 


For  magnetic  transitions; 

AH  =  NO  =>  radiation  wiD  have  odd  L. 

An  =  \’ES  =>  radiation  will  have  even  L. 

That  is,  if  T;  =  Ty,  (An  =  NO)  for  the  example  given  above  (I,  =  and 
ly  =  ^)  ,  then  the  radiation  will  be  Mi,  E2,  M3,  or  E4. 
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Table  3.1:  Classification  of  gamma  rays 


Classification 

Symbol 

L 

Parity  Change  in  Nucleus 

electric  dipole 

El 

1 

yes 

magnetic  dipole 

Ml 

no 

electric  quadrapole 

E2 

2 

no 

magnetic  quadrapole 

M2 

2 

yes 

electric  octupole 

E3 

3 

yes 

magnetic  octupole 

M3 

3 

no 

electric  2L-pole 

EL 

L 

no  for  L  even,  yes  for  L  odd 

magnetic  2L-pole 

ML 

L 

yes  for  L  even,  no  for  L  odd 

In  the  case  =  0,  the  selection  rules  predict  L  =  0  which  is  not  permitted 

in  radiative  transitions.  Internal  conversion  is  allowed,  however. 


3.5  Classification  of  Gamma-Rays 

Table  3.1  summarizes  the  classification  of  gamma  transitions  301.  The  classi¬ 
fication  is  based  on  the  angular  momentum  (L)  and  the  type  of  emission,  that  is. 
electric  or  magnetic. 

The  symbol  E  is  for  the  electric  transitions  and  the  symbol  M  is  for  the  mag¬ 
netic  transitions.  Note,  that  for  a  given  transition,  the  classification  of  the  transition 
(multipolarity  order  and  type)  is  dependent  on  the  parity  change  between  the  initial 
and  final  levels,  and  the  angular  momentum  difference  between  the  initial  and  final 
levels. 


3.6  Internal  Conversion  and  Internal  Conversion  Coefficients 

In  a  nuclear  transition,  the  energy  released  can  be  emitted  as  a  photon  with 
energy  E  =  hi/  or  it  can  be  transferred  to  an  orbital  electron  which  is  ejected  with 
an  energy  equal  to  hi/  -  E(,j^(j),  where  E^ridH)  is  equal  to  the  binding  energy  of  the 
electron  in  the  i**  shell.  This  process  of  ejecting  an  electron  in  order  to  emit  the 
transition  energy  is  called  interned  conversion.  For  all  electromagnetic  transitions. 
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except  for  the  L=0  to  L=0  transition,  gamma  emission  and  internal  conversion 
compete.  When  considering  gamma  decay  and  internal  conversion  as  the  only  decay 
processes,  the  decay  constant  may  be  written: 


(3.20) 


where  is  the  probability  for  emission  of  a  gamma  ray  per  second,  and  A^-  is  the 
probabihty  for  emission  of  a  conversion  electron  per  second. 

The  internal  conversion  coefficient  (a)  is  defined  as  the  ratio  of  the  probabihty 
of  electron  emission  to  the  probabihty  of  gamma  emission. 


a 


(3.21) 


This  ratio  is  a  measure  of  the  total  number  of  conversion  electrons  emitted  over  a 
time  penod  divided  by  the  number  of  gamma  rays  emitted  in  the  same  transition 
over  the  same  period  of  time.  The  internal  conversion  coefficient  is  dependent  on 
the  shell  from  which  the  electron  is  emitted.  The  total  probabihty  for  internal 
conversion,  (A,-),  can  be  written  ais  the  sum  of  the  partial  probabihties  for  ejection 
in  the  K.  L,  M....  shells. 


(3.22) 


The  partiai  coefficients  are  computed  for  each  sheU  and  are  summed  to  obtain  the 
toted  conversion  coefficient. 


Q  =  Ofc  ~  ai  ~  am  —  . . .  .  (3.23) 

Consideration  of  only  the  K-shell,  L-sheil,  and  M-shell  is  suffiaent,  since  the  partial 
coefficients  decrease  rapidly  with  the  increase  in  the  prinapal  quantum  number. 

In  addition  to  depending  on  the  shell  from  which  the  electron  is  emitted, 
conversion  coeffiaents  depend  on  the  atomic  number  (Z)  of  the  nucleus  and  on  the 
transition  energy  and  multipolarity  of  the  radiation  which  the  conversion  competes 
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wit'  The  conversion  coefficients  are  presented  in  tables  as  functions  of  atomic 
number,  electron  shell,  transition  energy,  and  multipolarity.  31 

3.7  Transition  Probabilities 

The  transition  probability  for  a  specific  type  of  transition  (e)  may  be  deter¬ 
mined  by  dividing  the  transition  rate  by  the  sum  of  all  the  transition  rates. 

probabihty  of  transitione  =  (3-24) 

This  is  simply  a  normalization  procedure  to  result  in  the  fraction  of  a  specific  type 
of  transition  occurring  between  two  states.  The  theoretical  estimates  for  the  transi¬ 
tion  rates  (disintegration  constants)  are  only  approximate.  Therefore  the  transition 
probabilities  are  also  only  approximate. 

As  previously  stated,  the  probability  of  a  nucleus  emitting  a  gamma- ray  de¬ 
creases  with  increasing  value  of  L.  This  means  that  2L-pole  radiation  is  more  prob¬ 
able  than  2L-f- 1-pole  radiation.  Also  the  theory  predicts  that  the  intensity  of  the 
electric  2L-pole  radiation  is  more  intense  than  the  magnetic  2L-pole  radiation.  How¬ 
ever,  this  is  not  necessarily  the  case  in  actual  decay  histories. 

3.8  Yrast  levels 

The  lowest  energy  level  for  each  spin  is  called  an  yrast  level.  The  term  was 
suggested  by  J.R.  Grover  and  is  a  Swedish  word  meaning  “dizziest”  22!.  Most  of 
the  energy  in  such  a  state  is  in  the  form  of  rotational  energy  and.  therefore,  the 
intrinsic  shape  of  the  nucleus  may  be  close  to  the  equilibrium  value  in  the  ground 
sta.e  region.  Because  of  this,  the  various  yrast  levels  of  a  nucleus  share  the  same 
intrinsic  state  cind  differ  only  in  their  rotational  motion. 

Because  of  the  special  position  of  the  yrast  levels,  it  is  common  for  a  highly 
excited  nucleus  to  decay  by  particle  emission  down  to  the  yrast  level.  From  the  yrast 
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level,  the  nucleus  will  decay  to  a  lower  state  by  gamma  emission  and  then  continue 
the  decay  process  by  gamma  decay,  through  lower  levels,  down  to  the  ground  state. 
•32! 

For  an  yrast  level  at  excitation  energy  E„,  the  maximum  kinetic  energy  of  an 
emitted  particle  is 

Epaw  =  Ec  -  -  Ej,,  (3.25) 

where 

•  Epart  =  the  maximum  kinetic  energy  of  the  particle  being  emitted. 

•  Ec  =  the  excitation  energy  of  the  compound  nucleus. 

•  Ej^ind  =  the  binding  energy  of  the  particle  being  emitted. 

•  Ej,  =  the  energy  of  the  yrast  level  under  consideration. 

It  is  possible  for  the  value  of  Ec  —  ^bind  ~  Ey  to  be  less  than  zero  while  Ec  —  E(,j,^ 
is  greater  than  zero.  In  this  case,  the  particle,  assuming  it  is  spinless  and  emitted 
with  zero  angular  momentum,  will  not  be  emitted  from  the  compound  nucleus  even 
though  the  excitation  energy  of  the  nucleus  is  greater  than  the  binding  energy  of  the 
particle  to  be  emitted.  The  quantity  Ey  in  effect  raises  the  level  at  which  gamma 
emission  will  compete  effectively  with  particle  emission.  The  increEise  is  above  the 
particle  emission  threshold  by  severed  million  electron  volts.  [32j 

3.9  Hindrance  Factors 

The  single- particle  model  assumes  that  the  nucleus  consists  of  an  inert  core 
made  up  of  completely  filled  shells.  In  this  model,  the  outer  shell  is  not  completely 
filled.  The  single- particle  theory  considers  the  inner  shells  to  move  in  independent 
orbits,  and  that  the  outer  shell  determines  all  of  the  low  energy  characteristics  of 
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the  nucleus.  The  interaction  between  a  given  nucleon  and  all  others  is  considered 
to  be  a  spherically  synunetric  potential. 

The  hindrance  factor  is  defined  as  the  ratio  of  the  experimentally  determined 
mean-life  (r„)  for  gamma  decay  to  the  calculated  mean-life  ("cic)  for  gamma  decay 
using  the  single- particle  model. 

Hindrance  Factor  =  — (3.26) 

^calc 

Since  the  single-particle  model  does  not  take  into  account  all  of  the  collective 
interactions  of  the  nucleons  in  the  nucleus,  the  experimentally  determined  lifetime 
of  the  nucleus  does  not  agree  with  the  computed  lifetime.  This  hindrance  factor 
adjusts  the  computed  lifetime  to  agree  with  the  experimentally  determined  lifetime. 
The  El  transition,  expected  to  be  the  most  probable  transition,  is  actually  very 
hindered.  Figure  3.2  (331  presents  the  hindrance  factors  for  the  El  transition  in 
heavy  nuclei.  In  heavy  nuclei,  the  El  transition  is  slowed  down  considerably  more 
than  in  the  light  nuclei.  In  the  larger  nuclei,  states  of  opposite  parity  differing  by 
one  unit  of  ang\ilar  momentum  are  more  hkely  to  be  verj'  different  in  structure. 
The  shell  model  does  not  provide  the  opportunity  for  the  occurrence  of  such  states 
differing  in  the  level  of  only  one  particle. 

In  light  nuclei,  the  Ml  transition  rate  agrees  quite  well  with  the  calculated 
single  particle  model  expectations.  In  the  heavy  nuclei,  the  Ml  trcinsition  is  fairly 
consistent  with  the  single  particle  model  calculated  values,  however  in  the  distorted 
nuclei,  the  Ml  transition  is  slightly  hindered.  In  contrast  with  other  multipoles, 
the  E2  transition  is  usually  slightly  enhanced  over  the  single  particle  values  in  both 
the  light  and  heavy  nuclei.  Figure  3.3  [33]  presents  the  hindrance  factors  for  the 
E2  transition  in  heavy  nuclei.  This  is  a  result  of  the  enhancement  of  the  collective 
quadrapole  moments  over  the  single  particle  quadrapole  moments.  .A.s  a  result  of  the 
enhancement  of  the  E2  rate,  it  is  expected  that  the  Ml  and  the  E2  transition  will 
favorably  compete  with  the  El  transition,  determined  by  the  single  particle  model 
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to  be  dominant.  '34l 

3.10  Angular  Distribution  of  Gamma  Rays 

In  a  gamma  cascade,  there  exists  a  correlation  between  the  direction  of  the 
first  gamma  emitted  and  the  direction  of  the  second  gamma  emitted.  The  angular 
distribution  for  the  second  emitted  gamma  relative  to  the  first  emitted  gamma  is  the 
same  for  electric  and  magnetic  radiation  of  the  same  multipole  order  (1)  and  same 
quantum  number  (m),  which  is  associated  with  the  z-component  of  the  angular 
momentum.  Therefore,  the  multipole  order,  but  not  the  parity  of  the  radiation, 
may  be  determined. 

The  angular-correlation  function  describes  the  distribution  of  angles  between 
successive  gammas  emitted  in  a  two  step  cascade.  The  angular-correlation  function 
W(^)  for  the  angle  9  between  two  consecutive  gamma  rays  in  the  cascade 

Jl  — +  J  — ♦  J2) 

can  be  shown  to  be: 

i=L 

W{e)dQ  =  Y,  ‘^2iP2i{cos9)dQ,  (3.27) 

i=0 

where  A2i  are  coefficients  which  depend  on  the  multipole  orders  of  the  two  gamma 
rays  and  is  given  by; 

Aj.  =  F2.(Ai, Ai,  Jl,  J)F2.(A2,  Az,  Jj,  J),  (3.28) 

where  F2i  are  certain  combinations  of  Clebsch- Gordon  coefficients  and  Racah  coef¬ 
ficients,  and  P2i(cos^)  are  the  even  Legendre  polynomials. 

One  may  analyze  the  angular  distribution  of  the  second  gamma  ray  in  a  two 
step  cascade  by  determining  the  distribution  of  the  energy  flow  from  a  radiating 
system.  The  Poynting  vector  S,  which  is  defined  by  the  relationship: 

S  =  —  E  X  H. 

47r 


(3.29) 
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describes  the  energy  flow  in  an  electromagnetic  wave.  If  a  large  sphere  is  drawn 
around  the  radiating  system,  then  at  any  point  on  the  sphv:r>>  S  is  the  energy 
escaping  per  square  centimeter  per  second. 


S:^— (3.30) 

47r  471 

The  energy  emitted  into  the  sohd  angle  dfi  around  the  direction  {9.  0)  is  the  product 
of  S  and  the  area  R^dfl  subtended  by  the  sohd  angle  at  the  sphere  surface. 

The  electric  multipole  field  of  order  l,m  at  point  r  =  (r.  9,  0)  is  given  by: 

Eg  = -curh-^^^-^X,,^(^,(p)i,  (3.31) 

K,  KT 


and 


He 


ur(r) 

KT 


(3.32) 


where  the  outgoing  wave  is  designated  by  u*.  The  magnetic  multipole  field  of  order 


I,m  is  given  by; 

E^  = 

KT 

Hm  =  -i.curlf^X,^(^,<p)| 

K  KT 

where  X(  „,(^.  d>)  are  the  vector  spherical  harmonics  defined  by: 


(«,</>) 


-zr  X  VYi_r„{9,0) 

- 1) 


(3.33) 

(3.34) 


(3.35) 


and  Y/m  are  the  scalar  spherical  harmonics.  The  angular  distribution  function 
Zi„i(9.<p)  is  defined  by: 


or 


ZU9,0}=X,,^(9.0)^X,,U9,0}, 


(3.36) 


m(m  -  1),  , 

1(1*1) 


Ir.  _  ni(m  -  1), 

2‘  1(1  *  1)  ' 


•) 

m 

1(1*1) 


(3.37) 
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The  angular  distribution  of  emitted  energy  is  the  same  for  electric  and  magnetic 
multipole  radiation  of  the  same  1  and  m  17  . 

For  m=0  and  1=  1,  2,  and  3,  plots  of  the  angular  distribution  function  (Z)  are 
shown  in  Figure  3.4,  Figure  3.5,  Figure  3.6,  respectively. 

The  angular  distribution  of  the  second  gamma  in  a  two  step  cascade  in  which 
the  first  gamma  ray  is  emitted  isotropically  is  plotted  in  Figure  3.7.  For  this  plot. 
10.000  isotropic  gammas  were  emitted,  followed  by  a  second  gamma,  distributed  by 
the  angular  distribution  function  (Z).  The  multipole  (1)  was  randomly  selected.  The 
angle  of  the  second  gamma  was  added  to  the  angle  of  the  first  gamma,  in  order  to 
orient  the  second  gamma  with  respect  to  the  0.00°  aixis.  The  gammas  were  grouped 
into  18  angle  bins. 

The  results  of  this  analysis  show  that  the  gammas  emitted  in  a  sample  in  which 
a  large  number  of  nuclei  are  undergoing  gamma- ray  transitions,  the  gammas  will  be 
emitted  in  isotropicsJly  in  the  laboratory  coordinate  system.  This  is  a  result  of  the 
random  orientation  of  the  nuclei  in  the  sample.  The  correlation  is  only  between  the 
two  gamma.s  emitted  relative  to  each  other,  and  therefore  for  the  randomly  oriented 
nuclei,  the  first  gamma  is  not  initially  oriented.  34l 

3.11  Expectations  of  the  Single-Particle  Model 

Based  on  the  single-particle  model,  a  few  expectations  of  some  consequence 
arise.  Briefly,  they  are: 

1.  The  lowest  permitted  multipole  will  dominate  the  gamma-ray  transitions. 

2.  A  comparis.n  of  an  electric  multipole  emission  of  order  L  to  a  magnetic 
multipole  emission  of  order  L  shows  that  the  electric  multipole  emission  is  more 
probable  by  a  factor  of  approximately  10^. 

3.  Emission  of  order  L— 1  radiation  is  less  probable  than  emission  of  order  L 
radiation  by  a  factor  of  approximately  10^. 


Figure  3.4:  Angular  distribution  of  second  gammas  relative  to  the  first 
gamma  (1=1). 


Figure  3.5:  Angular  distribution  of  second  gammas  relative  to  the  first 
gamma  (1=2). 
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Figure  3.7:  Angular  distribution  of  second  gammas  relative  to  the 
first  gamma,  with  the  multipole  (1)  randomly  selected  for 
10,000  isotropic  gammas. 
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These  expectations  are  based  on  the  very  general  approximations  of  the  single 
particle  model,  which  do  not  always  hold  true.  For  example,  it  is  sometimes  seen 
that  an  E2  transition  will  occur  more  frequently  than  an  Ml  transition.  This  does 
not  follow  the  expectations  stated  above.  34^ 


CHAPTER  4 

THE  GAMENT  CODE 

4.1  Introduction  to  the  GAMINT  Code 

A  gamma  cascade  code  (GAMINT  Code)  was  developed  to  examine  the  effect 
cascades  of  different  makeup  have  on  the  response  of  a  capture  tank.  The  GAMINT 
Code  generates  a  gamma-ray  energy  spectrum  for  the  neutron  capture  reaction 
■“U(n,7)^^®U.  The  code  also  determines  the  cascade  gamma-ray  vdelds,  energies, 
order  of  emission,  and  multiplicities  by  Monte  Carlo  techniques.  Included  in  this 
code  is  the  competing  process  of  internal  conversion.  This  code,  for  the  first  time, 
considers  the  internal  conversion  process  in  the  gamma  cascades  produced  in  the 
decay  of 

Electromagnetic  transitions  in  the  continuum  energy  levels,  transitions  from 
the  continuum  energy  levels  to  the  discrete  energy  levels,  and  transitions  in  the  dis¬ 
crete  levels  are  considered.  In  the  code,  electric  and  magnetic  multipole  transitions 
of  type  El  through  E4  and  Ml  through  M4  are  taken  into  account. 

Figure  4.1  through  Figure  4.7  are  flow  charts  of  the  G.AMINT  code.  These 
charts  summarize  the  main  steps  involved  in  the  computer  code  and  are  prepared 
to  aid  the  reader  in  understanding  the  general  strategy  of  the  code. 

4.2  Input  to  the  GAMINT  Code 

The  input  to  the  GAMINT  Code  is  performed  in  two  different  ways.  There 
is  an  external  input  file  and  an  internal  Block  Data  Subroutine.  The  external  file 
contains  those  variables  which  may  require  changes.  The  Block  Data  subroutine 
contains  those  variables  which  remain  fixed.  This  is  done  in  order  to  facilitate  the 
mechanics  of  changing  the  input  to  the  code. 
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Figure  4.1:  Flowchart  for  the  GAMINT  code  initialization. 
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Figure  4-2:  Flowchart  for  the  discrete  level  selection  criteria. 
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Figure  4.3:  Flowchart  for  the  discrete  level  transition  probabilities. 
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Figure  4.4:  Flowchart  for  the  continuum  selection  criteria. 
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Figure  4.5;  Flowchart  for  the  continuum  level  transitions. 
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Figure  4.6:  Flowchart  for  the  selection  of  the  transition. 
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Figure  4.7:  Flowchart  for  the  tallies  performed  in  the  GAMINT  code. 
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Table  4.1:  Input  variables  to  the  GAMINT  code 


Line  Nt  ,nbor 

Input  Variable 

TJ 

NL.NC.NE.NJ 

3 

XE,XJ,XP 

The  initial  input  into  the  GAMINT  code  consists  of  a  three  hne  input  file 
read  into  the  mciin  program  from  the  external  file  GAMINT5.  Line  1  of  the  input 
file  contains  the  job  title  (TJ)  which  may  be  conveniently  changed  to  keep  track 
of  the  current  apphcation.  Line  2  contains  the  number  of  multipoles  used  in  the 
calculation  (NL),  the  number  of  cascades  to  be  executed  (NC),  the  number  of  energy 
level  bands  to  be  used  in  the  energy  continuum  (NE),  and  the  value  of  the  largest 
level  spin  used  in  the  continuum  region  (NJ).  Line  3  contains  the  initial  nucleus 
energy  (XE),  the  initial  nucleus  spin  (XJ),  and  the  initial  nucleus  parity  (XP).  This 
line  i.s  specific  to  the  nucleus  under  consideration  and  facilitates  the  input  of  nuclei 
excited  above  4.803  MeV  neutron  separation  energy.  These  inputs  may  be  easily 
changed  in  order  to  perform  analysis  on  the  cascade  process. 

Table  4.1  summanzes  the  input  vanables  contained  in  the  GAMINT5  input 

file. 

Internal  to  the  GAMINT  Code  is  a  Block  Data  subroutine  which  is  used  to 
initialize  the  %'alues  of  the  variables  labeled  in  the  common  statements  in  the  code. 
This  subroutine  contains  the  following  input. 

•  Atomic  weight  of  the  compound  nucleus:  (.AW) 

•  Level  density  parameter:  (AE) 

•  Level  density  parameter  for  spin:  (SJ) 

•  Number  of  gamma  energy  groups:  (NG) 
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•  Number  of  discrete  levels:  (NVD) 

•  Number  of  values  for  intemai  conversion  coefficient  in  table  for  a  given  shell: 
(NV) 

•  Cutoff  energy  between  the  discrete  levels  and  the  continuum:  (CD) 

•  Upper  energy  of  the  IG-th  group:  EG(IG) 

•  Hindrance  factors  for  a  particular  multipole  transition  of  moment  IL,  for  elec¬ 
tric  transitions  (IP  =  1)  or  magnetic  transitions  (IP  =  2):  HF(IL,IP) 

•  Statistical  factor  for  the  particular  transition  of  moment  IL:  SL(IL) 

•  Spin  of  the  discrete  level-IV:  VJ(IV) 

•  Parity  of  the  discrete  level-FV':  VP(IV) 

•  Energy  of  the  discrete  level-FV:  VE(IV) 

•  Majdmum  transition  energy  in  given  shell-NS  for  intemai  conversion  coefficient 
table:  ETM(NS) 

•  Intemai  conversion  coefficient  for  sheU-NS,  multipole  type-IP.  moment-IL,  and 
I-th  energy  in  table:  AL(NS,IP,IL,I) 

4.3  Details  of  the  GAMINT  Code 

As  shown  in  the  flowchart  shown  in  Figure  4.1,  the  program  begins  by  reading 
in  the  input  file  (GAMINT5).  This  information  is  printed  out  in  the  output  file  in 
order  to  allow  the  user  to  check  the  accuracy  of  the  code  input. 

The  code  then  calculates  the  Weisskopf  factors  for  the  electric  and  magnetic 
dipoles.  These  values  are  stored  in  arrays  RP(1)  and  RP(2),  respectively.  If  RP(1) 
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as  shown  in  Equation  3.15  in  Section  3.3.  The  factor  (RP(IP))  is  the  relative  prob- 
abihty  of  being  an  electric  multipole  (IP=1)  versus  a  magnetic  multipole  (IP=2) 
transition.  The  magnetic  multipole  transition  is  less  probable  than  the  electric  mul¬ 
tipole  of  the  same  order  since  the  above  ratio  is  less  than  1.000. 

The  code  then  initializes  the  common  variables  through  the  use  of  the  Block 
Data  subroutine.  In  the  Block  Data  sub-routine,  the  known  discrete  levels  of 
are  read  into  the  program.  Table  4.2  [161  displays  the  energy  levels  and  the  spin- 
parity  assignments  of  the  excited  nucleus  used  in  this  program.  There  are  25 
well  known  levels  up  to  the  0.95  MeV  cutoff  energy  (CD). 

The  transition  rates  for  gamma-ray  transitions,  based  on  the  Weisskopf  model, 
are  then  determined  from  all  of  the  continuum  bands  and  discrete  levels  down  to 
all  lower  levels.  To  caJctilate  these  rates,  the  transition  factors  (RL(IL))  must  be 
calculated  for  the  all  multipole  moments  (IL).  The  transition  factor  consists  of  those 
terms  of  the  transition  rate  which  the  electric  and  magnetic  multipole  rates  have  in 
common. 


RL(IL)  =  SL(IL)  -  (AW  -  -(.66667  -  IL))  -  ^  («) 

The  transition  rate  for  a  transition  between  two  discrete  levels  is  then  the 
product  of  the  Weisskopf  factor  RP(IP),  the  transition  factor  RL(IL),  the  hindrance 
factor  HF(IL.IP),  and  the  transition  energy  (ET)  raised  to  the  2*IL— 1  power. 


PV(IL,IP,IV)  =  RP(IP)  *  RL(IL)  *  HF{IL,IP)  «  (ET  «  *(2  -  IL  *  1)).  (4.3) 


In  this  calculation,  the  El  transition  rate  is  hindered  by  multiplying  the  calculated 
transition  rate  by  a  hindrance  factor  of  1.0xl0~®  '34j.  The  E2  transition  is  enhanced 
by  multiplying  the  calculated  transition  rate  by  a  hindrance  factor  of  10.0  35 1.  .A.11 
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Table  4.3:  Transition  probabilities  (per  second)  for  discrete  level 
transitions 


Transition  Energy 

El 

Ml 

E2 

M2 

.1 

3.85x10® 

5.60x10'° 

1.08x10^ 

1.35x10" 

.2 

3.08x10^ 

4.48x10“ 

3.48x10® 

4.32x10® 

.3 

1.04x10® 

1.51x10“ 

2.63x10® 

3.28x10® 

.4 

2.47x10® 

3.58x10“ 

1.11x10“ 

1.38x10^ 

.5 

4.82x10® 

7.00x10“ 

3.39x10“ 

4.21x10' 

.6 

8.32x10® 

1,21x10“ 

8.43x10“ 

1.05x10® 

1.32x10® 

1.93x10“ 

1.82x10“ 

2.27x10® 

.8 

1.97x10® 

2.87x10“ 

3.55  X 10“ 

4.42x10® 

.9 

2.81x10® 

4.08x10“ 

6.40x10“ 

7.96x10® 

1.0 

3.85x10® 

5.60x10“ 

1.08x10“ 

1.35x10® 

other  transitions  have  a  hindrance  factor  set  equal  to  1.0.  There  exists  no  evidence 
for  hindrance  factors  of  other  multipolarities.  However,  these  other  multipolarities 
contribute  very  Uttle  to  the  process.  Table  4.3  shows  the  discrete  transition  proba¬ 
bilities  (sec"^)  for  the  El,  ML,  E2,  and  M2  transitions  for  selected  gamma  transition 
energies. 

The  transition  rate  between  a  continuum  level  band  and  a  lower  lying  level 
includes  a  spin  term.  Using  the  back-shifted  fermi  gas  model  in  Equation  2.17, 
Section  2.4  with  the  level  density  parameter  (a)  set  equal  to  9.20212,  the  fictive 
ground  state  position  (A)  set  equal  to  -0.58,  and  the  spin  parameter  set  equal 
to  3.7283,  the  transition  rate  for  the  continuum  levels  down  to  lower  lying  states 
becomes: 


PT(IL.IP.IJ,IE)  =  RP(IP)  *  RL(IL)  -  ET  -  «(2.  -  IL  -  1)  *  DE  *  (2.  -  RJ  *  1.) 

*  HF(IL,IP)  *  EXP(  -RJ  -  (RJ  -  1)/SJ)  (4.4) 


*  00010091  *  EXP(9.7494  *  SQRT(EE  -  .58))/(EE  -  .58)  -  -1.86000. 
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4.4  Treatment  of  Internal  Conversion 

The  internal  conversion  rate  for  each  of  the  possible  transitions  is  also  de¬ 
termined.  Since  internal  conversion  is  considered  an  additional  decay  process,  it 
increases  the  total  decay  rate.  The  total  transition  rate  A  is  the  sum  of  the  transi¬ 
tion  rate  due  to  a  gamma  transition  A.,  and  the  transition  rate  due  to  an  internal 
conversion  transition  A^. 

A  =  A.,  Aj.  (4.5) 

Once  the  transition  rates  PV(IL,IP,IV)  for  transitions  from  discrete  levels  and 
the  transition  rates  PT(IL,IP,IJ,IE)  for  transitions  from  the  continuum  bands  are 
known,  the  internal  conversion  coefficients  (ALT)  are  determined.  The  coefficients 
are  returned  from  the  FIT  subroutine,  and  by  using  the  definition  of  the  internal 
conversion  coeffiaent  (Equation  3.21,  Section  3.6),  together  with  the  known  gamma 
transition  rate,  the  internal  conversion  rate  is  computed. 

In  the  discrete  transitions; 

PE(IL,IP,IV)  =  PV(IL.IP,IV)  "  ALT.  (4.6) 

In  the  continuum  bands  to  lower  level  transitions; 

PC(IL,IP.IJ.IE)  =  PT(IL.IP,IJ,IE)  *  ALT.  (4.7) 

The  internal  conversion  coefficients  are  determined  by  a  procedure  published 
by  Hager  and  Seltzer  ;31].  The  conversion  coefficients  can  be  determined  for  the  four 
lowest  electric  and  magnetic  raultipole  moments  (El,...,E4,  Ml....,M4),  for  electrons 
in  the  K-,  L-,  or  M-shells.  The  table  values  for  internal  conversion  coefficients  are 
shown  in  Figure  4.8  through  Figure  4.10.  These  values  are  entered  into  the  code 
through  the  Block  Data  subroutine.  In  order  to  determine  the  interned  conversion 
coefficient  for  transition  energies  between  any  two  given  values,  a  cubic  sphne  in¬ 
terpolation  subroutine  (SPLINE)  provides  a  piece-wise  polynomial  fit  to  the  data. 
This  insures  that  the  accuracy  of  the  table  is  preserved. 


lb'  lo’ 

TRANSITION  ENERGY  -  KEV 


Figure  4.8:  Internal  conversion  coefficients  for  K-shell  conversion  used 
in  the  GAMINT  Code 
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Figure  4.10:  Internal  conversion  coefficients  for  M-shell  conversion  used 
in  the  GAMINT  Code 
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The  subroutine  FIT  serves  as  an  input-output  controlling  routine  and  this 
subroutine  calls  the  subroutine  SPLINE  in  which  the  actual  curve  interpolation  is 
performed.  In  the  main  program  it  is  necessary  to  determine  the  internal  conversion 
coefficient  for  each  gamma-ray  transition.  .A.  call  is  made  to  Subroutine  FIT  in  a  loop 
over  the  three  electron  shells.  Since  the  table  values  are  in  KeV.  the  transition  energy 
(in  MeV)  is  multiplied  by  a  factor  of  1000.0  and  this  value  is  used  as  the  energy 
interpolation  value.  Each  of  the  three  shells  have  a  different  number  of  known  values 
in  the  table.  Each  shell  has  associated  with  it  a  number  NV(NS)  to  indicate  the 
number  of  table  v<ilues.  When  the  call  to  FIT  is  made,  the  important  parameters 
passed  into  the  subroutine  include  the  electron  shell  (NS),  the  number  of  values 
in  table  NV(NS)  for  shell  NS,  and  the  transition  energy  (ETT).  The  important 
parameter  passed  back  to  the  mtiin  program  through  the  common  statement  is 
the  internal  conversion  coefficient  for  the  given  electron  shell  .\LV(NS).  A  value  is 
obtained  for  the  conversion  coefficient  for  each  of  the  three  shells  (K,L,  and  M)  by 
calling  the  FIT  subroutine  three  times,  once  for  each  shell.  These  three  values  are 
summed  to  give  the  total  interned  conversion  coefficient  a  (ALT). 

ALT  =  ALV(l)  -  ALV(2)  -  ALV(3),  (4.8) 

where  .ALV(NS),  (NS=1,2,3)  represents  the  partial  internal  conversion  coeffiaents. 
The  three  values  for  the  electron  shell  conversion  coefficient  are  summed  in  the  main 
program. 

4.5  Transition  Selection 

The  total  transition  rates  for  all  transitions  are  determined  by  summing  the 
transition  rate  for  gamma  decay  and  the  transition  rate  for  internal  conversion.  For 
a  discrete  to  discrete  level  transition,  the  total  transition  rate  is; 


TT(IL.IP.IV)  =  PVflL.IP.IV)  -  PEflL.IP.IV). 


(4.9) 
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and  for  a  continuum  band  to  lower  level  transition,  the  total  transition  rate  is; 

TU(IL,IP,IJ,IE)  =  PT(IL.IP,IJ,IE)  ^  PT{IL.IP.IJ.IE).  (4.10) 

The  total  transition  rates  are  summed  in  a  temporary  variable  (Tl)  and  each  in¬ 
dividual  rate  is  divided  by  Tl  in  order  to  normalize  the  total  transition  rate  for  a 
particular  transition. 


TT(IL,1P.IV)  =  TT(IL,IPJV)/T1.  (4.11) 

and 

TU(IL,IP,IJ,IE)  =  TU(IL,IP,IJ,IE)/T1 .  (4.12) 

A  cumulative  probability  distribution  (T2)  is  created  by  summing  the  normal¬ 
ized  rates  and  sampled  with  a  random  number  to  select  a  transition. 

T2=  Y.  TT(IL,IP,IV),  (4.13) 

IL.IPJV 

and 

T2  =  Y  TU(IL,IP,IJ,IE).  (4.14) 

IL.IP.IJJE 

A  simple  rejection  technique  is  used  to  select  the  transition.  The  random  number 
is  generated  in  the  GGUBS  subroutine  in  the  IMSLDPLIB  routine  library  on  MTS. 
This  random  number  is  then  compared  to  the  cumulative  probability,  as  the  normzd- 
ized  cumulative  distribution  is  built,  to  select  the  transition  (if  the  random  number  is 
less  than  the  cumulative  probability)  or  reject  the  transition  (if  the  random  number 
is  greater  than  the  cumulative  probability).  When  rejected,  the  transition  selection 
process  is  performed  again. 

Once  the  transition  is  selected,  another  random  number  selection  determines 
if  the  transition  occurs  by  gamma  decay  or  by  internal  conversion.  In  this  selection 
process,  the  gamma  transition  probability  is  divided  by  the  sum  of  the  gamma  and 
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internal  conversion  probabilities.  For  the  discrete  transitions,  the  gamma  transition 
probability  is: 

PG  =  PV(IL.IP,IV)/(PV(IL,IP,IV)  -  PE(ILJP,rV')).  (4.15) 

This  number  represents  the  gamma  transition  probability  for  a  particular  transition 
energy,  given  that  transition  is  selected.  The  probability  for  internal  conversion, 
given  a  particular  transition  occurs,  is  simply: 

PI  =  PE(IL,IP,r\0/(PV(IL,IP,IV)  ^  PE(IL,IP,IV)),  (4.16) 

The  sum  of  the  gamma  transition  probability  and  the  internal  conversion  probability 
is  unity.  The  above  procedure  is  also  performed  for  the  continuous  level  transitions. 
The  gamma  transition  probability,  given  a  transition  occurs,  is: 

PH  =  PT(IL,IP,IJ,IE)/(PT(IL,IP,IJ,IE)  -  PC(IL,IP,IJ,IE)) ,  (4.17) 

and  the  internal  conversion  transition  probability,  given  a  particular  transition  oc¬ 
curs,  is: 


PJ  =  PC(IL,IP,IJ.IE)/(PT(IL,IP.IJ.IE)  -  PC(IL.IP,IJ.IE)).  (4.18) 

The  new  random  number  is  generated  in  the  GGUBS  subroutine  and  compared  to 
the  gamma  transition  probability.  If  the  random  number  is  less  than  the  gamma 
transition  probability,  the  transition  proceeds  by  gamma  emission.  If  not.  internal 
conversion  occurs.  This  procedure  is  continued  until  the  ground  state  is  reached. 

4.6  Tallies  Performed  by  the  GAMINT  Code 

The  GAMINT  code  performs  as  many  iterations  els  programmed  into  the  input 
file.  Dunng  each  of  the  iterations,  information  is  tallied  for  analysis.  The  specific 
items  of  interest  which  the  code  can  keep  track  of  are 
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•  The  energy  group  of  the  transition. 

•  The  type  of  the  transition  (gamma  ray  or  internal  conversion). 

•  The  multipolarity  of  the  transition. 

•  The  total  gamma-ray  energy  for  each  group. 

•  The  order  of  the  gamma-ray  emitted  in  the  cascade. 

This  information  serves  as  the  basis  for  the  analysis  of  the  gamma  cascades  in  the 
decay  of  the  excited  nucleus. 


CHAPTER  5 
RESULTS 


5.1  Results  of  the  GAMINT  Code 

The  GAMINT  Code  was  implemented  with  the  foUowmg  inp'’t  data  in  the 
input  file  GAMINT5. 

•  Job  title  (TJ):  GAMMA  SPEC  CALCULATION  W/I.C. 

•  Number  of  multipoles  used  in  calculation  (NL):  3 

•  Number  of  cascades  executed  (NC):  1000 

•  Number  of  energy  level  bands  to  be  used  in  the  energy  continuum  (NE):  100 

•  Largest  level  spin  used  in  the  continuum  region  (NJ):  10 

•  Initial  nucleus  energy  (XE):  4.8030  <  XE  <  7.3000 

•  Initial  nucleus  spin  (XJ):  0.50000 

•  Initial  nucleus  parity  (XP);  1.0000 

The  number  of  cascades  INC)  in  the  input  file  determines  how  many  decays 
to  ground  state  (program  loops)  will  be  considered.  The  initial  excitation  energy 
(XE)  may  be  changed  to  allow  for  the  generation  of  a  spectrum  based  on  a  different 
initial  excitation  energy.  The  number  of  level  bands  considered  in  the  continuum 
may  be  varied,  however,  100  proved  to  be  a  sufficient  number  of  bands  to  use. 

The  GAMINT  Code  determines  the  energy,  multiplicity,  ajid  order  of  each 
gamma  ray  in  the  gamma  cascades  emitted  in  the  decay  of  the  excited  nucleus. 
The  code  output  also  contains  useful  information  such  as  the  number  of  conversion 
electrons  per  cascade  for  each  of  the  selected  energy  groups.  The  same  input  file 
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was  run  through  the  FGAMA  Code  which  produces  a  spectrum  without  any  consid¬ 
eration  for  internal  conversion  process.  The  spectra  from  the  two  codes  were  drawn 
on  the  same  graph  using  the  DISSPLA  Graphics  Package  on  the  MTS  computer 
system  at  Rensselaer.  This  plot  is  shown  in  Figure  5.1. 

comparison  shows  that  in  the  low  energy  region,  the  number  of  gamma  rays 
emitted  is  greatly  suppressed  when  internal  conversion  is  considered.  This  is  because 
the  low  energy  transitions  m  a.  high-Z  nucleus  are  more  hkely  to  occur  through  the 
internal  conversion  process  than  the  gamma-ray  emission  process.  This  is  expected 
because  the  internal  conversion  coefficient  decreases  rapidly  with  energy.  In  the 
higher  energy  region,  internal  conversion  has  httle  effect  on  the  number  of  gamma- 
rays  produced.  It  is  seen  however  that  the  spectrum  with  internal  conversion  is 
slightly  higher  than  the  spectrum  without  internal  conversion. 

A  gamma- ray  emission  ratio  was  determined  by  dividing  the  gamma- ray  mul- 
tiphcity  by  the  sum  of  the  gamma-ray  multiphcity  and  the  internal  conversion  mul- 
tiphcity.  A  value  of  unity  for  this  ratio  indicates  that  the  internal  conversion  process 
has  no  impact  on  the  number  of  gamma  rays  emitted  in  a  particular  energy  group. 
A  \’alue  for  the  ratio  of  less  than  unity  indicates  that  the  internal  conversion  process 
decreases  the  number  of  gamma  rays  emitted  in  a  particular  energy'  group.  The 
plot,  shown  in  Figure  5.2,  reveals  that  the  lower  energy  gamma-rays  are  suppressed 
as  a  result  of  the  favorable  competition  of  the  internal  conversion  process. 

Figure  5.3  displays  the  fraction  of  total  gamma  ray  energy  released  per  energy 
group.  For  this  plot,  the  gamma  ray  energies  were  summed  into  their  respective 
energy  bins  and  the  result  was  divided  by  the  total  gamma  ray  energy  released. 
This  results  in  a  normalized  spectrum. 

Figure  5.4  displays  the  multipolarity  of  the  gamma  rays  emitted  in  the  decay 
of  The  Ml  transition  is  enhanced  by  a  factor  of  10  in  this  decay.  The  El 

transition  is  hindered  by  a  factor  of  a  rmllion.  This  leads  to  the  majonty  of  the 
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transitions  being  magnetic  dipole  (Ml). 

The  GAMINT  code  allows  one  to  determine  the  order  in  which  the  gamma 
rays  are  emitted  from  the  compound  nucleus.  The  plot  shown  in  Figure  5.5 
was  constructed  showing  the  frequency  with  which  the  n-th  gamma  ray  (where  n  = 
1,2,. ...7)  is  emitted  in  the  cascade  for  a  given  energy.  The  frequency  is  shown  on  a 
per  cascade  basis.  The  first  gamma  ray  emitted  in  a  cascade  is  larger  on  the  average 
than  the  other  gamma  rays,  with  a  peak  around  1.9  MeV.  The  second  gamma  ray 
is  smaller  on  the  average  than  the  first  gamma  ray,  with  a  peak  around  1.3  MeV. 
The  third  gamma  ray  has  a  peak  around  1.2  MeV.  The  fourth  gamma  ray  has  a 
peak  around  0.8  MeV.  For  the  few  cascades  w’lth  more  than  four  gamma  rays  in  the 
cascade,  the  energy  of  these  gamma  rays  is  m  the  0.0  to  0.5  MeV  region. 

Figure  5.6  through  Figure  5.9  are  scatter  plots  of  the  first  gamma  ray  emitted 
in  a  cascade  versus  the  n‘^  gamma  (for  n=2.3,4.5).  .A.s  the  value  of  n  gets  larger, 
there  are  fewer  cascades  with  that  number  of  gamma  rays  in  it.  Therefore,  as  n 
increases,  more  cascades  have  veilues  of  zero  for  the  n-th  gamma  ray  and  there  are 
more  points  lying  on  the  first  gamma  ray  axis,  indicating  a  value  of  zero  for  the  n-th 
gamma  ray.  These  plots  contain  information  on  1000  cascades. 

Figure  5.10  is  a  scatter  plot  of  the  first  gamma  emitted  in  the  cascade  versus 
the  last  gamma  emitted  in  the  same  cascade.  The  lone  point  in  the  upper  nght 
hand  comer  of  the  plot  is  the  superposition  of  the  cascades  with  only  one  gamma 
ray  from  the  4.803  MeV  level  down  to  ground  state.  Therefore  the  first  gamma  ray 
is  the  last  gamma  ray  for  these  cascades. 

A  3-dimensional  plot  was  constructed  to  show  the  gamma-ray  spectmm  for 
different  excitation  energies.  In  Figure  5.11,  the  gamma  ray  .ipectrum  is  displayed 
for  initicd  excitation  energy  from  4.803  MeV,  up  to  7.3000  MeV.  The  cutoff  level  of 
7.3000  .MeV  was  selected  because  of  consideration  for  the  yrast  level  in  the  excited 
nucleus.  For  exatation  energies  above  7.3000  MeV,  particle  emission  will  occur. 
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Figure  5.4:  Gamma-ray  frequency  versus  energy  presented  by  multi¬ 
polarity 
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Figure  5.5:  Gamma- ray  frequency  versus  energy  presented  by  gamma- 
ray  order. 
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Figure  5.7:  Gamma  Ray  Energy  for  the  First  and  Third  Gamma  Ravs 
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Figure  5.8:  Gamma  Ray  Energy  for  the  First  and  Fourth  Gamma  Rays 
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Figure  5.9:  Gamma  Ray  Energy  for  the  First  and  Fifth  Gamma  Rays 
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Figure  5.10:  Gamma  Ray  Energy  for  the  First  and  Last  Gamma  Rays 
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This  regime  is  not  considered  in  this  work.  The  higher  excitation  energy  nuclei 
exhibit  spectrum  shapes  similar  to  the  lower  energy  spectrum,  however  there  is  a 
slight  shift  to  the  higher  energy  gamma  rays.  This  is  to  be  expected  as  a  result  of 
the  increased  initial  excitation  energy 

Figure  5.12  through  Figure  5.15  show  the  angular  distribution  of  the  n‘^ 
gamma  ray  emitted  in  the  decay  of  formed  in  the  thermal  neutron  capture 
reaction  ■^®U(n.7)^^®U.  In  this  analysis,  1000  cascades  were  considered.  The  angu¬ 
lar  distribution  probabihties  of  the  second  gamma  ray  relative  to  the  first  gamma 
ray  was  computed  using  a  Monte  Carlo  program  based  on  the  angular  distribution 
function  shown  in  Equation  3.37  in  Section  3.10. 

Each  of  the  plots  shown  in  Figure  5.12  through  Figure  5.15,  are  based  on  the 
first  gamma  ray  being  emitted  at  a  preferred  direction  of  0.00°.  Figure  5.12  shows 
the  direction  of  che  second  gamma  ray  with  respect  to  the  first  gamma  ray.  Note, 
aince  the  Ml  transition  is  enhanced  in  the  decay  of  the  angular  distribution 
of  the  second  gamma  ray,  relative  to  the  first  gamma  ray  emitted  in  a  preferred 
direction,  is  a  dipole.  We  see,  however,  that  some  of  the  gamma  rays  fall  outside  of 
the  dipole  distribution.  This  is  a  result  of  other  than  dipole  transitions  occurring  in 
the  1000  cascades  examined. 

Figure  5.13  shows  the  third  gamma  ray  angular  distribution,  relative  to  the 
first  gamma  ray  emitted  along  the  preferred  0.00°  axis.  The  probability  in  this  case 
is  determined  by  multiplying  the  probability  of  the  second  gamma  in  the  cascade  by 
the  probability  of  the  third  gamma  in  the  cascade.  This  is  because  the  third  gamma 
ray  direction  and  probabihty  is  dependent  on  the  direction  of  the  second  gamma  ray 
in  the  cascade,  which  in  turn  is  dependent  on  the  first  gamma-ray  direction.  The  first 
gamma  ray  is  assumed  to  be  emitted  in  the  preferred  direction  with  a  probability  of 
unity.  .Again,  due  to  other  than  Ml  transitions,  there  are  some  gamma  rays  falling 
outside  of  a  dipole  emission  signature. 
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As  the  number  n  (n=2.3.4,5)  of  the  gamma  ray  in  the  cascade  increases,  the 
physical  number  of  points  plotted  in  the  distribution  decreases.  This  occurs  because 
the  number  of  gamma  rays  in  the  cascades  varies. 

5.2  Comparison  to  Published  Results 

The  spectrum  generated  by  the  GAMINT  code  is  compared  to  the  calculated 
spectrum  of  Takahashi,  to  the  experimental  spectrum  of  John  et  al,,  and  to  the 
experimental  data  of  Campion,  in  Figure  5.16.  The  spectrum  determined  by  John 
et  al.  is  only  reported  for  energies  above  1.00  MeV.  The  results  of  the  GAMINT  code 
are  presented  in  terms  of  gamma- ray  intensity  in  photons  per  100  capture  events 
per  0.25  MeV  bin.  The  Takahashi  spectrum  is  reported  in  0.25  MeV  bins.  The 
data  of  Campion  is  plotted  for  several  energies  in  the  0.5  MeV  to  0.6  MeV  region. 
Campion’s  data  is  also  given  for  3.58  MeV,  3.66  MeV,  and  4.06  MeV. 

The  similarities  in  the  spectra  occur  in  the  gross  structure.  In  Takahashi’s 
spectrum,  there  is  a  spike  below  0.2  MeV.  In  the  GAMINT  code  spectrum  this 
spike  is  absent,  a  direct  result  of  internal  conversion  competition.  Around  0.3  MeV 
the  spectrum  drops  very  abruptly.  This  absence  of  gamma  rays  in  the  0.3  MeV 
region  is  due  to  the  lack  of  transitions  that  allow  for  the  selection  rule  criteria  for 
energy  transition  of  0.2  to  0.3  MeV.  The  preponderance  of  the  low  energy  gamma 
rays  emitted  from  the  exated  nucleus  are  emitted  in  the  lower  lying  discrete 
level  transitions. 

The  GAMINT  code  predicts  a  much  lower  gamma  ray  intensity  in  the  0  to  1.0 
MeV  region  and  a  higher  gamma  ray  intensity  in  the  1.0  to  1.5  MeV  region  than 
does  TakaJhashi.  The  spectrum  predicted  by  the  GAMINT  code  is  depressed  at  lower 
energies,  as  expected,  because  of  the  strong  competition  from  intemtd  conversion. 
Internal  conversion  was  not  considered  in  the  Takahashi  calculation. 

The  transitions  in  the  continuum  region  contribute  to  a  large  peak  in  the  1.2 
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Figure  5.12:  Angular  distribution  of  second  gamma  relative  to  first 
gamma. 
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Figure  5.14;  Angular  distribution  of  fourth  gamma  relative  to  first 
gamma. 
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Figure  5.15:  Angular  distribution  of  fifth  gamma  relative  to  first 
gamma. 
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MeV  region.  In  both  calculated  spectra,  the  peak  of  the  gamma  ray  spectrum  is 
shifted  to  lower  energies  than  the  experimentally  determined  spectrum  of  John  et 
al.  The  effect  of  the  single  particle  model  on  the  average  energy  of  the  spectrum  is 
to  shift  it  to  lower  energy,  since  the  model  only  considers  vibrational  spectra  and 
neglects  the  rotational  spectra  found  in  the  actual  nucleus.  The  collective  motion, 
in  the  form  of  vibration  and  rotation,  in  a  large  nuclei  such  as  will  lead  to  a 
gamma  rays  of  higher  average  energy,  due  to  the  larger  transitions  in  the  rotational 


bands.  [361 


CHAPTER  6 

SIMULATION  OF  CAPTURE  TANK  RESPONSE  USING 

MICROSHIELD 

Neutron  transmission  and  self-indication  ratio  measurements  have  been  performed 
_2i  for  ~^U{n,~f)  at  the  Rensselaer  Polytechnic  Institute  Linear  Accelerator  and  re¬ 
cently  elsewhere  [3]  ,4i  in  order  to  determine  resonance  self- shielding.  A  preliminary 
analysis  of  the  data  of  Byoun  et  al.  [5]  by  Harris  et  al.  i6i  examined  the  capture 
tank  efficiencies  for  average  s-wave  and  p-wave  neutron  capture.  The  self-indication 
ratio  depends  on  the  gamma  cascade  multiphdties,  as  well  as  average  spectra,  how¬ 
ever.  Because  the  capture  tank  responds  to  the  combined  effect  of  the  gamma  rays 
of  various  energies  from  a  cascade,  as  transported  through  intervening  material,  the 
cascade  structure  has  a  profound  effect  on  the  response.  This  work  for  the  first  time 
analyzes  the  multiple  gammas  in  the  individual  cascades,  the  transport,  and  the 
resulting  energy  deposition  in  the  capture  tank.  The  GAMINT  Code  (gamma  cas¬ 
cade  model)  was  developed  to  simulate  the  deexcitation  of  the  excited  nucleus 
provides  the  gamma  source  for  this  transport  model.  The  GAMINT  Code  predicts 
the  energy  and  order  of  emission  of  each  gamma  in  a  cascade.  This  is  important  in 
the  study  of  the  energy  deposition  for  the  individual  cascades,  since  the  percentage 
of  gamma- ray  energy  reaching  the  capture  tank  is  highly  dependent  on  the  energy 
of  the  individual  gamma  rays  in  the  cascade.  Therefore  the  cascade  makeup  will 
have  a  profound  effect  on  the  energy  deposited  by  the  individual  cascades. 

6.1  Introduction  to  the  Microshield  Program 

The  Microshield  Program  :37j  was  used  to  perform  gamma  transport  calcula¬ 
tions  between  the  capture  sample  and  the  capture  tank.  This  program  wais  selected 
because  of  its  incorporation  of  buildup  factors,  its  capability  to  include  user  defined 
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materials,  and  its  capability  to  define  energy  groups  for  the  gamma- ray  source  spec¬ 
trum.  This  particular  version  of  Microshield  operates  on  an  IBM  personal  computer. 

Shield  materials  determine  the  radiation  attenuation  and  buildup  character¬ 
istics  for  the  Microshield  calculation  of  the  dose  rate  at  a  chosen  point.  Materials 
are  specified  for  each  of  the  shield  regions  for  the  particular  geometrj'  selected.  The 
Microshield  program  provides  a  listing  of  built-in  materials,  or  the  user  may  create 
a  custom  material  for  each  of  the  shield  regions.  The  custom  materials  may  be  en¬ 
tered  by  selecting  one  of  the  built  in  materials  and  entering  a  default  density,  or  by 
entering  a  default  name,  default  density,  and  the  number  of  atoms  of  each  element 
in  the  material. 

The  buildup  factors  which  may  be  selected  are  Taylor  formulation,  geometric 
progression,  Berger  formulation,  or  no  buildup.  The  buildup  factor  corrects  for  the 
contribution  of  photon  scattering  at  the  dose  point.  In  the  Microshield  program, 
the  buildup  factor  for  each  differential  source  volume  is  calculated  during  the  kernel 
integration.  Values  of  the  hne-of-sight  distance  between  the  source  kernel  and  the 
dose  point,  are  calculated  for  each  new  source  point.  The  attenuation  mean  free 
paths  are  determined  and  used  for  the  build  up  factor  calculation.  This  shielding 
mean  free  path  is  used  for  all  of  the  materials  between  the  source  and  the  dose 
point.  Interpolation  for  energy  is  done  for  the  buildup  formula  cEdculations.  A 
linear  interpolation  scheme  is  used  for  the  Taylor  and  Berger  correlations,  and  a 
logarithmic  interpolation  is  used  for  the  geometric  progression  method. 

For  mixed  materials  and  several  shields,  the  combined  attenuation  exponent 
is  determined  by  the  relationship: 

ghieids  materialg 

Mr=  Z!  Ux  Z  (61) 

fc=i  1=1 

where: 


•  /xr  =  number  of  shielding  mean  free  paths. 
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•  [fij p)i  =  energy  dependent  mass  attenuation  coefficient  of  material  i  in  shield 
k. 

•  Pi  =  density  of  material  i  in  shield  k. 

•  t*  =  thickness  of  shield  k  between  source  and  dose  point. 

Microshield  looks  up  values  for  mass  attenuation  coefficients  by  material  type, 
interpolates  for  energy,  multiphes  by  the  individual  material  density,  and  sums  to 
obtain  the  linear  attenuation  coefficient  for  each  of  the  shic.’d  regions.  The  hnear 
attenuation  coefficients  are  then  multiplied  by  the  path  length  through  each  shield 
and  summed  over  all  shields. 

The  custom  material  densities  are  determined  by  an  atomic  weight  ratio  for 
the  atoms  of  a  particular  species  within  the  custom  material. 


Pi  =  Ptmlk  “ 


N,  -  A, 


(6.2) 


where  the  bulk  density  is  the  density  for  the  overall  custom  material,  N,  is  the 
number  density  of  the  i‘^  species,  and  A^  is  the  atomic  weight  of  the  i‘^  species. 

The  effective  atomic  weights  (EAWj)  of  the  custom  materials  is  the  total  molec¬ 
ular  weight  divided  by  the  number  of  atoms. 


EAW,  = 


E.  Ni  *  A. 


(6.3) 


EiN,  ■ 

The  mass  attenuation  coefficient  for  the  custom  material  is  the  sum  of  the 


density- weighted  individual  species. 

PIP  =  Y^[{plp)i  “  {pilpbrdk)]- 


(6.4) 


The  effective  atomic  number  is  used  for  the  interpolation  of  buildup  factors. 
Microshield  uses  an  effective  atomic  number  for  the  custom  materials. 


EAN 


custom 


E^' 


(6.5) 
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where  Z,  is  the  atomic  number  of  the  atoms  of  species  i.  For  mixed  materials,  the 
effective  atomic  number  is  determined  by: 

V  d.-EAN. 

’  .  (6.6) 

where  d,  is  the  density  of  material  in  shield  i,  and  the  subscript  i  refers  to  the  mixed 
materials.  '37^ 

6.2  Description  of  the  Physical  System 

The  physical  system  modelled  in  this  analysis  consists  of  a  1.25  m.  diameter 
capture  tank  filled  with  a  xylene  based  scintillator.  A  neutron-beam  tube  in  the 
tank  contains  the  capture  sample  which  is  a  0.0792  cm.  thick.  7.8  cm.  diameter 
disk  of  .  Surrounding  the  capture  sample  is  4.95  cm.  thick  annulus  of  ®LiH 
and  the  beam  tube  wall. 

6.3  Microshield  Model 

In  order  to  simplify  the  geometry  for  the  physical  set  up  described  above,  the 
cyhndrical  geometry  has  been  converted  to  a  spherical  geometry  by  equating  the 
chord  length  of  the  cylinder  (Cc^/)  to  that  of  the  sphere  (c^ere)-  This  will  give 
equivalent  results  in  the  Monte  Carlo  transport  through  the  two  geometries.  [38] 


Ccy/  =  2  »  t. 


(6.9) 
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For  a  sphere  of  diameter  d.  the  mean  chord  length  is  the  radius.  Therefore  the 
equivalent  cylinder  may  by  approximated  as  a  sphere  with  a  radius  equeil  to  twice 
the  thickness  of  the  cylinder. 

For  this  Microshield  work,  four  cases  are  analyzed.  The  cases  considered  are 

•  Off-resonant  capture  without  internal  conversion 

•  Resonant  capture  without  internal  conversion 

•  Off- resonant  capture  with  internal  conversion 

•  Resonant  capture  with  internal  conversion 

Setup  1:  In  the  off-resonant  neutron  capture  case,  the  captures  are  uniformly 
distributed  throughout  the  cyUnder.  This  cylinder  has  a  thickness  of  0.0790  cm. 
Therefore  the  equivalent  sphere  has  a  radius  of  0.1584  cm.  The  gamma  rays  originate 
uniformly  throughout  the  sphere. 

Setup  2:  In  the  resonant  neutron  capture  case,  the  captures  occur  pnmarily 
in  the  front  face  of  the  cylinder,  since  there  is  a  high  probability  for  capture.  The 
equivalent  sphere  is  made  to  be  the  same  size  as  in  setup  1,  however,  in  order  to 
model  this  situation  on  Microshield,  the  following  procedure  is  performed.  First  the 
gammas  emitted  into  the  disk  from  the  active  region  of  the  face  are  considered.  In 
this  case  the  active  source  region  is  made  much  smaller  and  located  in  the  center 
of  the  sphere.  The  surrounding  material  is  also  and  this  is  equivalent  to  the 
inactive  region  of  the  cylinder.  In  addition,  the  gammas  emitted  away  from  the  face 
are  modelled  by  a  small  sphere.  The  surrounding  material  is  replaced  by  a  v'acuum 
in  order  to  maintain  the  model  integrity.  The  two  answers  cu^e  added,  resulting  in  a 
description  of  the  resonant  capture  case. 
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6.4  Results  of  the  Microshield  Model 

Table  6.1  shows  the  probabilities  of  transport  of  gammas  of  various  energies 
from  capture  sample  to  capture  tank  for  the  case  of  resonant  capture.  Table  6.2 
shows  the  probabilities  of  transport  of  gammas  for  the  case  of  off-resonant  capture. 
Note  that  lower  energy  gamma  rays  are  not  likely  to  reach  the  capture  tank.  Internal 
conversion  diminishes  the  low  energy  gamma  yield  and  results  in  a  total  gamma 
energy  which  differs  from  cascade  to  cascade. 

Figure  6.1  displays  a  plot  for  the  case  of  resonant  capture  without  internal 
conversion  for  100  representative  cascades.  Since  no  internal  conversion  is  considered 
in  this  case,  the  total  available  gamma  energy  for  each  cascade  is  4.803  MeV.  The 
mean  gamma  energy  deposited  by  a  cascade  in  the  capture  tank  is  3.17  MeV. 

Figure  6.2  displays  a  plot  for  the  case  of  off-resonant  capture  without  internal 
conversion.  As  in  the  previous  case,  no  internal  conversion  is  considered  so  the  toted 
available  gamma  energy  for  each  caiscade  is  4.803  MeV.  The  mean  gamma  energy 
deposited  by  a  cascade  in  the  capture  tank  is  3.08  MeV. 

Figure  6  3  displays  a  plot  for  the  case  of  resonant  capture  with  internal  con¬ 
version.  Since  intemad  conversion  is  considered,  the  total  available  gamma  energy 
for  each  cascade  is  less  than  4.803  MeV.  For  the  100  representative  gamma  cascades, 
the  average  available  energy  is  4.56  MeV.  The  mean  gamma  energy  deposited  in  the 
capture  tank  is  3.15  MeV. 

Figure  6.4  displays  a  plot  for  the  case  of  off-resonant  capture  with  internal 
conversion.  .As  in  the  previous  case,  the  total  available  gamma  energy  for  each 
cascade  is  less  than  4.803  MeV.  For  this  case  the  average  available  energy  per  gamma 
cascade  is  also  4.56  MeV.  The  mean  gamma  energy  deposited  in  the  capture  tank 
is  3.07  MeV. 

Table  6.3  displays  the  average  energy  that  reached  the  tank  for  a  cascade  for 
the  cases  considered.  It  is  shown  in  this  analysis  that  for  either  the  resonant  or  the 
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Table  6.2:  Probability  of  transport  of  gammas  of  various  energies 
from  capture  sample  to  capture  tank  for  off  -resonant  cap¬ 
ture 
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Table  6.3:  Energy  Deposited  in  Tank  by  a  Cascade 
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off- resonant  case,  the  average  amount  of  energy  reaching  the  lank  is  only  slightly 
different  between  internal  conversion  and  the  no  internal  conversion  cascades.  There 
is  a  larger  difference  between  the  resonant  and  off  -resonant  cases  if  one  considers  ei¬ 
ther  the  internal  conversion  case  or  the  no  interned  conversion  case.  The  off-resonant 
case  has  a  greater  amount  of  uranium  to  transport  through,  since  the  capture  is  uni¬ 
formly  distributed  throughout  the  disk.  The  result  is  a  slight  decrease  in  the  net 
amount  of  energy  reaching  the  tank.  The  cases  in  which  internal  conversion  is  con¬ 
sidered  have  fewer  low  energy  gamma  rays  to  start  with  than  the  cases  in  which 
no  internal  conversion  is  considered.  However,  this  difference  has  less  effect  on  the 
energy  deposited  in  the  tank  than  does  the  position  of  the  capture  in  the  disk. 

The  spectra  of  the  cascade  energy  reaching  the  tank  for  the  cases  considered 
are  shown  in  Figure  6.5  through  Figure  6.8. 
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CASCADE  NUnBER 


Figure  6.1:  Total  gamma  cascade  energy  deposited  in  capture  tank  for 
100  representative  cascades  (resonant  capture,  no  internal 
conversion) 
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Figure  6.2:  Total  gamma  cascade  energy  deposited  in  capture  tank 
for  100  representative  cascades  (ofF-resonant  capture,  no 
internal  conversion) 
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rigure  6.4:  Total  gamma  cascade  energy  deposited  in  capture  tank  for 
100  representative  cascades  (oiF-resonant  capture.,  internal 
conversion) 
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Figure  6.6:  Spectrum  of  cascade  gamma  energy  deposited  in  capture 
tank  (off-resonant  capture  ) 


RELftTlVE-rREQUENCy 

0.0  0.2  0.4  0.6  O.B 


9 


CHAPTER  7 

DISCUSSION  AND  CONCLUSION 


The  work  done  in  this  thesis  is  the  development  of  a  gamma  ray  cascade  code  for 
the  ““U(n,7)“®'U  reaction  for  the  determination  of  the  response  of  a  capture  tank 
to  gamma  ray  cascades.  The  results  of  this  thesis  may  be  extended  to  capture  tanks 
other  than  the  one  used  at  Rensselaer  Pol)rtechnic  Institute,  as  well  as  other  neutron 
capture  detectors  i39j.  In  this  work,  the  following  has  been  accomphshed. 

(1)  A  gamma  cascade  model  has  been  created  which  predicts  the  energy  of 
the  gamma  rays  emitted  in  a  cascade,  the  order  of  emission  of  the  gamma  rays 
in  a  cascade,  the  multipolarity  of  the  transitions  in  the  decay  of  the  nucleus, 
the  multiphcity  of  the  gamma  rays  in  each  cascade,  the  angular  distribution  of  the 
gamma  rays  relative  to  the  first  emitted  gamma  ray,  and  predicts  the  spectrum  of 
neutron  capture  reactions  in  excited  up  to  the  yrast  level. 

(2)  For  the  first  time,  the  effect  of  internal  conversion  on  the  gamma  cascade 
energy  has  been  analyzed  in  a  cascade  simulation  and  apphed  to  an  apphcations 
problem.  Internal  conversion  dominates  the  low  energy  portion  of  the  spectrum, 
affecting  the  number  of  low  energy  gamma  rays  emitted  from  the  source.  The 
gross  structure  of  the  spectrum  determined  in  this  work  agrees  well  with  other 
calculated  and  experimental  spectra.  However,  due  to  the  inclusion  of  internal 
conversion,  the  low  energy  gamma  transitions  are  significantly  diminished  in  the 
present  work.  The  low  energy  gamma  rays  are  favorably  attenuated  in  the  transport 
of  the  source  spectrum  through  the  material  in  a  capture  tank.  It  is  found  that  the 
energy  deposition  in  the  tank  is  not  contingent  on  whether  or  not  internal  conversion 
is  considered  at  the  source. 

(3)  The  response  of  the  capture  tank  used  in  experiments  at  Rensselaer  has 
been  ancilyzed.  The  cases  examined  are  resonant  capture  with  and  without  internal 
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conversion,  and  off-resonant  capture  with  and  without  internal  conversion.  .4s  dis¬ 
cussed  above,  the  internal  conversion  had  negligible  effect  on  the  average  cascade 
energy  deposited  in  the  tank.  The  off-resonant  capture  occurs  randomly  in  the  cap¬ 
ture  sample,  and  due  to  self-shielding,  the  energy  deposition  in  the  tank  was  slightly 
lower  than  for  resonant  capture.  In  the  resonant  capture  case,  the  gamma  rays  orig¬ 
inate  near  the  ^^ace  of  the  capture  sample,  decreasing  the  amount  of  self- shielding 
in  the  sample.  The  source  spectrum  for  s-wave  and  p-wave  neutrons  was  compared 
and  found  to  be  virtually  the  same.  Therefore  the  position  of  the  capture  reaction 
in  the  sample  has  more  effect  on  the  energy  deposited  in  the  capture  tank  than  the 
other  factors  examined. 

(4)  The  methodology  and  the  findings  of  this  work  are  also  relevant  to  any 
capture  tank  or  other  gamma-ray  detectors.  By  using  the  results  of  the  gamma- 
ray  cascade  transport  calculations,  one  may  determine  an  optimum  discriminator 
level  by  maximizing  the  number  of  cascades  which  have  an  energy  exceeding  the 
discriminator  setting.  This  will  result  in  more  light  collection  for  a  given  number  of 
captures. 

In  order  to  determine  if  the  single  particle  model  is  an  acceptable  model  to  use 
in  describing  the  decay  of  the  excited  nucleus,  one  must  know  how  deformed 

the  nucleus  is.  Examination  of  the  Nilsson  diagram  for  shows  that  the  nucleus 
is  only  slightly  deformed,  and  therefore  the  single  particle  model  gives  an  acceptable 
prediction  of  the  gamma  decay  of  this  nucleus.  However,  because  rotational  bands 
are  actually  present,  the  single  particle  model  will  predict  an  average  gamma-ray  en¬ 
ergy  lower  than  the  collective  model  and  this  is  seen  in  the  shifting  of  the  calculated 
spectrum  to  lower  energy.  A  gajnma  Cciscade  code  based  on  the  collective  model  is 
only  possible  if  the  level  spadngs,  with  the  rotational  bands  included,  are  known. 
This  is  not  the  case,  so  further  experimental  and  theoretical  work  in  determining 
the  collective  model  level  spaangs  would  allow  for  a  more  accurate  cascade  model. 


Future  work  on  this  gamma  cascade  simulation  should  also  incorporate  particle 
emission  for  those  excitation  energies  above  the  yrast  level.  Further  improvements 
may  be  made  by  expanding  the  number  of  nuchdes  which  may  be  examined. 
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APPENDIX  A 

GAMINT  SOURCE  CODE 

C  PROGRAM  GAMINT3  (INPUT. OUTPUT) 

C  •GAMINT  COMPUTES  GAMMA  ENERGY  SPECTRUM  AND  MULTIPLICITY 
C  FROM  NUCLEAR  DEEXCITATION  WITH  INTERNAL  CONVERSION 
C  INCLUDED. 

C  ALSO  COMPUTES  PULSE  HEIGHT  DISTRIBUTION  IN  A  CAPTURE  TANK. 

C 

C 

REAL  El(3) .E2(3) ,AL(3.2.4.23) .ALV(3) .ALT,CG(40) .CI(40) 

1, TJ(20) ,PV(3,2,25) .PT(3.2.40.1000) .SG(40) ,EO(1000,10) 

2, VJ(30) .VP(30) ,VE(30) ,EG(40) ,DP(20) .SL(3) ,CD 

3, RP(2) ,AW,AE.SJ.GE,HF(3,2) ,RL(3) .PS(50) ,PF(5) ,SA(40) 

4, PE(3,2.2S) ,71(3,2,25) .PC(3.2.40,1000) ,GN(10,40) 

5, TU(3,2,40.1000) ,SE(40) .SB(40) ,ETM(3) .TL(IOOO.IO) 

6, TP(1000.10) ,EINT(40) ,EINTN(40) .ETS(IOOO) ,PGE(40) 

7, TH(1000) .ETH(IOOO) .TG(1000,40) ,TI(1000,40) .CM(40,2.3) 
INTEGER  NS.IP,IC,NA,NP,NV(3) .NG.NVD.II 
COMMON/PASS/TJ,PV,PT.SG,EO 

1 , E1 .E2 .NS ,IP , IL , IC .NA, AL ,NP ,ETT ,ALV , ALT .CG ,CI 

2, VJ,VP.VE,EG.DP,SL,NG,NV.CD,NVD 

3, RP,AW.AE,SJ,GE,HF.RL,PS.PF.SA 

4, PE,TT.PC,GN.II.PGE 

5, TU,SE,SB,ETM,TL,TP 

6 ,  EINT , EIHTN , ETS , TH , ETH , TG , TI , CM 
DIMENSION  RN(IO) 

DOUBLE  PRECISION  DSEED 
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C  DATA 

READ  (5,10)  TJ 

WRITE(6.20)  TJ 

READ  (5,12)  NL.NC.WE.NJ 

WRITE(6,22)  NL.NC.NE.NJ 

READ  (5,11)  XE.XJ.XP 

SEED  =7.0 

DSEED=SEED 

WRITE(6.21)  IE, IJ, IP. SEED 

C  DETERMINE  WEISSKOPF  FACTORS  FOR  ELECTRIC  (1)  AND  MAGNETIC  (2) 
C  DIPOLES 

RP(1)=1. 

RP(2)=l./(4.364*AWt<f.666667) 

C  DETERMINE  TRANSITION  FACTORS  FOR  2*IL  POLE 
DO  80  IL=1.NL 

RL ( IL ) »SL ( IL ) « ( AW*  * ( . 666667*  IL ) ) • ( 1 . 82307E2 1 ) / 

1(164, 443**(2*IL+1)) 

80  CONTINUE 
C  DATA  OUTPUT 

WRITE(6,21)  AW.AE.SJ.GE.CD 
WRITE(6,22)  NG.NVD 
WRITE(6,21)  (EG(IG) .IG=1,NG) 

WRITE(6,21)  (DP(IG) .IG=1,NG) 

WRITE(6,21)  (SL(IL) .IL=1,NL) 

WRITE(6,21)  (RL(IL),IL=1.NL) 

WRITE(6,21)  (RP(IP),IP=1,2) 

DO  82  IP=1.2 

82  WRITE(6,21)  (HF(IL,IP) .IL=1 .NL) 

DO  90  rV=l,NVD 


90  WRITE(6,21)  VJ(IV) .VP(IV) .VE(IV) 
WRITE(6,21)  (PF(IU).IU=1,5) 
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WRITE(6.26) 

C  INITIALIZE 

LA=1111 
LB=2222 
LC=3333 
LD=4444 
LE=5555 
ZA=1. 1111111 
ZB=2 . 2222222 
2C=3. 3333333 
ZD=4. 4444444 
ZE»5. 5555555 
11=0 
12=0 
GNX=0 . 0 
SJ=SJ*SJ=2. 
PGE(21)=.80 
PGE(22)=.79 
PGE(23)=.79 
PGE(24)=.78 
PGE(25)=.78 
PGE(26)=.77 
PGE(27)=.76 
PGE(28)=.75 
PGE(29)=.74 
PGE(30)=.73 
PGE(31)=.72 


119 


PGE(32)=.72 
PGE(33)=.71 
PGE(34)=.69 
PGE(35)=.67 
PGE(36)=.63 
PGE(37)=.59 
PGE(38)=.51 
PGE(39)=.39 
PGE(40)=.16 
DO  92  1=1,20 

92  PGE(I)=0.0 
DO  96  INT=1.NG 
IGT=NG-INT+1 
EINT(IGT)=0.0 
EINTN(IGT)=0.0 

96  CONTINUE 

DO  97  IC=1,NC 
ETH(IC)=0.0 
TH(IC)=0.0 
DO  97  IG=1,NG 

97  TG(IC,IG)=0.0 
DO  100  IG=1,NG 
CG(IG)=0.0 
CI(IG)=0.0 
SE(IG)=0.0 
SG(IG)=0.0 

DO  100  11=1,10 
100  GN(II,IG)=0.0 
DO  101  ID=1,50 
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101  PS(ID)=0. 

DO  103  IC=1,NC 
DO  103  11=1,10 
TL(IC,II)=0.0 
TP(IC.II)=0.0 
103  EO(IC,II)=0.0 
EF=0. 

DO  98  IPT=1,2 
DO  98  ILT=1,NL 
DO  98  IGT=1.NG 
98  CM(IGT.IPT.ILT)=0.0 

C  CASCADES 

DO  160  IC=1,NC 
11=1 
CE*IE 
CJ=XJ 
CP=IP 
IT=0 
ED=0. 

C  COMPUTE  RATES.  HEISSKOPF  RATES  WITH  HINDRANCE  FACTORS  HF 
112  T1=0, 

DO  140  IL=1,NL 
DO  140  IP=1.2 

C  TRANSITION  RATES  TO  DISCRETE  LEVELS 
DO  110  IV=1,NVD 
PV(IL,IP,IV)=0. 

PE(IL,IP,IV)=0. 

ALT=0.0 


ETT=0.0 
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TT(IL.IP.IV)=0. 

IFCCE.LE.VEdV))  GO  TO  107 
IF(CP*(-1.)**(IP+IL+1) .NE.VP(IV))  GO  TO  107 
I1=ABS(CJ-VJ(IV)) 

I2=CJ+VJ(IV) 

102  IF(IL.EQ.Il)  GO  TO  104 
11=11+1 

IF (II. GT.  12)  GO  TO  107 
GO  TO  102 

104  P V  ( IL .  IP ,  IV)  =RP  ( IP )  *111  ( IL )  •  (CE- VE ( IV) )  ••  ( 2* IL+ 1 ) 
PV(IL,IP.IV)=PV(IL,IP.IV)*HF(IL,IP) 

IF(IL.EQ.l)  PV(IL,IP.IV)=PV(IL,IP.rV)*(CE-VE(IV)) 

C  WRITE(6,510)IL.IP.IV,PV(IL,IP.IV) 

C  510  FORMATCIH  ,4HIL  =,I4,21.4HIP  =.I4,2X.4HIV  =,I4,2I, 

C  14aPV  =,E12.5) 

C  COMPOTE  THE  INTERNAL  CONVERSION  COEFFICIENTS 
C  INTITIALIZE  THE  VALUE  OF  THE  TOTAL  ICC. 

C  THE  VALUE  FOR  NS  WILL  BE  LOOPED  FROM  1  TO  3,  THE  NUMBER 
C  OF  ATOMIC  SHELLS  FOR  CONVERSION  (I.E.  K.L.M) 

DO  415  NS=1,3 
NP=1 

ETT= (CE-VE( IV) ) •1000 . 0 
IF(ETT.Eq.O.O)GO  TO  107 
C  WRITE(6.413)IL.NS.ETT 

C  413  FORMATdH  ,5HIL  =  ,  14 , 31 . 8HSHELL  =  ,I4.3I.7HT.E.  =  ,E12.5) 
C  NV(NS)  IS  THE  NUMBER  OF  VALUES  IN  TABLE  FOR  A  GIVEN  NS.IP.IL. 
NA=NV(NS) 

ALV(NS)=0.0 


CALL  FIT 
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IF(ALV{NS)  .LT.OJGO  TO  415 
ALT=ALT+ALV(NS) 

C  IF(IC.GT.3)G0  TO  415 

C  WRITE(6,414)IP.IL.NS.ETT.ALV(NS) 

C  414  FORMATClH  .2(13,11) .4HNS  =.I3.2I.4HE  = 

C  l.E12.5,3X,6HICC  =  .E12.5) 

GO  TO  415 

C  417  'rfRITE(6.418)ETT 

C  418  FORMATCIH  .21,4HE  =  .E12.5.1I.4HMEV  .5X,12H0UT  OF  RANGE) 

415  CONTINUE 

C  WRITE (6, 419) ETT, ALT 

C  419  FORMATCIH  .2X,4HE  =  ,E12.5.1X.4HMEV  .5X,12HICC  TOTAL  =  .E12.5) 

C  COMPUTE  THE  TRANSITION  PROBABILITY  FOR  AN 
C  INTERNAL  CONVERSION  EVENT. 

410  PE(IL.IP.n)=PV(IL,IP,IV)*ALT 
C  COMPUTE  THE  TOTAL  TRANSITION  PROBABILITY  (SUM  THE 
C  GAMMA  AND  THE  ICC  PROBABILITIES). 

TT(IL,IP.IV)=PV(IL,IP.IV)+PE(IL,IP,IV) 

C  ACCUMULATE  THE  TOTAL  TRANSITION  PROBABILITY 
T1=T1+TT(IL.IP,IV) 

107  CONTINUE 

C  IFdC.EQ.l)  WRITE(6,511)  IL, IP. IV. II, 12. LA 

C  511  FORMATdfl  .4HIL  =,I4,2X,4HIP  =,I4.2X.4HIV  =14, 2X. 

C  14HI1  =,I4,2X.4HI2  =,I4,5X,I4) 

C  IF(IC.EQ.l)  ’WRrrE(6.512)  PV(IL.IP ,IV) ,PE(IL.IP .IV) . 

C  1ALT,TT(IL,IP.IV).T1 

C  512  FORMATdH  ,4HPV  =,E12.5,2X,4HPE  =.E12.5.2X,5HALT  =. 

C  1E12.5, 

C  12X.4HTT  =,E12.5.2X,4HT1  =,E12.5) 
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110  CONTINUE 
C  CONTINUUM  LEVELS 

C  TRANSITION  RATES  TO  CONTINUUM  LEVELS 
DO  120  IE=1.NE 
DO  120  IJ=1,NJ 
PT(IL,IP.IJ.IE)=0. 

PC(IL.IP.IJ.IE)=0. 

ALT=0 . 0 
ETT=0.0 

TU(IL.IP.IJ.IE)=0. 

IF(CE.LE.CD)  GO  TO  128 

DE=»<CE-CD)/NE 

EE=CD+(IE-.5)»DE 

RJ=IJ-.5 

I1*ABS(CJ-RJ) 

I2=CJ+RJ 

122  IF(IL.EQ.Il)  GO  TO  124 
11*11+1 

IF(I1.GT.I2)  GO  TO  128 
GO  TO  122 

124  PT(IL,IP,IJ,IE)=RP(IP)*RL(IL)*(CE-EE)**(2.*IL+1) 
1*DE*(2.*RJ+1.)*EIP(-RJ*(RJ+1)/SJ) 

2* . 00010091*EXP (9 . 7494*SqRT (EE+ . 58) ) / (EE+ . 58) •• 1 . 86000 
PT(IL.IP.IJ,IE)=PT(IL.IP.IJ,IE)*HF(IL.IP) 

IF(IL.EQ.l)  PT(IL.IP,IJ,IE)=PT(IL.IP,IJ.IE)*(CE-EE) 

C  WRITE(6.520)IL,IP.IJ,IE,PT(IL,IP,IJ.IE) 

C  520  FORMATCIH  ,4flIL  =.I4,2X,4HIP  =,I4,2I.4HIJ  =.I4,2X. 

C  14HIE  =,I4,2I,4HPT  *.E12.5) 

C  COMPUTE  THE  INTERNAL  CONVERSION  COEFFICIENTS 


C  INTITIALIZE  THE  VALUE  OF  THE  TOTAL  ICC. 

C  THE  VALUE  FOR  NS  WILL  BE  LOOPED  FROM  1  TO  3,  THE  NUMBER 
C  OF  ATOMIC  SHELLS  FOR  CONVERSION  (I.E.  K,L,M) 

DO  15  NS*1,3 
NP=1 

ETT=ABS ( (CE-EE) ) * 1000 . 0 
IF(ETT.EQ.O.O)GO  TO  15 
C  WRITE(6.13)IL,NS,ETT 

C  13  FORMATCIH  .5HIL  =  . 14 . 31 , 8HSHELL  =  ,I4.3I.7HT.E.  =  ,E12.5) 

C  NV(NS)  IS  THE  NUMBER  OF  VALUES  IN  TABLE  FOR  A  GIVEN  NS.IP.IL. 
NA=NV(NS) 

ALV(NS)=0.0 

CALL  FIT 

IF(ALV(NS) .LT.O.)GO  TO  15 
ALT=ALT+ALV(NS) 

C  IF(IC.GT.3)G0  TO  15 

C  WRITE(6,14)IP,IL.NS.ETT,ALV(NS) 

C  14  FORMATdH  ,2(13 , 11)  .4HNS  =.I3.21,4HE  = 

C  1,E12.5,3I,6HICC  =  ,E12.5) 

GO  TO  15 

C  17  WRITE(6,18)ETT 

C  18  FORMATdH  ,2X,4HE  =  ,E12.5,1I,4HMEV  ,5I,12H0UT  OF  RANGE) 

15  CONTINUE 

C  WRITE(6,19)ETT,ALT 

C  19  FORMATdH  ,2I,4HE  *  ,E12.5,1I,4HMEV  .5I,12HICC  TOTAL  =  .E12. 
C  COMPUTE  THE  TRANSITION  PROBABILITT  FOR  AN 
C  INTERNAL  CONVERSION  EVENT. 

PC(IL.IP.IJ,IE)=PT(IL,IP.IJ,IE)*ALT 
C  WRITE ( 6 . 530 ) ALT , PT ( IL , IP . IJ , IE) . PC ( IL . IP . IJ . IE) 
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C  530  FORMATdH  , SHALT  =.E12.5,2I.4HPT  =.E12.5,2I.4flPC 
C  1=.E12.5) 

C  COMPUTE  THE  TOTAL  TRANSITION  PROBABILITY  (SUM  THE 
C  GAMMA  AND  THE  ICC  PROBABILITIES). 

TU(IL.IP.IJ.IE)=PT(IL,IP.IJ.IE)+PC(IL.IP.IJ.IE) 

C  ACCUMULATE  THE  TOTAL  TRANSITION  PROBABILITY 
T1=T1+TU(IL,IP.IJ,IE) 

128  CONTINUE 

C  IF(IC.Et).l)  WRITE(6,22)  IL, IP. IE. U. 11.12. LB 

C  IF(IC.EQ.l)  WRITE(6.21)  PT(IL .IP .IJ .IE) .PC(IL . IP . IJ . IE) . 

C  lALT.TUdL.IP.IJ.IE)  .T1 

120  CONTINUE 
140  CONTINUE 

C  NORMALIZE  AND  SELECT  TRANSITION 
C  T1  IS  SUM  OF  RATES.  T2  IS  CUMULATIVE  PROBABILITY 
CALL  GGUBSCDSEED.l.RN) 

T2=0. 

DO  150  IL=1,NL 

DO  150  IP=1.2 

DO  152  IV=1,NVD 

TT(IL.IP.IV)=TT(IL.IP.IV)/T1 

T2=T2+TT(IL.IP,IV) 

C  WRITE(6,540)TT(IL.IP.IV).Tl.T2.RNd).IL.IP.IV 
C  540  F0RMAT(4HTT  =,E12.5.2I.4HT1  = ,E12.5.2I.4flT2  =,E12.5.2I, 
C  lYHRNCl)  =,E12.4.3X,3(I4.2I)) 

IF(RNd)  .LT.T2)  GO  TO  164 
152  CONTINUE 

DO  154  IE=1,NE 


DO  154  IJ=1,NJ 


Ta(IL,IP.IJ,IE)=TU(IL,IP.IJ.IE)/Tl 

T2=T2+TU(IL.IP,IJ,IE) 

C  WRITE(6,541)TU(IL.IP.IJ.IE) .T1,T2,RN(1) .IL,IP,IJ,IE 
C  541  F0RMAT(4HTU  =.E12.5.2I,4HT1  =.E12.5.21.4HT2  =,E12.5, 
C  17HRN(1)  =, £12.5.31,4(14, 21)) 

IF(RN(1) .LE.T2)  GO  TO  162 
154  CONTINUE 
150  CONTINUE 

C  TRANSITION  TO  DISCRETE  LEVEL 
164  ET=CE-VE(IV) 

IT=IT+1 

CE=VE(IV) 

CJ=VJ(IV) 

CP=VP(IV) 

C  COMPUTE  PROBABILITY  OF  GAMMA  TRANSITION  GIVEN 
C  A  TRANSITION  OCCURS. 

IF(TT(IL.IP.IV) .EQ.O.)  GO  TO  220 
PG=PV(IL,IP,IV)/(PV(IL,IP.IV)+PE(IL,IP.IV)) 

C  COMPUTE  PROBABILITY  OF  INTERNAL  CONVERSION 
C  GIVEN  A  TRANSITION  OCCURS. 

PI=PE(IL,IP,IV)/(PV(IL,IP.IV)+PE(IL.IP,IV)) 

C  IF(IC.GT.5)G0  TO  222 

C  WRITE(6.550)RN(1) .PG,PI,IL,IP,IV 

C  550  FORMATClH  ,E12.5,1X,4HPG  =,E12.5,2I,4HPI 
C  1=, £12.5,31,3(14, 21)) 

GO  TO  222 
220  PG=0. 

PI=0. 

222  CALL  GGUBS(DSEED,1.RN) 
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C  IF(IC.EQ.l)  WRITE(6,22)  IC.IT.IL.IF.IV.LC 

C  IF(IC.EQ.l)  WRITE(6.21)  ET.CE.CJ.CP 

C  IF(IC.EQ.l)  WRITE(6,21)  PV(IL,IP.IV) .PE(IL,IP,IV) , 

C  ITT (IL. IP. IV) 

C  IF(IC.EQ.l)  WRITE(6.21)  PG,PI.RN(1) 

IFCPG.GT.RNCD)  GO  TO  165 
IZ=2 

GO  TO  166 
165  IZ=1 

GO  TO  166 

C  TRANSITION  TO  CONTINUUM  LEVEL 
162  ET=CE-(CD+(IE-.5)fDE) 

IT=IT+1 

CE=CE-ET 

CJ=IJ-.5 

CP=CP*(-1.)*»(IP+IL+1) 

C  COMPUTE  PROBABILITY  OF  GAMMA  TRANSITION  GIVEN 
C  A  TRANSITION  OCCURS. 

IF(TU(IL.IP,IJ,IE) .Eq.O.)GO  TO  320 
PH=PT(IL.IP.IJ..IE)/(PT(IL.IP.IJ,IE)+PC(IL,IP,IJ.IE)) 
C  COMPUTE  PROBABILITY  OF  INTERNAL  CONVERSION 
C  GIVEN  A  TRANSITION  OCCURS. 

PJ=PC (IL . IP . IJ . IE) / (PT(IL . IP , IJ . IE) +PC (IL , IP . IJ . IE) ) 
C  IF(1C.GT.5)G0  TO  330 

C  WRITE(6,560)RN(1).PH.PJ,IL.IP,IJ,IE 

C  560  FORMATdH  ,E12.5,1I.4HPH  =.E12.5,2I,4HPJ 
C  1=.E12.5,3I.4(I4.2I)) 

GO  TO  330 


320  PH=0. 


128 


PJ=0. 

330  CALL  GGUBSCDSEED.l.RN) 

C  IF(IC.LE.IO)  WRITE(6.21)  PH.PJ.RN(l) ,ZA 
IFCPH.GT.RNd))  GO  TO  161 
IZ=2 

GO  TO  169 
161  IZ=1 

169  IF(IC.EQ.l)  WRITE(6,22)  IC.IT.IL.IP.U.IE.LD 

C  IF(IC.EQ.l)  WRITE(6.21)  ET.CE.CJ.CP 

C  IF(IC.EQ.l)  WRITE(6.21)  PT(IL,IP.IJ.IE) .PC(IL.IP.IJ,IE) , 

c  itucil.ip.ij.ie) 

C  IFdC.EQ.l)  WRITE(6,21)  PB.PJ.RNd) 

C  IF(IC.Eq.l)  WRITE(6,22)  IZ 

166  CONTINUE 
C  TALLIES 

C  DETERMINE  THE  ENERGY  GROUP  OF  THE  TRANSITION. 

DO  170  11=1, NG 
IG=NG-I1+1 

C  WRITE(6. 570)11, IG,ET,EG(IG) 

C  570  FORMATClH  ,4HI1  =,I4.2I,4HIG  =,I4,2I.4HET  =,£12.5,21, 

C  14HEG  =,E12.5) 

IF  (ET.LT.EG(IG))  GO  TO  172 

170  CONTINUE 
IG=1 

C  SEPERATE  THE  GAMMA  TRANSITIONS  FROM  THE  INTERNAL 
C  CONVERSION  TRANSITIONS. 

172  IF(IZ.EQ.2)  GO  TO  176 
SG(IG)=SG(IG)+1. 


EO(IC,II)=ET 


TL(IC,II)=IL 

TP(IC,II)=IP 

C  Sim  THE  GAMMAS  FOR  EACH  INDIVIDUAL  MULTIPOLARITY 
CM(IG,IP,IL)=CM(IG,IP,IL)+1.0 
C  DETERMINE  THE  TOTAL  GAMMA  ENERGY  FOR  EACH  BIN 
IF(I2.EQ.2)  GO  TO  184 
DO  182  INT=1.NG 
IGT=NG-INT+1 

IF(ET.LE.EG(IGT))  GO  TO  183 

182  CONTINUE 

183  EINT(IGT)=EINT(IGT)+ET 
ETS ( IC) =ETS ( IC) +EO ( IC , II ) 

ETH  ( IC  )  =ETH  ( IC ) +EO  ( IC ,  II )  ■♦PGE  ( IG ) 

WRITE (6.333)IG,E0(IC,II),PGE(IG) 

333  FORMATdH  ,4HIG  *.I4.2X,4HEa  =  ,E12.5 ,21,5HPGE  =,E12.5) 

184  CONTINUE 
CG(IG)=CG(IG)+1. 

TG(IC,IG)=TG(IC,IG)+1 

GN(II,IG)=GN(II,IG)+1 

11=11+1 

TET=TET+ET 

C  WRnE(6,900)IZ.IG,SG(IG)  ,SE(IG)  .ET.TET 
GO  TO  177 

176  SE(IG)=SE(IG)+1. 

CI(IG)=CI(IG)+1. 

TI(IC,IG)=TI(IC,IG)+1 

TET=TET+ET 

C  WRITE(6,900)IZ,IG,SG(IG) ,SE(IG) .ET.TET 
C  900  FORMATdH  ,4HIZ  =.I4.2X.4HIG  =.I4.2I.4HSG  =.E12.5.2X 


130 


C  1.4HSE  =.£12.5,21 

C  2 , iSHTRlNSmON  ENERGY  =.E12.5,2I.14flT0TAL  ENERGY  =.E12.5) 

GO  TO  199 

177  CALL  GGUBS(DSEED,2.RN) 

EU=ET-ET*ET/24. 

DU=EU/5. 

DO  171  ITJ=1.5 

IF(RN(1) .LT.PF(IU))  GO  TO  173 
171  CONTINUE 
IU=5 

173  EH=DU*(IU-1)+RN(2)*DU 
EB=Efl*(.95+.006*ET) 

CR=DP(IG) 

IF(ET.GT.l.Ol)  CR=CR+.0215>»(1.-1./ET) 

EH=EH*CR 

ED*ED+EH 

IF(IC.LE.IOO)  WRITE(6.27)  IC.ED.ET.EU.EH.ZB 
199  IF(CE.GT.VE(1))  GO  TO  112 
C  END  CASCADE  AT  GROUND  STATE 
DO  174  ID=1,25 
IFCED.LT. .2*ID)  GO  TO  175 

174  CONTINUE 
ID=25 

175  PS(ID)=PS(ID)+1. 

EF=EF+ED 

160  CONTINUE 

C  DETERMINE  THE  ORDER  OF  GAMMAS  EMITTED 
C  DO  724  11=1,10 

C  DO  724  13=1, NG 
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C  IG=NG-I3+1 

C  GNX=GNX+GN(II.IG) 

C  WRITE(6. 726)11, IG.GN(II.IG),GNX 

C  726  FORMATdH  .15H0RDER  OF  GAMMl:  .1X.I4,2X, 

C  16HGROOP ; , IX , 14 , 2X . 17HNUMBER  OF 

C  IGAMMAS: .E12.5.4X,E12.5) 

C  724  CONTINUE 

DO  728  11=1,7 
DO  728  I4=1,NG 
IG=NG-I4+1 

GN(II,IG)=GN(II,IG)/NC 
WRrrE(6, 727)11, EG(IG) ,GN(II,IG) 

727  FORMlTClH  ,I4,2X.E12.5,2X,E12.5) 

728  CONTINUE 

C  PRINT  OUT  THE  GAMMAS  EMITTED  IN  GIVEN  CASCADE 
DO  820  IC=1,NC 
WRITE(6,826)IC 

826  FORMAT (14) 

DO  820  11=1,7 

WRITE(6, 827)11, EOCIC, II) ,TP(IC,II) ,TL(IC,II) 

827  FORMAT (14 ,2X ,E12 . 5 , 2X ,E12 . 5 ,2X ,E12 . 5) 

820  CONTINUE 

C  PRINT  OUT  THE  FRACTION  OF  GAMMA  ENERGY  IN  EACH  BIN 
DO  284  INT=1,NG 
IGT=NG-INT+1 

EINTN  ( IGT) =EINT ( IGT) /TET 
284  WRITE(6,901)IGT,EINTN(IGT) ,EG(IGT) 

901  F0RMAT(I4,21,E12.5,2X,E12.5) 

C  PRINT  OUT  THE  INDIVIDUAL  MULTIPOLARITT  OUTPUT 
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DO  950  IP=1.2 
DO  950  IL=1.NL 
DO  950  IG=1,NG 

950  WRITE(6,952)EG(IG) .CMCIG.IP.IL) 

952  F0RMAT(E12.5,21,E12.5) 

C  CASCADE  AVERAGES 

C  THE  STANDARD  DEVIATION  AND  FRACTIONAL  UNCERTAITY  IN 
C  THE  GAMMA  MULTIPLICITY  ARE  DETERMINED. 

TE=0. 

TY=0. 

TZ=0. 

YT=0. 

ZT=0. 

DO  180  IG=1,MG 
SD«0. 

SI=0. 

RD=0. 

RI=0. 

WRITE(6.21)  SG(IG),SD.SE(IG),SI,ZC 
IF(SG(IG) .GT.O.)  SD*1./S0RT(SG(IG)) 
IF(SG(IG).GT.O.)  RD*SD/SG(IG) 

SG(IG)=SG(IG)/MC 

TY=TY+SG(IG) 

IF(SE(IG).GT.O.)  SI=1./SQRT(SE(IG)) 
IF(SE(IG).GT.O.)  RI=SI/SE{IG) 

SE(IG)=SE(IG)/NC 

TZ=TZ+SE(IG) 

IF(IG.LT.NG)  DE=EG(IG)-EG(IG+1) 

IF(IG.EQ.NG)  DE=EG(NG) 
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WRITE(6,21)  DE.SD.RP  SI.RI.ZD 

SA(IG)=SG(IG)/DE 

SB(IG)=SE(IG)/DE 

IF  (IG.LT.NG)  EA=(EG(IG)+EG(IG+1))*.5 
IF  (IG.EQ.NG)  EA=EG(NG)*.5 
EY=EA*SG(IG) 

EZ=EA*SE(IG) 

YT=YT+EY 

ZT=ZT+EZ 

WRITE(6,28)  EG(IG) ,SG(IG) ,SD.SA(IG) .EY.YT.ZE 
WRITE(6,28)  EG(IG) ,SE(IG) .SI.SB(IG) .EZ.TET.ZT 
180  CONTINUE 
EF=EF/NC 
TM=TY+TZ 

WRITE(6,21)  EF,TY,TZ,TM 

DO  185  IG=1,NG 

SG(IG)=SG(IG)/TM 

SA(IG)=SA(IG)/TM 

SE(IG)=SE(IG)/TM 

SB(IG)=SB(IG)/TH 

185  WRITE(6,21)  EG(IG) .SG(IG) .SA(IG) ,SE(IG) .SB(IG) 

DO  190  ID=1.25 
SP=0. 

IF(PS(ID).GT.O.)  SP=1./SQRT(PS(ID)) 

PS(ID)=PS(ID)/NC 

T4=.2*ID 

190  WRITE(6.21)  T4,PS(ID).SP 

C  PRINT  OUT  THE  NUMBER  OF  GAMMAS  AND  INTERNAL  CONVERSIONS  PER 
C  ENERGY  GROUP 
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DO  992  IG=1.NG 

WRITE(6,21)EG(IG) .CG(IG) .CI(IG) 

992  CONTINUE 

DO  986  IC=1,NC 

986  WRITE(6.987)IC,ETS(IC) .ETH(IC) 

987  FQRHAT(I4.2I,E12.5,2I.E12.5) 

C  FORMATS 

10  FORMAT (20A4) 

11  FORMAT (5E12. 5) 

12  FORMAT (1014) 

20  FORMAT (IH  .20A4) 

22  FORMAT (IH  .1014) 

21  FORMATClH  ,6E12.5) 

24  FORMATClH  ,3I4,3E12.5) 

25  FORMATClH  .4I3,3E12.5) 

26  FORMATClH  ) 

27  FORMATClH  ,I4,5E12.5) 

28  FORMATClH  ,7E12.5) 

STOP 

END 

C 

C  THIS  SUBROUTINE  FITS  THE  POINTS  TO  A  CUBIC  SPLINE  FIT. 

SUBROUTINE  FH 

REAL  A(23)  .8(4) .C(23) .D(4,23) .EC23) .F(23) .H(23) .S(23) . 

1T(23) ,1(23) ,YC23) .ST(23) 

COMMON/Cl/A.B.C,D.E.F.H.S,T,I,y,ST 

REAL  El(3)  ,E2C3) .AL(3,2.4,23) ,ALV(3) .ALT,CG(40) .CI(40) 

1. TJ(20) ,PV(3.2.25) ,PTC3,2.40 ,1000) .SG(40) ,£0(1000,10) 

2. VJ(30) ,VPC30) .VE(30) .EG(40) .DP(20) .SL(3) ,CD 
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3, RP(2) ,AW,AE,SJ,GE.flF(3,2) .RL(3) .PS(50) ,PF(5) ,SA(40) 

4, PE(3,2.25) .TT(3.2,25) .PC(3.2.40.1000) ,GN(10,40) 

5, TU(3,2.40,1000) ,SE(40) .SB(40) ,ETM{3) .TL(IOOO.IO) 

6, TP(1000.10) .EINT(40) ,EIirrN(40) .ETS(IOOO) .PGE(40) 

7, TH(1000) .ETHdOOO) .TG(1000.40) ,TI{1000.40) .CM(40,2,3) 

INTEGER  NS,IP,IC,NA,NP,NV(3).NG,NVD.II 

COMMON/P ASS/T J .PV .PT.SG .EO 

1 . El ,E2 .NS , IP . IL , IC .NA , AL ,NP .ETT . ALV . ALT. CG . Cl 

2. VJ.VP.VE,EG.DP.SL.NG.NV.CD.NVD 

3. aP.AW.AE.SJ.GE.HF.RL.PS.PF,SA 

4. PE.TT.PC,GN.II.PGE 

5. TU.SE.SB.ETM.TL.TP 

6 .  EINT . EINTN . ETS . TH . ETH , TG , TI . CM 
DIMENSION  ERG(23.3).EN(23) 

C  THE  FOLLOWING  ENERGY  VALUES  SPAN  THE  RANGE  OF  THE  TABLES. 

DATA  ERG/1.,1.7.3.0.6.2.10.,14. ,20. ,28. ,40. ,53. ,70. 

1.83. . 100. .123. .150. .215. .300. .390. .500. .730. .1000. .1250. .1500. 

2.1. . 1.7.3.0.6.2.10. .14. .20.. 28. .40. .53. .70. .83. .100. 

3.123. . 150. .215. .300. .390. .500. .730. .1000. .1500. .0.0 

4.1. . 2. .4. .8. .15. .25. .40. .52.. 70. .103. .150. .280. .500. 

5,10*0./ 

C  WRITE (6, 29) 

El(l)=116.0 

El(2)=-1.0 

El(3)=-1.0 

E2(l)=116.0 

E2(2)=35.0 

E2(3)=29.0 


ES=E1(NS) 
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ETM(1)=1600.0 
ETM (2) =1500.0 
ETM(3)=500.0 

C  WRITE(6,335)NS.E1(NS).ETT 

C  WRITE(6,337)E2(NS) ,ETM(NS) 

IF(ETT.GT.ETM(NS))GO  TO  51 
IF(ETT.LT.E2(NS))G0  TO  50 
IF (NS . EQ . 1 ) E2 (NS) =E1 (NS) 
E1(NS)=-E1(NS)+1. 
E2(NS)=-E2(NS)+1. 
NE2=-E2(NS)+.S01 
NN=NA 

C  WRITE(6,338)E1(NS) ,E2(NS) ,ME2,NN 

IF(ES.LT.0.)NN=NA+5 
DO  15  1=6, NN 
M=ERG(I,NS)+.01 
J=M+NE2 
IJ=9 

C  WRITE (6, 339) IJ 

C  WRITE(6,336)J,M,I.NS.ERG(I.NS) 

IF(M.LT.150)G0  TO  13 
IF(M.GE.390)G0  TO  11 
J=5*((10fJ+25)/50) 

IJ=10 

C  WRITE(6,339)IJ 

GO  TO  14 

11  IF(M.GE.1000)GO  TO  112 
IJ=11 

C  WRITE(6,339)IJ 
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J=10*((J+5)/10) 

GU  TO  14 

112  J=50*((J+25)/50) 

IJ=12 

C  WRITE(6,339)IJ 
GO  TO  14 

13  J=(10*J+5)/10 
IJ=13 

C  WRITE(6.339)IJ 

14  EN(I)=J 

15  C(I)=EN(I)+E1(NS) 

C  WRITE(6,340)C(I) ,ES 
IF(ES.LT.O.)GO  TO  30 
DO  16  1*1,5 
C(I)=ERG(I,WS) 

C  WRITE(6,341)I,NS,ERG(I,NS) 

16  EN(I)*ERG(I.IIS)-E1(NS) 

DO  17  1=1, NA 

IZZ=17 

C  WRITE(6, 330)1, NA,C(I),EN(I),IZZ 

17  C (I) =ALOG(SqRT ( (1021 . 952+C (I) ) *C(I) ) ) 
DO  18  1=1, NP 

IZZ=18 

W=ETT+E1(NS) 

C  WRITE(6, 331)1, W,IZZ 

18  X(I)=ALOG(SQRT((1021.952+W)*W)) 
H=2*IL+1 

DO  19  1=1, NA 


IZZ=19 
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Z=AL(NS.IP,IL.I)*EN(I)**M 
C  WRITE(6.332)I,M,EN(I).2.IZZ 

19  A(I)=ALQG(AL(NS,IP.IL,I)*EN(I)**M) 

CALL  SPLINE (KA.NP) 

DO  20  1=1, NP 
ZZ=EIP(Y(I))/ETr**M 
C  WRITE (6, 342 )ZZ. NS 

20  ALV(NS)=EXP(Y(I))/ETT**M 
GO  TO  33 

30  DO  31  1=1, NA 

C  WRITECe, 333)1, EN(I+5) .ALCNS.IP.IL.I) .ETT.ES.EN(I) 
C(I)=AL0G(EN(I+5)) 

31  A(I)=ALOG(AL(NS.IP,IL,I)) 

DO  35  1=1, NP 

C  WRITE(6,343)I,ETT 
35  X(I)=ALOG(ETT) 

C  WRITE(6,334) 

CALL  SPLINE(NA,NP) 

ALV(NS)=EIP(Y(I)) 

C  WRITE(6,344)Y(I) 

33  CONTINUE 
L1=0 

IF(ES.LT.0.)L1=5 

M1=NA+L1 

C  WRITE(6,345)ES,NA,L1,M1 
DO  37  1=1, NP 

C  WRITE(6,346)ETT,EN(L1+1) ,EN(M1) 
IF(ETT.LT.EN(L1+1))ALV(NS)=-1. 

37  IF(ETT.GT.EN(M1))ALV(NS)=-1. 


1 
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GO  TO  47 

C  50  WRITE(6,347) 

50  ALV(NS)=AL(NS,IP.IL.l) 

GO  TO  47 

C  51  WRITE(6.348) 

51  ALV(NS)=AL(NS,IP,IL.WV(NS)) 

C  29  FORMAT (IH  .24HBEGINMING  OF  FIT  ROUTINE) 

C  330  FORMATCIH  .3HI  =,I4.2X.4HNA  = .14.21, 6HC(I)  =,E12.5 
C  1.2I,7HEN(I)  =.E12.5,2I.12HLOOP  NUMBER  ,14) 

C  331  FORMATClH  ,3HI  =.I4,3X.3HW  =,E12.5.2X,12HLOOP  NUMBER  .14) 

C  332  FORMATClH  .3HI  =.I4,2X,3HM  =.I4.2X,7HEN(I)  =.E12.5, 

C  121. 3HZ  =.E12.5.2X,12HLQ0P  NUMBER  ,14) 

C  333  FORMATClH  ,3HI  =.I4.2X,4HEN  =,E12.5,2X.4HAL  =,E12.5. 

C  12X.SHETT  *.E12.5.2X,4HES  *.£12.5,2I,4HEN  =,E12.5,2I,7HL00P  31) 

C  334  FORMATClH  ,14HCALLING  SPLINE) 

C  335  FORMATClH  ,4HNS  =,I4,2X,8HElCNS)  =,E12.5,2X5HETT  =,£12.5) 

C  336  FORMATClH  ,3HJ  =,I4,2X,3HM  =,I4.2X,3HI  =.I4.21,4HNS  =,I4, 

C  12X.11HERGCI,NS)  =,E12.5) 

C  337  FORMATClH  ,8HE2CNS)  =.E12.5,2X,9HETMCNS)  =,E12.5) 

C  338  FORMATClH  ,8HElCNS)  = .E12 . 5 ,2X8HE2CNS)  =.E12.5,2I,5HNE2  = 

C  i.I4,2X4HNN  =,I4) 

C  339  FORMATClH  .19HUSED  THE  J  IN  LINE  .14) 

C  340  FORMATClH  .6HCCI)  = ,E12.5 ,2X,4HES  =.£12.5) 

C  341  FORMATClH  ,3HI  =,I4,2X,4HNS  =,I4,2I,11HERGCI.NS)  =,E12.5. 

C  120HAND  IS  SET  =  TO  CCD) 

C  342  FORMATClH  .7HL00P  20,2I.4HZZ  =,E12.5,2I,9H=  ALVCNS) . 

C  12X,4HNS  =,I4) 

C  343  FORMATClH  .7HL00P  35,2I,3HI  = . 14 , 21 , SHETT  =.E12.5) 

C  344  FORMATClH  ,6HYCI)  =,E12.5) 
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C  345  FORMATClH  ,4HES  =,E12.5.2X.4m  =.I4,21,4HL1  =,I4,2I, 

C  14HM1  =,I4) 

C  346  FORMATClH  ,5HETT  =,E12.5,2I.10HEN(Ll+l)  =.E12.5,2I. 

C  ISHENCHl)  =,E12.5) 

C  347  FORMATClH  ,5HAT  50) 

C  348  FORMATClH  ,5HAT  51) 

C  WRITEC6,349)NS,ALVCNS) 

C  349  FORMATClH  .ISHEND  OF  FIT  LOOP .2I.4H1IS  =.I4,2X,9HALVCNS)  =, 
C  1E12.5) 

47  RETURN 
END 

C  THIS  IS  THE  CUBIC  SPLINE  SUROUTINE 
SUBROUTINE  SPLINE CNA.NP) 

REAL  AC23) ,BC4) ,CC23) ,DC4,23) ,EC23) ,FC23) ,HC23) .SC23) , 

1TC23) ,XC23) ,YC23) ,STC23) 

COMMON/Cl/A,B,C,D.E,F,H,S,T,X.y.ST 
C  WRITEC6.338) 

338  FORMATClH  , 27HBEGINNING  OF  SPLINE  ROUTINE) 

M=NA-1 

DO  7  1=1, M 

ECI)=CCI+1)-CCI) 

7  FCI)=ACI+1)-ACI) 

MS=NA-2 

DO  8  1=1, MS 
DC1.I)=ECI+1) 

DC2,I)=2.*CECI)+ECI+1)) 

DC3,I)=ECI) 

8  DC4,I)=3.*CECI+1)*FCI)/ECI)+ECI)*FCI+1)/ECI+1)) 


SC1)=-1. 
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T(1)=2.*F(1)/E(1) 

DO  9  1=2, M 
L=I-1 

S(I)=-(D(1.L)+D(2.L)*S(L))/(D(3.L)*S(L)) 

9  T(I)  =  (-Ta)*(D(2.L)+D(3,L)*S(I))+D(4.L))/D(3.L) 
AA=2.*F(H)/E(H) 

H(NA)=(T(M)+AA*S(M))/(S(H)+1.) 

DO  10  1=1, M 
K=NA-I 

10  H(K)=(H(K+1)-T(K))/S(K) 

DO  12  1=1, NP 
DO  13  K=1,M 

IF(X(I).LT.C(K+1))  GO  TO  14 

13  CONTINUE 

14  K1=K 

Z*(C(Kl+l)-C(Kl))/2. 

B(3)=(H(K1+1)-H(K1))/(4.*Z) 

B (4) = . 25* ( (H(K1+1) +H(K1) ) *Z-A(K1+1) +A (Kl) ) /Z**3 

B(l)=(A(Kl+l)+A(Kl)-2.*B(3)*Z**2)/2. 

B(2)=(H(Kl+l)+H(Kl)-6.*B(4)*Z**2)/2. 

ZP=X(I)-(C(Kl+l)+C(Kl))/2. 

Y(I)=BCl)+ZP*(B(2)+ZP*(B(3)+ZP*B(4))) 

12  CONTINUE 
RETURN 
END 

BLOCK  DATA 

REAL  El(3)  .E2(3) .AL(3,2,4,23) .ALV(3) ,ALT.CG(40) ,CI(40) 

1, TJ(20) ,PV(3,2.25) ,PTC3.2.40,1000) .SG(40) .£0(1000,10) 

2, VJ(30) ,VP(30) .VE(30) ,EG(40) .DP(20) .SL(3) .CD 


3, RP(2),AW.AE.SJ,GE,HF(3.2).RL(3),PS(50),PF(5) .SA(40) 

4, PE(3.2,25) ,TT(3.2.25) .PC(3.2,40,1000) ,GN(10,40) 

5 , TU(3 ,2 ,40 . 1000) ,SE(40) ,SB(40) ,ETM(3) ,TL(1000 . 10) 

6, TP(1000.10)  .EIirr(40)  ,EmTN{40)  .ETS(IOOO)  .PGE(40) 

7,  TH(1000),ETH(1000).TG(1000, 40), 11(1000,40) ,CM(40,2.3) 
INTEGER  NS.IP,IC,NA.NP.NV(3),NG,NVD.II 
COMMON/PASS/TJ , PV , PT . SG ,E0 

1, El,E2,NS,IP.IL,IC.NA,AL.NP.ETr.ALV.ALT.CG,CI 

2, VJ,VP,VE,EG.DP,SL,NG,NV,CD,NVD 

3, RP,AW,AE,SJ,GE,HF,RL,PS,PF,SA 

4, PE,TT,PC,GN,II,PGE 

5, TU,SE,SB,ETM.TL,TP 

6 ,  EINT , EINTN , ETS , TH , ETH , TG , TI ,  CM 

DATA  AW,AE,SJ,GE,CD/239.2,9.302,3.702, .04948, .95/ 

DATA  NG,NVD/40,25/ 

DATA  NV(1),NV(2),NV(3)/23,17,8/ 

DATA  EG(1) ,EG(2) ,EG(3)/10. ,9.75,9.50/ 

DATA  EG(4) ,EG(5) ,EG(6)/9.25,9.0,8.75/ 

DATA  EG(7) ,EG(8) ,EG(9)/8 .5,8.25,8.0/ 

DATA  EG(IO) ,EG(11) ,EG(12)/7.75,7.5.7.25/ 

DATA  EG(13) ,EG(14) ,EG(15)/7.0,6.75,6.50/ 

DATA  EG(16) ,EG(17) ,EG(18)/6.25,6.0,5.75/ 

DATA  EG(19) ,EG(20) ,EG(21)/5.5,5.25,5.0/ 

DATA  EG(22) ,EG(23) ,EG(24)/4.75,4.5,4.25/ 

DATA  EG<25) ,EG(26) ,EG(27)/4. ,3.75,3.50/ 

DATA  EG(28) ,EG(29) ,EG(30)/3.25.3.0,2.75/ 

DATA  EG(31) ,EG(32) ,EG(33)/2.5,2.25,2.0/ 

DATA  EG(34) ,EG(35) ,EG(36)/1.75,1 .5,1 .25/ 

DATA  EG(37) ,EG{38) .EG(39)/1 . ,0.75,0 .5/ 
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DATA  EG (40) /O. 25/ 

DATA  DP(1) .DP(2) .DP (3)/ .533, .565, .602/ 

DATA  DP(4) ,DP(5) .DP(6)/ .646, .697, .746/ 

DATA  DP(7) .DP(8) .DP(9)/.782, .801, .808/ 

DATA  DP (10) .DP (11) .DP (12)/. 800, .780, .752/ 
DATA  DP (13) .DP (14) .DP (15)/. 677, .549, .341/ 
DATA  DP(16) .DP(17) .DP(18)/.122. .004, .000/ 
DATA  HF(l.l) .HF(2.1) .HF(3.1)/. 000001,1. .1./ 
DATA  HF(1.2),HF(2.2).HF(3,2)/10.,1.,1./ 

DATA  SL(1) .SL(2) ,SL(3)/ .25 . .0048, .000060469/ 
DATA  VJ(1).VP(1),VE(1)/2.5.1.,.0/ 

DATA  VJ(2).VP(2).VE(2)/3.5.1.,.043/ 

DATA  VJ(3) ,VP(3) ,VE(3)/4.5.1. , .097/ 

DATA  VJ(4),VP(4),VE(4)/.5,1.,.133/ 

DATA  VJ(5),7P(S),VE(5)/1.5,1.,.144/ 

DATA  VJ(6) ,VP(6) ,VE(6)/5.5,1. , .165/ 

DATA  VJ(7) ,VP(7) ,VE(7)/3.5,1. , .173/ 

DATA  VJ(8) ,VP(8) .VE(8)/2.5,1. , .191/ 

DATA  VJ(9) ,VP(9) ,VE(9)/3.5,1. . .220/ 

DATA  VJ(IO) ,VP(10) ,VE(10)/4.5,1. . .229/ 

DATA  VJ(ll) .VP(ll) ,VE(11)/5.5,1. , .300/ 

DATA  VJ(12) ,VP(12) ,VE(12)/4.5,1. . .307/ 

DATA  VJ(13) .VP(13) ,VE(13)/6.5,1. , .373/ 

DATA  VJ(14) ,VP(14) ,VE(14)/ .5.1. , .696/ 

DATA  VJ(15) ,VP(15) ,VE(15)/1.5.1. , .715/ 

DATA  VJ(16) .VP(16) ,VE(16)/1.5,1. . .735/ 

DATA  VJ(17) ,VP(17) .VE(17)/1 .5,-1. . .742/ 

DATA  VJ(18) ,VP(18) ,VE(18)/.5.-l. , .754/ 

DATA  VJ(19) ,VP(19) ,VE(19)/3.5,-l. , .799/ 


DATA  VJ(20) ,VP(20) .VE(20)/.5.-l. . .812/ 

DATA  VJ(21) ,VP(21) , VE(21)/1 .5.-1 . . .821/ 

DATA  VJ(22) ,VP(22) ,VE(22)/2.5.-l. . .854/ 

DATA  VJ(23) ,VP(23) .VE(23)/5.5.-l. . .894/ 

DATA  VJ(24) .VP(24) ,VE(24)/2.5.-l. , .930/ 

DATA  VJ(25) .VP (25) .VE(25)/4.5.-l. . .943/ 

DATA  PF(1).PF(2).PF(3)/.01..04..1/ 

DATA  PF(4).PF(5)/.25.1./ 

DATA  ETM(1),ETM(2).ETM(3)/1600. 0,1500. 0,500.0/ 

DATA  AL(l,l,l,l),AL(l,l,l,2)/.273,.269/ 

DATA  Aid. 1,1, 3) ,AL(l.l,l,4),AL(l,l,l,5)/.262, .247, .230/ 

DATA  AL(1, 1,1, 6), AL(1, 1,1, 7), Aid,  1,1, 8)/. 212,  .191,  .167/ 

DATA  ALd,l,l,9)  ,ALd ,  1 ,1 ,10)  ,AL(1 ,1 , 1 ,11) /  .  139 ,  .  115 ,  .0923/ 
DATA  ALd,l,1.12)  .ALd,l,l,13)  ,ALd,  1,1, 14)/. 0792,  .0655,  .0521/ 
DATA  ALd,l,l,15)  ,ALd,l,1.16)  ,ALd,l,l,17)/.04l3,  .0256,  .0159/ 
DATA  ALd,l,l,18),ALd,l,l,19),ALd,l,l,20)/.0105,  .00725,. 0040 
18/ 

DATA  ALd,l,l,21)  ,ALd  ,1 ,1 ,22)  ,ALd ,  1 , 1 .23)/ . 00260 ,. 00184 .. 001 
139/ 

DATA  ALd.l,2,l),ALd,1.2,2)/.188..191/ 

DATA  ALd.l,2.3).ALd,l,2.4),ALd,l,2,5)/.198.  .210,  .220/ 

DATA  ALd.l,2,6),ALd,l,2,7),ALd.l,2.8)/.228,.231..229/ 

DATA  ALd,l,2,9),ALd,1.2,10)  ,AL(1,1,2,11)/.217.  .200,  .176/ 

DATA  ALd,1.2,12)  ,ALd,1.2.13)  ,AL(1,1,2,14)/.158.  .138,  .116/ 
DATA  ALd,l,2,15).ALd,1.2.16),ALd.l,2,17)/.096,  .0636,. 0413/ 
DATA  ALd,1.2,18),ALd.l,2.19),ALd,l,2,20)/.0282..0198,  .0113/ 
DATA  ALd.l,2,21).ALd,l,2,22),ALd.l,2.23)/.00717..00501. 
1.00317/ 

DATA  ALCl, 1,3,1) .ALCl ,1 ,3,2)/ .0369, .0441/ 


DATA  AL(1.1.3,3) ,AL(1.1.3,4) .AL(1. 1.3. 5)/. 058, .0927, .133/ 

DATA  AL(1,1,3,6) ,AL(1 .1 .3.7) .AL(1 ,1 ,3,8)/ . 177 , .229, .282/ 

DATA  AL(1,1,3.9) ,AL(1.1.3,10) ,AL(1.1.3.11)/ .332, .354, .353/ 

DATA  AL(1,1,3.12) .AL(1 .1 .3,13) .AL(1 .1,3,14)/ .340 , .316 , .279/ 
DATA  AL(1.1,3,15) .AL(1,1,3,16) .AL(1. 1.3. 17)/. 240, .164, .106/ 
DATA  AL(1.1,3.18) ,AL(1.1,3,19) .AL(1, 1.3, 20)/ .0708, .0480, .0257/ 
DATA  AL(1.1.3,21) ,AL(1.1.3.22) ,AL(1. 1,3, 23)/. 0156, .0105, .00759 
1/ 

DATA  AL(1, 1,4,1) ,AL(1 .1.4.2)/ .00166, .00352/ 

DATA  AL(1.1,4.3) .AL(1,1,4,4) .AL(1, 1,4, 5)/. 00847. .0289, .0668/ 
DATA  AL(1.1.4,6),AL(1,1.4,7),AL(1,1,4.8)/.125,.218..348/ 

DATA  AL(1.1,4,9) .AL(1.1.4,10) .AL(1,1.4,11)/.519. .649, .736/ 

DATA  AL(1,1,4.12) ,AL(1 ,1 .4.13) ,AL(l,1.4.14)/.755, .741, .685/ 
DATA  AL(1,1,4,15) ,AL(1 ,1 .4,16) ,AL(l,1.4,17)/.603, .418, .262/ 
DATA  AL<1, 1,4, 18), AL(1, 1,4, 19), Aid.  1,4, 20)/.  166,.  107,. 0526/ 
DATA  AL(1.1.4,21) ,AL(1 ,1 ,4.22) . AL(1, 1,4, 23)/ .0298, .0191 , .0134/ 
DATA  AL(1.2,1.1) .AL(1,2,1,2)/12.1,11.9/ 

DATA  Aid, 2, 1.3)  .ALd,2,l,4)  ,ALd,2.1,5)/11.6,10.7,9.8/ 

DATA  ALd.2.1,6)  .ALd,2.1,7)  ,ALd.2,l,8)/8.88.7.80.6.63/ 

DATA  ALd.2,1.9).ALd,2,l,10),ALd.2,l.ll)/5.28.4.21,3.21/ 

DATA  ALd.2,1.12),ALd,2,1.13).ALd,2,1.14)/2.66,2.11.1.59/ 
DATA  Aid, 2. 1,15)  .Aid, 2, 1,16)  .Aid. 2,1.17)/!.  19,. 653.  .350/ 
DATA  Aid, 2,1, 18)  .ALCl, 2, 1,19)  .Aid, 2. 1,20)/. 201,.  119,  .0519/ 
DATA  Aid. 2. 1,21), ALCl, 2, 1,22)  .ALCl, 2, 1,23)/.  0264,.  0154,  .0098 
14/ 

DATA  ALCl, 2. 2,1) .ALCl ,2,2.2)/55.9,54.7/ 

DATA  ALd,2.2,3)  .ALd,2,2.4)  ,ALd.2.2.5)/52.8.48.3.43.5/ 

DATA  ALCl, 2. 2. 6)  .ALd,2.2.7)  .ALd,2.2.8)/38.8.33.4.27.6/ 

DATA  ALCl, 2. 2, 9) . ALCl, 2, 2. 10)  .ATd. 2. 2. 11) /21. 1,16. 1,11. 7/ 
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DATA  AL(1,2.2,12).AL(1.2.2,13),AL(1.2.2.14)/9.39,7.17,5.15/ 
DATA  AL(1,2.2,15) ,AL(1.2.2,16) ,AL(1.2,2.17)/3.96,1.84. .907/ 
DATA  AL(1.2,2.18) .AL(1.2.2,19) ,AL(l.2.2,20)/.490. .278. .115/ 
DATA  AL(1.2.2.21) .AL(1.2,2,22) .AL(1.2.2.23)/.0576, .0338, .0219/ 
DATA  AL(1.2.3,1) .AL(l,2.3.2)/52.7.53.2/ 

DATA  AL(1,2,3,3).AL(1,2,3.4).AL{1,2.3.5)/54.2.S5.5.55.7/ 

DATA  AL(l,2,3.6),AL(1.2.3.7).AL(1.2.3.8)/54.7,52.2,47.7/ 

DATA  Aid, 2. 3, 9) , AL(1 .2 .3 , 10)  .AL(1 ,2 ,3. ll)/40 . 5 .33 .3 .25 . 6/ 

DATA  AL(1.2.3.12),AL(1,2,3.13).AL(1.2.3.14)/21.0.16.3.11.9/ 
DATA  AL(1.2.3.15),AL(1,2,3.16).AL(1,2.3.17)/8.49.4.12,1.93/ 
DATA  AL(l,2.3.18).AL(l,2.3.19),AL(1.2,3.20)/.989..533..205/ 
DATA  AL(1,2,3.21) ,AL(1,2.3.22) .AL(1. 2. 3. 23)/ .0974. .0553. .0351/ 
DATA  AL(1,2.4.1) ,AL(1 .2 .4.2) /II .5,13 .6/ 

DATA  AL(1,2.4.3) ,AL(1.2.4,4) ,AL(1,2,4.5)/17.6,26.9.36.8/ 

DATA  AL(1,2,4,6),AI(1.2,4.7).AL(1,2,4,8)/4€.1,56.2,61.8/ 

DATA  AL(1,2,4,9) .AL(1,2,4,10) ,AL(l,2,4,ll;/63.5,59.3,50.4/ 

DATA  AL(1,2,4,12) ,AL(1,2,4,13) ,AL(l,2,4,14)/43.3,35.1.26.2/ 
DATA  AL(1.2,4.15).AL(1,2.4.16),AL{1,2.4.17)/19.0,9.09,4.08/ 
DATA  AL(1,2,4,18) ,AL(1 ,2 .4 ,19) ,AL(1 ,2,4.20)/l . 98 , 1 . 01 , .355/ 
DATA  AL(1,2.4,21) ,AL(1 ,2,4.22) .AL(1 ,2,4,23) /. 158 ,. 0862 ,. 0530/ 
DATA  AL(2,1.1.1).AL(2,1.1.2)/1.51..997/ 

DATA  AL(2,1.1,3).AL(2.1,1.4).AL(2,1,1.5)/.624..350..210/ 

DATA  AL(2,1.1,6),AL(2,1.1,7),AL(2,1.1.8)/.122. .0859, .0579/ 

DATA  AL(2.1,1.9).AL(2.1,1.10).AL(2,l.l.ll)/.0369..0242. .0109/ 
DATA  AL(2,1,1,12) , AL(2. 1,1, 13) .AL(2. 1,1, 14)/. 00525, .00302, .00 
1182/ 

DATA  AL(2,1,1.15) , AL(2, 1,1, 16) .AL(2. 1,1, 17)/. 000885, .000519, 

1 . 000255/ 

DATA  AL(2, 1,2,1) ,AL(2,1.2,2)/1530. 0.706.0/ 


DATA  AL(2.1.2,3) ,AL(2.1,2,4) , AL(2. 1,2, 5)/296. 0,103. 0,40. 8/ 
DATA  AL(2,1,2,6) ,AL(2,1,2,7) ,AL(2,1,2,8)/15.3.8.22,4.09/ 

DATA  AL(2,1,2,9) ,AL(2,1,2,10) ,AL(2,1,2.11)/1.86, .890, .224/ 

DATA  AL(2.1,2,12) ,AL(2, 1,2,13) ,AL(2, 1,2, 14)/. 0661, .0269, .0123/ 
DATA  AL(2.1,2,15) ,AL(2,1,2,16) ,AL(2, 1,2, 17)/. 00428, .00290, 
1.000851/ 

DATA  AL(2. 1,3,1) ,AL(2,1, 3, 2)/95600. 0,38500.0/ 

DATA  AL(2,1,3.3) ,AL(2,1,3,4) , AL(2,1, 3, 5)/13900. 0,3980. 0,1330. 
10/ 

DATA  AL(2.1,3,6) ,AL(2,1,3,7) , AL{2, 1,3, 8)/418. 0,199. 0,86. 6/ 

DATA  AL(2.1,3,9) ,AL(2, 1 ,3,10) ,AL(2,1,3,11)/33.7,13.9,2.62/ 

DATA  AL(2.1,3,12) ,AL(2,1,3,13) ,AL(2,l,3,14)/.580, .188, .0690/ 
DATA  AL(2,1,3,15) ,AL(2,1,3,16) ,AL(2,1,3,17)/.0175, .00691, .002 
125/ 

DATA  AL(2, 1, 4, 1).AL(2,1, 4, 2)/2120000. 0,839000.0/ 

DATA  AL(2,1,4,3) ,AL(2,1,4,4) ,AL(2, 1,4, 5)/286000. 0,74500.0, 
122500.0/ 

DATA  AL(2,1,4,6) ,AL(2,1,4,7) , AL(2, 1,4, 8) /6220. 0,2720. 0,1070.0/ 
DATA  AL(2,1,4,9) , AL(2 , 1 ,4 , 10) ,AL(2 , 1 ,4 , 11) /368 .0 , 135 .0 ,20 .2/ 
DATA  AL(2,1,4,12) ,AL(2,1,4,13) ,AL(2,l,4,14)/3.85, .968, .300/ 
DATA  AL(2,1,4,15) ,AL(2,1,4,16) ,AL(2, 1,4, 17)/. 0588, .0193, .00504 
1/ 

DATA  AL(2,2,1,1) ,AL(2,2,l,2)/80.5,50.4/ 

DATA  AL(2,2,1,3) ,AL(2,2,1,4) ,AL(2,2,1,5)/29.8,15.7,8.90/ 

DATA  AL(2,2,1,6) ,AL(2,2,1,7) ,AL(2,2,1,8)/4.88,3.31,2.14/ 

DATA  AL(2,2,1,9) ,AL(2,2,1,10) ,AL(2,2,1.11)/1.30, .810, .325/ 

DATA  AL(2,2,1,12),AL(2,2,1,13) , AL(2, 2,1, 14)/ . 138 , .0697 , .0366/ 
DATA  AL(2,2,1,15),AL(2,2,1,16) ,AL(2, 2, 1,17)/. 0137, .00639, 
1.00219/ 


148 


DATA  AL(2,2.2.1).AL(2.2,2,2)/6090.0.2930.0/ 

DATA  AL(2,2,2,3) .AL(2.2.2,4) .AL{2, 2.2,5) /1290. 0,483. 0,205.0/ 
DATA  AL(2,2,2,6),AL(2.2.2,7).AL(2.2,2,8)/83.9,47.6,25.3/ 
data  AL(2,2,2,9),AL(2,2,2,10),AL(2,2,2,11)/12.4,6.41,1.85/ 

DATA  AL(2,2,2.12) ,AL(2,2,2,13) .AL(2, 2, 2, 14)/. 603, .256, .117/ 
DATA  AL(2,2,2,15) ,AL(2,2,2,16) ,AL(2,2,2,17)/.0375, .0162, .00527 
1/ 

data  AL(2,2,3,1),AL(2,2,3,2)/547000.0,200000.0/ 
data  AL(2,2,3,3) , AL(2, 2, 3, 4) ,AL(2, 2,3, 5)/65600. 0,17100. 0,5390 
1.0/ 

DATA  AL{2.2,3,6) ,AL(2,2,3,7) ,AL(2, 2,3, 8)/1620. 0,759. 0,328.0/ 
DATA  AL(2,2,3,9),AL(2,2,3,10),AL(2,2,3,11)/128.0,54.1,10.9/ 
DATA  AL(2,2,3.12) ,AL(2,2,3,13) ,AL(2,2,3,14)/2.63, .910, .349/ 
DATA  AL(2,2,3,15) ,AL(2,2,3,16) ,AL(2,2,3,17)/ .0897 , .0339, .00 
1966/ 

DATA  AL(2,2,4,1) ,AL(2. 2, 4, 2)/100000000. 0,23700000.0/ 

DATA  AL(2,2,4,3) ,AL(2,2,4,4) , AL(2, 2, 4, 5) /2040000. 0,426000.0 

1,110000.0/ 

DATA  AL(2,2.4,6) . AL(2 , 2 .4 ,7) . AL(2,2,4,8) /26600 . 0 , 10900 . 0 .4010 

1.0/ 

DATA  AL(2,2,4,9) ,AL(2,2,4,10) .AL(2,2.4,11)/1310.0,466.0,67.9/ 
DATA  AL(2,2,4,12) .AL(2,2,4,13) ,AL(2,2,4.14)/12.3,3.42,1.09/ 
DATA  AL(2.2,4,15) ,AL(2,2.4.16) ,AL(2, 2, 4, 17)/. 219, .0707, .0170/ 
DATA  AL(3,l,l,l),AL(3,l,l,2)/.630,.206/ 

DATA  AL(3,1,1,3) ,AL(3,1,1,4) ,AL(3,1,1,5)/.108, .0515, .0194/ 

DATA  AL(3,1,1,6) ,AL(3 , 1 ,1 ,7) .AL(3, 1,1 ,8)/ .00743, .00165 , .00047 
13/ 

DATA  AL(3,1,2,1),AL(3,1,2,2)/1060.0,138.0/ 

DATA  AL(3,1,2,3),AL(3,1,2,4),AL(3.1.2.5)/42.9,11.3,2.0/ 


DATA  AL(3,1,2,6) ,AL(3,1.2.7) ,AL(3, 1,2, 8)/. 370, .0280, .00365/ 
DATA  AL(3, 1, 3, 1),AL(3, 1,3,2)7120000. 0,9380.0/ 

DATA  AL(3. 1, 3, 3) , AL(3, 1,3, 4) ,AL(3, 1,3, 5)72220.0,432. 0,52. 4/ 
DATA  AL(3,1,3,6) ,AL(3,1,3,7) ,AL(3,l,3,8)/6.75, .288, .0222/ 

DATA  AL(3, 1,4,1) ,AL(3, 1,4, 2)710400000.0,444000.0/ 

DATA  AL(3,1,4,3) ,AL(3,1,4,4) ,AL(3, 1,4, 5)776900.0,10800. 0,887. 
10/ 

DATA  AL(3, 1, 4, 6), AL(3, 1,4, 7) ,AL(3, 1,4, 8)780.6,2.06, .103/ 

DATA  AL(3,2,l,l),AL(3,2,l,2)/34.0,9.92/ 

DATA  AL(3,2,1,3) ,A1(3,2,1,4) ,AL(3,2,1,5)/4.88,2.16, .739/ 

DATA  AL(3,2,1,6) ,AL(3,2,1,7) ,AL(3,2,l,8)/.255, .0458, .00981/ 
DATA  AL(3,2,2,1) ,AL(3, 2, 2, 2)74290.0,608.0/ 

DATA  AL(3,2,2,3),AL(3,2,2,4),AL(3,2,2,5)/201.0,57.4,11.5/ 

DATA  AL(3,2,2,6) ,AL(3,2,2,7) ,AL(3,2,2,8)/2.47, .236, .0335/ 

DATA  AL(3, 2, 3, 1),AL(3, 2, 3, 2)7667000. 0,45300.0/ 

DATA  AL(3,2,3,3),AL(3,2,3,4),AL(3,2,3,5)/9960.0,1810.0,208.0/ 
DATA  AL(3,2,3,6) ,AL(3,2,3,7) ,AL(3, 2,3,8)726.7,1. 26, .109/ 

DATA  AL(3,2,4,1) ,AL(3, 2, 4, 2)761100000.0, 2330000.0/ 

DATA  AL(3,2,4,3) ,AL(3,2,4,4) , AL(3, 2, 4, 5)/374000. 0,47700. 0,350 

10.0/ 

DATA  AL(3,2,4,6) ,AL(3,2,4,7) ,AL(3,2,4,8)/292.0,7.16, .372/ 

END 

C  VARIABLES 

C  AE=LEVEL  DENSITY  ENERGY  PARAMETER 

C  AL(IL,IV)=INTERNAL  CONVERSION  COEFFICIENT 

C  AW=ATOMIC  WEIGHT 

C  CJ,CP,CE=CURRENT  LEVEL  SPIN,  PARITY,  AND  ENERGY 

C  DP(IG)=FOIL  ESCAPE  PROBABILITY  FOR  GAMMA  OF  GROUP  IG 
C  El (NS)  =  LOWEST  TABULATED  ENERGY  OF  SUBSHELL  TO  BE 
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C  INTERP.  EICEPT  FOR  L  AND  M  SHELL  WHICH  HAVE  El  (NS)  =  -1.0 
C  E2(NS)  =  LOWEST  VALUE  OF  THE  NS  SUBSHELL 
C  ED=ENERGY  DEPOSITED  IN  TANK  BY  CASCADE 

C  EF=AVERAGE  ENERGY  DEPOSITED  IN  DETECTOR  BY  CASCADES 
C  EO(IC,II)=ENERGY  OF  II-TH  GAMMA  IN  IC-TH  CASCADE 
C  ET=TRANSmON  ENERGY 

C  EU=MAI  ENERGY  DEPOSITED  IN  TANK  BY  GAMMA 
C  EG(IG)=UPPER  ENERGY  OF  GROUP  IG 

C  GN(II.IG)=GROUP  NUMBER  OF  II-TE  GAMMA  IN  CASCADE 
C  HF(IL,IP)=HINDRANCE  FACTOR  FOR  IL.IP 

C  HG(IL,IP)=HINDRANCE  FACTOR  FOR  IL.IP  TRANSITIONS 

C  IC=CASCADE  INDEX 

C  IE=E1CITATI0K  ENERGY  INDEX 

C  IJsFINAL  LEVEL  SPIN  PLUS  1/2 

C  IL=L  FOR  GAMMA  TRANSITION 

C  IP*1  FOR  ELECTRIC  OR  2  FOR  MAGNETIC  TRANSITION 
C  IT=TRANSITION  INDEX  WITHIN  CASCADE 

C  NC=NUMBER  OF  CASCADES 

C  NE=NUMBER  OF  INTERVALS  IN  ENERGY  CONTINUUM 

C  NG=NUMBER  TALLY  ENERGY  GROUPS 

C  NJ=LARGEST  LEVEL  SPIN  USED  IN  CONTINUUM  REGION 

C  NL=NUMBER  OF  MULTIPOLES  USED  IN  CALCULATION 

C  NP  =  THE  NUMBER  OF  INTERNAL  CONVERSION  COEFFICIENTS  TO 
C  INTERPOLATED.  THIS  WILL  BE  1,  FOR  THIS  PROGRAM. 

C  NR=INTEGER  PSEUDORANDOM  NUMBER 

C  NV(NS)=NUMBER  OF  VALUES  IN  TABLE  FOR  SHELL  NS 

C  NVD=NUMBER  OF  DISCRETE  LEVELS 

C  PC(IL,IP,IJ,IE)=INTERNAL  CONVERSION  PROBABILITY  IN 


C 


CONTINUUM 
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C  PE(IL,IP,IV)=INTERNAL  CONVERSION  COEFFICIENT  IN 

C  DISCRETE  LEVELS 

C  PF(IU)=PROBABILrn  GAMMA  DEPOSITS  ENERGY  IN  TANK  IN  INTERVAL  lU 
C  PS(ID)=PROBABILITY  PULSE  IS  IN  INTERVAL  ID 

C  PT(IL.IP.IJ.IE)=TRANSITION  PROBABILITY  FOR  CONTINUUM 

C  PV(IL,IP,IV)=TRANSITION  PROBABILITY  TO  DISCRETE  LEVEL  IV 
C  RL(IL)=TRANS1TI0N  FACTOR  FOR  2*IL  POLE 

C  RJ=LJ-.5 

C  RN=PSEUDORANDOM  NUMBER  ON  (0.,1.) 

C  RP(l)=WEISSKOPF  FACTOR  FOR  ELECTRIC  MULTIPOLE  RADIATION 

C  RP(2)=VfEISSKOPF  FACTOR  FOR  MAGNETIC  MULTIPOLE  RADIATION 

C  SG(IG)=GAMMA  MULTIPLICITY  IN  GROUP  IG 

C  SL(IL)=COMBINATORIAL  FACTOR  FOR  2.*IL-P0LE 

C  TE=T0TAL  ENERGY  RELEASE  PER  CASCADE 

C  TET=SUM  OF  ENERGIES  IN  A  CASCADE 

C  SJ=LEVEL  DENSITY  PARAMETER  FOR  SPIN 

C  TT(IL,IP,IV)=TOTAL  PROBABILITY  OF  TRANSITION  IN 

C  DISCRETE  LEVEL 

C  TU(IL,IP,IJ.IE')=TOTAL  TRANSITION  PROBABILITY  IN  THE 

C  CONTINUUM 

C  Tl,T2,T3,Il,I2=TEMPORARIES 

C  TJ=ONE  LINE  JOB  TITLE 

C  TL(IC,II)=L  VALUE  FOR  THE  II-TB  GAMMA  IN  IC-TH  CASCADE 
C  TP(IC,II)=IP  VALUE  FOR  THE  II-TH  GAMMA  IN  IC-TH  CASCADE 
C  TY=TOTAL  GAMMA  YIELD  (MULTIPLICITY)  PER  CASCADE 

C  VJ(IV),VP(IV),VE(IV)=SPIN.  PARITY,  ENERGY  OF  DISCRETE  LEVEL  IV 
C  IE,XJ,IP=INITIAL  NUCLEUS  ENERGY.  SPIN,  PARITY 


APPENDIX  B 

SAMPLE  OF  REPRESENTATIVE  CASCADES 


This  appendix  contains  a  listing  of  100  cascades. 
The  gamma  rays  in  each  cascade  are  listed  in  order 
of  emission  and  in  units  of  MeV. 


CASCADE  1 
0.30631E+01 
0.49367E+00 
0 . 59800E+00 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  5 
0 . 82840E+00 
0.14972E+01 
0.78662E+00 
0.16478E+01 
O.OOOOOE+00 

CASCADE  9 
0 . 98252E+00 
0.10797E+01 
0.62081E+00 
0.40736E+00 
0.39146E+00 


CASCADE  2 
0.27934E+01 
0 . 73640E+00 
0.11402E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  6 
0.15605E+01 
0.19830E+01 
0.83426E+00 
0.17300E+00 
0 . OOOOOE+00 

CASCADE  10 
0.15605E+01 
0.12265E+01 
0.12087E+01 
0 . OOOOOE+00 
O.OOOOOE+00 


CASCADE  3 
0 . 19843E+01 
0.98107E+00 
0.70567E+00 
0.91196E+00 
0 . OOOOOE+00 

CASCADE  7 
0.13678E+01 
0 . 85739E+00 
0.15383E+01 
0 . 86653E+00 
0 . OOOOOE+00 

CASCADE  11 
0.13293E+01 
0.18045E+01 
0.97326E+00 
0.56300E+00 
0 . OOOOOE+00 


CASCADE  4 
0.20228E+01 
0 . 14184E+01 
0.92989E+00 
O.OOOOOE+00 
0 . OOOOOE+OO 

CASCADE  8 
0.17146E+01 
0.17642E+01 
0.62822E+00 
0 . OOOOOE+OO 
O.OOOCOE+00 

CASCADE  12 
0.75134E+00 
0.12872E+01 
0.17328E+01 
0.17300E+00 
0. OOOOOE+OO 
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CASCADE  13 
0.24081E+01 
0 . 88860E+00 
0.15063E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  17 
0 . 19075E+01 
0 . 74634E+00 
0.13978E+01 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  21 
0 . 14449E+01 
0.22757E+01 
0 . 94945E+00 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  25 
0.15605E+01 
0.14328E+01 
0.79522E+00 
0.10145E+01 
O.OOOOOE+00 

CASCADE  29 


CASCADE  14 
0.26778E+01 
0.10988E+01 
0.10264E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  18 
0 . 90546E+00 
0.15475E+01 
0.62304E+00 
0.71877E+00 
0 . 86428E+00 

CASCADE  22 
0.94399E+00 
0.11200E+01 
0.61722E+00 
0.55662E+00 
0  35374E+00 

CASCADE  26 
0.10596E+01 
0.17459E+01 
0 . 65830E+00 
0.99634E+00 
0 . OOOOOE+00 

CASCADE  30 


CASCADE  15 
0.36604E+00 
0.34347E+01 
0 . 86930E+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  19 
0.22155E+01 
0 . 12527E+01 
0 . 13348E+01 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  23 
0.82840E+00 
0.22231E+01 
0 . 16075E+01 
0 . 14400E+00 
0 . OOOOOE+00 

CASCADE  27 
0 . 23696E+01 
0 . 10755E+01 
0.81575E+00 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  31 


CASCADE  16 
0.71281E+00 
0.19312E+01 
0.87651E+00 
0.10915E+01 
O.OOOOOE+00 

CASCADE  20 
0.26008E+01 
0.65742E+00 
0.14008E+01 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  24 
0.19072E+01 
0.10410E+01 
0 . 10693E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  28 
0.12137E+01 
0.24677E+01 
0.97755E+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  32 


0.24852E+01 
0.97799E+00 
0.11958E+01 
0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  33 
0.18302E+01 
0.11631E+01 
0.25361E+00 
0.15131E+01 
O.OOOOOE+00 

CASCADE  37 
0.22155E+01 
0.12200E+01 
0.11946E+01 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  41 
0.19458E+01 
0.16688E+01 
0.11205E+01 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  45 
0.13678E+01 


0.15990E+01 
0.16567E+01 
0 . 53459E+00 
0.10127E+01 
0 . OOOOOE+00 

CASCADE  34 
0.24081E+01 
0.72966E+00 
0.95021E+00 
0.58200E+00 
O.OOOOOE+00 

CASCADE  38 
0.55869E+00 
0.28825E+01 
0.13'')18E+01 

0.0000.';  +00 

O.OOOOOE+00 

CASCADE  42 
0.17531E+01 
0.19844E+01 
0 . 87449E+00 
0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  46 
0.33328E+01 


0.20614E+01 
0 . 16214E+01 
0.11202E+01 
0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  35 
0.11366E+01 
0 . 13446E+01 
0 . 63787E+00 
0 . 10008E+01 
0 . OOOOOE+00 

CASCADE  39 
0.23311E+01 
0.96643E+00 
0.16411E+00 
0 . 96688E+00 
0 . OOOOOE+00 

CASCADE  43 
0.22925E+01 
0 . 14746E+01 
0 . 10358E+01 
O.OOOOOE+00 
0 . OOOOOE+OO 

CASCADE  47 
0.20999E+01 


0.29090E+01 
0 . 17030E+01 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 

CASCADE  36 
0.32173E+01 
0 . 15857E+01 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 

CASCADE  40 
0.12522E+01 
0.11313E+01 
0.78615E+00 
0 . 65372E+00 
0.56300E+00 

CASCADE  44 
0.23696E+01 
0.89746E+00 
0 . 54200E+00 
0.86095E+00 
0 . OOOOOE+OO 

CASCADE  48 
0.31787E+01 


0.73313E+00 

0.67454E+00 

0.64112E+00 

0.13864E+01 

CASCADE  49 
0.25622E+01 
0.12198E+01 
0.88799E+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  53 
0.35640E+01 
0.10950E+01 
O.OOOOOE+00 
0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  57 
0.36796E+01 
0.97938E+00 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  61 
0.10596E+01 
0.11313E+01 


0.10917E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  50 
0.15990E+01 
0.20624E+01 
0 . 99759E+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  54 
0.28705E+01 
0 . 77396E+00 
0.19100E+00 
0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  58 
0.10596E+01 
0.13548E+01 
0.14026E+01 
0.84197E+00 
0 . OOOOOE+00 

CASCADE  62 
0.15605E+01 
0.15245E+01 


0.16742E+01 
0 . 8D5fts’^+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  51 
0 . 20999E+01 
0 . 15866E+01 
0.98355E+00 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  55 
0.33714E+01 
0 . 14316E+01 
O.OOOOOE+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  59 
0.23696E+01 
0 . 77068E+00 
0.19100E+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  63 
0 . 24852E+01 
0.12379E+01 


0.14913E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  52 
0.82840E+00 
0.17391E+01 
0 . 70057E+00 
0.35970E+00 
0.88143E+00 

CASCADE  56 
0.10210E+01 
0.14868E+01 
0.10291E+01 
0.10751E+01 
0 . 14800E+00 

CASCADE  60 
0.20228E+01 
0.16929E+01 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  64 
0.20228E+01 
0.51561E+00 


0.15374E+01 
0 . 10747E+01 
O.OOOOOE+00 

CASCADE  65 
0.29090E+01 
0.33512E+00 
0.22000E+00 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  69 
0.11752E+01 
0.24502E+01 
0.80566E+00 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  73 
0.18302E+01 
0.47536E+00 
0.14778E+01 
0.32364E+00 
0 . OOOOOE+00 

CASCADE  77 
0.94399E+00 
0.13818E+01 
0.14585E+01 


0.39972E+00 

0.77784E+00 

0.19100E+00 

CASCADE  66 
0.20999E+01 
0.10080E+01 
0.71154E+00 
0.62100E+00 
0 . OOOOOE+00 

CASCADE  70 
0.63575E+00 
0.22038E+01 
0.88676E+00 
0.91650E+00 
O.OOOOOE+00 

CASCADE  74 
0.19843E+01 
0.96238E+00 
0.11449E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  78 
0 . 12522E+01 
0.87126E+00 
0 . 95988E+00 


0.93594E+00 

O.OOOOOE+00 

O.OOOOOE+00 

CASCADE  67 
0 . 16375E+01 
0 . 16284E+01 
0.95571E+00 
0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  71 
0 . 21384E+01 
0.67726E+00 
0.95952E+00 
0.88380E+00 
0 . OOOOOE+00 

CASCADE  75 
0.22155E+01 
0.10398E+01 
0 . 53494E+00 
0.82176E+00 
O.OOOOOE+00 

CASCADE  79 
0.75134E+00 
0 . 24658E+01 
0.15858E+01 


O.OOOOOE+00 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  68 
0.26393E+01 
0.96489E+00 
0 . 10658E+01 
0 . OOOOOE+00 
0 . OOOOOE+OO 

CASCADE  72 
0.28320E+01 
0 . 87299E+00 
0.95405E*00 
0. OOOOOE+OO 
0. OOOOOE+OO 

CASCADE  76 
0.11366E+01 
0.14533E+01 
0.58735E+00 
0.54399E+00 
0.94877E+00 

CASCADE  80 
0.32943E+01 
0.13177E+01 
0 . OOOOOE+OO 


0 . 87472E+00 
O.OOOOOE+00 

CASCADE  81 
0 . 12522E+01 
0.20416E+01 
0.56300E+00 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  85 
0 . 12908E+01 
0.22420E+01 
0 . 12273E+01 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  89 
0.20614E+01 
0.14781E+01 
0 . 12635E+01 
O.OOOOOE+00 
O.OOOOOE+00 

CASCAD':  93 
0.10520t.  01 
0.88750E+00 
0.10626E+01 
0.10881E+01 


0 . 66573E+00 
0.90990E+00 

CASCADE  82 
0 . 13293E+01 
0.19306E+01 
0.32916E+00 
0.10229E+01 
0 . OOOOOE+00 

CASCADE  86 
0.21769E+01 
0.16174E+01 
0.81766E+00 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  90 
0.16375E+01 
0.18278E+01 
0.11937E+01 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  94 
0.13678E+01 
0.19260E+01 
0.12392E+01 
0 . OOOOOE+00 


0 . OOOOOE+00 
O.OOOOOE+00 

CASCADE  83 
0.27549E+01 
0.76318E+00 
0 . 11409E+01 
O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  87 
0.24852E+01 
0.70442E+00 
0.62027E+00 
0 . 99312E+00 
O.OOOOOE+00 

CASCADE  91 
0.24852E+01 
0.12516E+01 
0 . 10663E+01 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  95 
0 . 24467E+01 
0 . 72427E+00 
0.16321E+01 
0 . OOOOOE+00 


O.OOOOOE+00 
0 . OOOOOE+00 

CASCADE  84 
0.21384E+01 
0.14145E+01 
0.11061E+01 
0 . OOOOOE+00 
0 . OOOOOE+00 

CASCADE  88 
0 . 14449E+01 
0.22757E+01 
0.89145E+00 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  92 
0.21769E+01 
0.82965E+00 
0.16054E-+01 
O.OOOOOE+00 
O.OOOOOE+00 

CASCADE  96 
0.12137E+01 
0.16232E+01 
0.93992E+00 
0.83521E+00 


O.OOOOOE+00 

CASCADE  97 
0.14449E+01 
0.49367E+00 
0 . 83279E+00 
0.38399E+00 
0.14567E+01 


0 . OOOOOE+00 

CASCADE  98 
0.22925E+01 
0.77243E+00 
0 . 46100E+00 
0.10860E+01 
0 . OOOOOE+00 


0 . OOOOOE+00 

CASCADE  99 
0 . 16375E+01 
0 . 15841E+01 
0 . 33149E+00 
0.10769E+01 
0 .  OOOOOE+00 


0 . OOOOOE+00 

CASCADE  100 
0.90546E+00 
0.16654E+01 
0 . 83983E+00 
0.42687E+00 
0 . 96548E+00 


